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Entwicklung neuartiger Mess- und Auswertungsstrategien für 
ein scannendes Wolkenradar und deren Anwendungsbereiche 



Kurzbeschreibung 

Die Beobachtung wolkenphysikalischer Prozesse wird erst durch räumlich und zeitlich hoch 
aufgelöste Messungen mit Wolkenradaren möglich. Nachdem vertikal zeigende Geräte weit 
verbreitet sind, erweitem scannende Wolkenradare die Messmöglichkeiten, denn sie erlauben 
es, einen bestimmten Bereich der Wolke mehrfach nacheinander zu erfassen. Nur so kann die 
zeitliche Entwicklung innerhalb der Wolke getrennt von der räumlichen Variabilität erfasst wer- 
den. 

In der vorliegenden Arbeit werden Messungen mit dem scannenden 35.5 GHz Wolkenradar 
MIRA36-S vorgestellt. Der Schwerpunkt liegt dabei auf der Entwicklung geeigneter Scan- und 
Auswerteverfahren, die von der neuen technischen Möglichkeit Gebrauch machen. Bei den mei- 
sten hier vorgestellten Messungen handelt es sich um eine schnelle Abfolge von RHl-Scans, die 
in Richtung des vorherrschenden Windes durchgeführt wurden. Die Auswerteverfahren sind 
so konzipiert, dass sie auch automatisiert werden können, also für den operationeilen Einsatz 
geeignet sind. 

Ein Algorithmus zur Detektion und Analyse der Schmelzschicht wird vorgestellt und disku- 
tiert. Neben der mittleren Höhe der Schmelzschicht können ihre Ober- und Unterkante identi- 
fiziert werden. Durch den Einsatz von drei verschiedenen Messgrößen (Reflektivität, LDR und 
vertikale Streuergeschwindigkeit) lassen sich verschiedene Teilprozesse in der Schmelzschicht 
räumlich getrennt identifizieren. Unter anderem kann gezeigt werden, dass die Schmelzdauer 
ab einer minimalen Schmelzschichtdicke nahezu unabhängig von den Eigenschaften der Hy- 
drometeore direkt oberhalb der Schmelzschicht bei etwa 150 s liegt. 

Ein weiterer Algorithmus zerlegt die gemessenen radialen Geschwindigkeiten in großräumige 
Anteile und einen Restanteil, der als residuale Geschwindigkeit definiert wird. Die in dieser 
Arbeit eingeführte residuale Geschwindigkeit ist ein beachtliches Werkzeug, das interne Fluk- 
tuationen aufzeigt, mit denen zusätzliche Erkenntnisse über die interne Dynamik der Wolke 
ermöglicht werden. Mit Hilfe der residualen Geschwindigkeiten sind zudem Wellen (z.B. in- 
terne Schwerewellen oder Kelvin-Helmholtz-Instabilitäten) zu erkennen und zu analysieren. 
Eigenschaften wie Amplitude, Wellenlänge und Phasengeschwindigkeit lassen sich daraus ab- 
leiten. 



Development of new measurement and processing strategies 
for a scanning cloud radar and their applications 



Abstract 

The Observation of cloud physical processes is only possible tiirough spatially and temporally 
resolved measurements witii cloud radars. Although vertically pointing devices are widespread, 
scanning cloud radars expand the measurement capabilities, because tiiey allow observations of 
a certain area of the cloud several times in succession. This aUows to record the development 
within the cloud over time separately from the spatial variability. 

In this werk measurements are presented with the scanning 35.5 GHz cloud radar MIRA36- 
S. The focus is on the development of appropriate scanning and evaluation methods that avail 
the new technical possibility. Most measurements presented here are a rapid succession of RHl 
Scans, performed in the direction of the prevailing wind. The evaluation methods are designed 
to be automated for operational use. 

An algorithm for detection and analysis of the melting layer is presented and discussed. In 
addition to the mean height of the melting layer, its upper and lower limits are identified. By 
using three different measures (Reflectivity, LDR and vertical scatterer velocity) various sub- 
processes in the melting layer can be identified in spatially Separation. Furthermore, it can be 
shown that the melting time from a minimum melting layer thickness is nearly independent of 
the properties of the hydrometeors directly above the melting layer, and amounts to approxima- 
tely 150 s. 

Another algorithm decomposes the measured radial velocities in large-scale portions and a re- 
maining portion, which is defined as the residual velocity. The residual velocity introduced in 
this publication is a remarkable tool revealing internal fluctuations which increase the knowled- 
ge of the internal dynamics of the cloud. By using the residual velocities, waves (e.g. internal 
gravity waves or Kelvin-Helmholtz instabilities) can also be detected and analyzed. Characteri- 
stics such as amplitude, wavelength and phase velocity can be derived from it. 
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1 Einleitung 

Wolken spielen eine entscheidende Rolle im Strahlungs- und Wasserhaushalt der Erde. Sie ent- 
stehen durch Kondensation oder Deposition von Wasser, meist an Aerosolen und absorbieren, 
reflektieren und emittieren Strahlung. Wolken bestehen aus Eispartikeln und Tropfen, so ge- 
nannten Hydrometeoren. Durch Kondensation, bzw. Deposition und Stoßprozesse nehmen die- 
se an Größe zu und beginnen zu faUen. Niederschlag in Form von Regen, Schnee, Graupel, 
Hagel ist die Folge. Welche thermodynamischen und dynamischen Prozesse in welchem Ma- 
ße die Bildung von Niederschlag beeinflussen, ist bisher noch nicht ausreichend verstanden und 
quantifiziert worden und Gegenstand der Forschung. In dieser Arbeit soll mit Hilfe von Messun- 
gen mit einem Wolkenradar das Verständnis von dynamischen Vorgängen in Wolken erweitert 
werden. 

Radargeräte werden seit 1941 zur meteorologischen Femerkundung eingesetzt. Damit werden 
erstmals flächige Messungen von Niederschlags- und Wolkenbereichen ermöglicht. So wurde 
bereits 1941 von J.W. Ryde bei General Electric das erste Wetterradar entwickelt und genutzt. 
Das erste Radar kurzer Wellenlänge (1.25 cm), das auch Wolkenpartikel detektieren konnte, 
wurde in den 50er Jahren von Plank et al. (1955) vorgestellt. Die kurze Wellenlänge eines 
Wolkenradars ermöglicht die Detektion kleinster Wassertropfen und Eispartikel, die ein Radar 
größerer Wellenlänge (z.B. Niederschlagsradare) nicht erfasst. Wolkenradare helfen ebenfalls 
bei der Abschätzung von mikrophysikalischen und dynamischen Parametern in Wolken, wie 
beispielsweise der Größenverteilung der Hydrometeore und von Turbulenzparametern. Um die- 
se Parameter abzuleiten, wurden meist Messungen vertikal ausgerichteter Wolkenradare ver- 
wendet. Die mehrdimensionale Dynamik von Strukturen in Wolken in Kombination mit der 
zeitlichen Entwicklung dieser Strukturen wurde bisher kaum untersucht und ist Gegenstand 
dieser Arbeit. Auf der Basis von schnellen zweidimensionalen Messungen, so genannten RHI- 
Scans \ werden neuartige Verfahren entwickelt, gezeigt und angewendet: 

• Eine automatisierte Detektion der Schmelzschicht, um Bereiche zu separieren, in denen 
die Hydrometeore weitgehend in flüssiger, fester oder gemischter Phase vorliegen. 

• Die Visualisierung von Auf- und Abwindbereichen, um dynamische Transportvorgänge 
detektieren, unterscheiden und verfolgen zu können. 

Radare mit Wellenlängen im Millimeterbereich werden heutzutage weltweit eingesetzt. Diese 
sind in der Regel vertikal ausgerichtet und meist nicht mobil. Beispielsweise im ARM^ Pro- 
gramm, gefördert durch das amerikanische Ministerium für Energie (DOE^) werden mehrere 

'Range-Height-Indicator-Scans, die Radarantenne wird nur in der Elevation bewegt, siehe Kap. 3.2 
^Atmospheric Radiation Measurement 
^United States Department of Energy 



6 



1 Einleitung 



stationäre Anlagen in Alaska, den Great Plains und im tropischen Westpazifik sowie eine mo- 
bile Anlage betrieben. Damit werden kontinuierliche Messungen von Wolken und Strahlung 
durchgeführt (Stokes und Schwartz, 1994; Ackerman und Stokes, 2003). Zum Einsatz kommen 
mehrere vertikal zeigende, polarimetrische 35-GHz-Doppler-Radare, zwei vertikal zeigende po- 
larimetrische 95-GHz-Doppler-Radare und ein scannendes polarimetrisches 95-GHz-Doppler- 
Radar. Diese haben das Ziel, die Wolkenobergrenzen, die Reflektivitäten, die vertikalen (bzw. 
die radialen) Partikelgeschwindigkeiten, sowie die Doppler-Streubreiten in Wolken zu messen. 
Dabei liegt der Schwerpunkt auf der kontinuierlichen Erfassung von Wolkenstrukturen in Ver- 
bindung mit Strahlungsmessungen und der Bereitstellung von Datenprodukten, die zur Weiter- 
entwicklung von Klima-Modellen beitragen. 

CloudNet, ein Forschungsprojekt, das von der europäischen Kommission gefördert wurde, be- 
treibt ein Netzwerk von drei Wolkenfemerkundungsstationen in den Niederlanden, Frankreich 
und England (Illingworth et al., 2007). Dabei werden vertikal ausgerichtete 95-GHz-Doppler- 
Radare und polarimetrische 35-GHz-Doppler-Radare als Teil ihrer Messgeräteausstattung ver- 
wendet, um die Repräsentation von Wolken in europäischen Wettervorhersagemodellen zu eva- 
luieren. Diese bodengestützten Wolkenradare messen die vertikale Verteilung von Wolken- 
schichten und ermöglichen die Anfertigung von Statistiken, beispielsweise des Wolkenbede- 
ckungsgrades mit der zugehörigen Höhe der Wolken, der Höhen der Wolkengrenzen und der 
Radar-Reflektivitäten für verschiedene Wolkentypen. Wenn ein Radarsystem an einem festen 
Standort über viele Jahre kontinuierlich betrieben wird, decken die Daten den vollständigen Be- 
reich der dynamischen und mikrophysikalischen Systeme ab, die für diesen Einsatzort typisch 
sind. Diese Datensätze können auch genutzt werden, um Informationen der Wolkenschichtung 
in den Kontext von großskaligen oder mesoskaligen dynamischen Systemen zu bringen. Solche 
Informationen können dann verwendet werden, um Wechselwirkungen zwischen dynamischen 
und mikrophysikalischen Prozessen zu studieren und zur Verbesserung der Fähigkeit von Mo- 
dellen beitragen, diese Interaktionen zu simulieren. 

Wie auch in den gerade beschriebenen Fällen, wurden meist rein vertikal zeigende Wolkenra- 
dare eingesetzt. Damit kann nur eine Zeitreihe am Ort der Messung dokumentiert werden. Die 
räumliche Änderung der Wolkenstruktur ist damit nicht beobachtbar. Abhilfe schaffen scan- 
nende Wolkenradare. Das erste wissenschaftlich genutzte, scannende Wolkenradar wurde in 
den frühen 80er Jahren beim Environmental Technology Laboratory (ETL) der amerikanischen 
Wetter- und Ozeanografiebehörde (NOAA) entwickelt. Dieses Radar ist noch immer im Einsatz. 

Nachdem Wolken in der Klimaforschung zunehmend an Bedeutung gewannen (siehe z.B. Ste- 
phens und Webster, 1981), wurden in den 1990er Jahren dank zahlreicher ingenieurtechnischer 
Neuerungen neu konstruierte bzw. umfangreich erneuerte Wolkenradare in verschiedenen For- 
schungsprojekten eingesetzt (z.B. Pazmany et al., 1994; Clothiaux et al., 1995). Die Geschichte 
der angewandten Radarmeteorologie am Institut für Meteorologie und Klimaforschung (IMK) 
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des KIT' beginnt 1993 mit der Inbetriebnahme eines C-Band-Radar. Aus dessen Daten wird 
beispielsweise die flächendeckende Niederschlagsverteilung im Umkreis von 120 km abgelei- 
tet (Hannesen, 1998; Dotzek, 1999; Peters, 2008). Die Messungen mit einem C-Band-Radar 
sind darauf beschränkt, große Niederschlagspartikel zu delektierten. Wolken, die aus kleinen 
Hydrometeoren bestehen, sind mit diesem Radar nicht messbar. Um auch die Vorgeschichte 
der Niederschlagsentwicklung, also die Frühphase von niederschlagsbildenden Wolken erfas- 
sen zu können, wurde im Jahr 2005 ein scannendes Wolkenradar in Betrieb genommen, dessen 
Messungen die Grundlage für diese Arbeit sind. Dieses Wolkenradar ist ein mobiles 35.5-GHz- 
Wolkenradar mit der Bezeichnung MIRA36-S und besitzt eine bewegliche Antenne. Es ist ei- 
ne Weiterentwicklung des rein vertikal abtastenden Wolkenradars MIRA36. MIRA36-S wurde 
auf Basis von Vorgaben des meteorologischen Instituts des KIT von der Herstellerfirma Me- 
tek GmbH entwickelt und besitzt eine weltweit selten realisierte Scan-Einheit, die neben dem 
gesamten Azimutbereich auch einen Elevationsbereich von 45 ° bis zum Zenit überdeckt. 

MIRA36-S ist seit 2005 im Einsatz. Nachdem erste Tests auf der Messwiese des KIT-Campus- 
Nord erfolgreich waren, wurden RHI-Messungen von MIRA36-S mit Vertikalmessungen von 
MIRA36 in Bezug gesetzt (Handwerker und Görsdorf, 2006). Das Resultat zeigte, dass Zeit- 
Höhen-Schnitte nur in manchen Situationen die räumliche Struktur von Wolken wiedergeben 
können. Die Notwendigkeit von abtastenden Wolkenradaren wurde damit unterstrichen. 

Der Vergleich von Messungen des Wolkenradars M1RA36-S (vertikal ausgerichtete Antenne) 
mit einem 95 -GHz- Wolkenradar (W-Band) des ARM-Programms lieferte eine qualitativ gute 
Übereinstimmung (Handwerkerund Miller, 2008). Das W-Band-Radarmaß Reflektivitäten, die 
im Mittel 3 dB über denen von MIRA36-S liegen. Oberhalb etwa -10 dBz werden Mie-Effekte 
sichtbar und das Ka-Band-Radar MIRA36-S zeigt mit zunehmender Reflektivität höhere Werte 
als das W-Band-Radar. Das Wolkenradar MIRA36-S deckt dabei einen, um ca. 15 dB größeren, 
Bereich der Reflektivität ab, von —65 dBz bis 20 dBz. Die Messungen der Vertikalgeschwin- 
digkeiten zeigen eine gute Übereinstimmung. Fallen die Streuer schneller als mit 3 m/s, so sind 
die Werte der gemessenen Vertikalgeschwindigkeiten des W-Band-Radars niedriger. Aufgrund 
der größer werdenden Hydrometeore und der niedrigeren Wellenlänge kommt zunehmend die 
Mie-Streuung zum Tragen. Weiterhin ist MIRA36-S sensibler gegenüber nichtmeteorologi- 
schen Streuern, wie beispielsweise Insekten und Staub. Durch eine bessere Separation des co- 
polaren und cross-polaren Signals ist eine zuverlässigere Messung des Imearen Depolarisations- 
Verhältnisses (LDR, Kap. 3.4.2) gewährleistet. Dadurch können Wolken innerhalb der Grenz- 
schicht besser von nichtmeteorologischen Streuern unterschieden werden. 

Gemeinsame Messungen von MIRA36-S zusammen mit einem Lidar- Windpro filer (Wellenlänge 
A = 2023 nm) zeigten gute Übereinstimmungen für Profile des Horizontalwinds in gleichen 

'Karlsruher Institut für Technologie 



8 



1 Einleitung 



Messbereichen (Handwerker et al., 2008), die aus so genannten PPI-Scans^ abgeleitet und mit- 
tels VAD^-Algorithmus berechnet wurde. Unterschiede der gemessenen Vertikalgeschwindig- 
keiten beider Messgeräte konnten dazu genutzt werden, die Tropfengrößenverteilung im leich- 
ten Niederschlag abzuschätzen (Träumner et al., 2010). 

Damit wurde gezeigt, dass MIRA36-S, verglichen mit anderen Wolkenradaren, herausragen- 
de Messeigenschaften aufweist, die für zweidimensionale Messungen genutzt werden können 
(Grenzhäuser und Handwerker, 2010), welche Gegenstand dieser Arbeit sind. 



Plan-Position-Indicator-Scans, die Radarantenne wird nur im Azimut bewegt, siehe Kap. 3.2 
^Velocity Azimuth Display, siehe Kap. 2.4.2 
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2 Grundlagen der Radarmeteorologie 

Radare sind gegenwärtig in der Meteorologie wichtige Messgeräte der Femerkundung von Nie- 
derschlag und WoUcen. In erster Linie werden Radargeräte in der Meteorologie eingesetzt, 
um Niederschlag zu messen. Dabei dienen Messungen mittels Radar unter anderem zur Vi- 
sualisierung und Abschätzung der Entwicklung von Niederschlagsgebieten. Daraus lassen sich 
beispielsweise Klimatologien von Niederschlagsepisoden erstellen (Carbone et al., 2002). Die 
aus Messungen abgeleiteten Reflektivitäten und Dopplergeschwindigkeiten liefern wichtige Pa- 
rameter, die von Wetterdiensten und meteorologischen Warndiensten zunehmend zur Kurz- 
fristvorhersage verwendet werden. Auch für die Hydrologie sind Radarmessungen von großer 
Bedeutung, da Niederschlagsraten direkt abgeleitet (Beschreibungen verschiedener Methoden 
finden sich u.a. in Beheng (1998)) und akkumulierte Niederschlagsmengen zur Vorhersage 
von Abflussraten herangezogen werden können. Weiterhin werden Radarmessungen zur Ve- 
rifikation modellierter Niederschlagsparameter (z.B Pfeifer et al., 2008) und Wolkenparame- 
ter (z.B Illingworth et al., 2007) verwendet. Die Anwendungsbereiche von Radardaten sind 
somit sehr vielfältig und die gemessenen Daten werden sogar dazu genutzt, um beispielswei- 
se Erwärmungsraten während der Niederschlagsbildung abzuschätzen („latent heat nudging": 
Leuenberger und Rossa, 2007). 

Niederschlagsradargeräte werden überwiegend operationeil von Wetterdiensten betrieben. In 
vielen europäischen Ländern werden C-Band Radare eingesetzt. Der Deutsche Wetterdienst bei- 
spielsweise verwendet in seinem Radarverbund 16 C-Band Radare, die mit einer Wellenlänge 
von ca. 5 cm und einer Frequenz von 5.6 GHz eirbeiten. Der amerikanische und der französische 
Wetterdienst setzen auch S-Band Radare mit einer Wellenlänge von 10 cm (3 GHz) im opera- 
tionellen Betrieb ein. Im Vergleich zu Messungen von Niederschlag werden Messungen von 
Wolken mit Radargeräten geringer Wellenlänge, so genannten WoUcenradaren, durchgeführt. 
Die typischen Frequenzen von Wolkenradaren sind 35 GHz (0.9cm) und 95 GHz (0.3 cm). Im 
Übergangsbereich zwischen Wolken- und Niederschlagsradaren werden X-Band-Radare (3 cm, 
10 GHz) genutzt, etwa für mobile Einsätze. Die Wahl der Wellenlänge und anderer technischer 
Spezifikationen (z.B. Pulswiederholfrequenz) ist immer ein Kompromiss zwischen der Detekti- 
on von Streuern in einem bestimmten Größenbereich, der Dämpfung durch Streuer, der Reich- 
weite, sowie der zeitlichen und räumlichen Auflösung der Messung. 

Das allgemeine Radarprinzip besteht im Aussenden von Mikrowellenstrahlung und dem an- 
schließenden Empfangen der rückgestreuten Strahlung. Für meteorologische Radare werden 

meist kurze Pulse elektromagnetischer Wellen mit Pulslängen im Bereich einer Mikrosekunde 
und bei einem Strahlöffnungswinkel im Bereich von einem Grad gesendet. Aus der Laufzeit der 
empfangenen Welle und der Antennenposition wird der Ort des empfangenen Signals bestimmt. 
Die empfangene Leistung (proportional zum Quadrat der AmpHtude der empfangenen Welle) 
hängt von der Entfernung der Streuer, von radarspezifischen Größen, wie der Antennenform, 
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der Niederschlagsdämpfung und vom Radarreflektivitätsfaktor ab. Der Radarreflektivitätsfaktor 
beschreibt U.A. die Rückstreueigenschaften von Streuern im Strahlvolumen. Die Radialkompo- 
nente des Geschwindigkeitsvektors kann bei Doppler-fähigen Radaren aus der Veränderung 
der Phasendifferenz zwischen abgestrahlter und empfangener Welle ermittelt werden. Wird 
zusätzlich die Änderung der Polarisation von gesendeter und empfangener Strahlung gemes- 
sen, lassen sich daraus Schlüsse über die Form (Kugelform, asphärische Partikel) und die Art 
(flüssig, eisförmig, schmelzend) ableiten. 

In den nachfolgenden Unterkapiteln folgen einige allgemeine Definitionen und Zusammenhän- 
ge, die für viele Radartypen gelten (siehe dazu auch Sauvageot, 1992; Doviac und Zrnic, 1984; 
Skolnik, 1 990), bevor in Kapitel 3 auf das Wolkenradar MIR A36-S, mit dem die in dieser Arbeit 
analysierten Radarmessungen durchgeführt wurden, eingegangen wird. Dabei wird zuerst die 
Streuung und Dämpfung der elektromagnetischen Wellen in der Atmosphäre beschrieben, bevor 
die Radargleichung vorgestellt und die Geschwindigkeit von Streuelementen besprochen wird. 

2.1 Streuung 

Breitet sich eine elektromagnetische Welle in einem Medium aus, so führen Änderungen des 
Brechungsindex zur Streuung der Welle. Für die Radarstrahlung wichtig ist die elastische Streu- 
ung (z.B. Rayleigh-Streuung), die sich dadurch auszeichnet, dass die Summe der kinetischen 
Energien nach der Wechselwirkung mit Streukörpern genauso groß ist wie zuvor. Dabei werden 
elektrische und magnetische Multipolschwingungcn im Streuer angeregt, die ihrerseits wie- 
der als Sender wirken. Streu-Elemente in der Atmosphäre sind hauptsächlich sphärische und 
abgeplattete Wassertropfen sowie sphärische und teilweise unregelmäßig geformte Eispartikel. 
Neben der Streuung an Hydrometeoren gibt es auch Streuung an nicht meteorologischen Objek- 
ten. Abhängig vom Radarstandort und von den meteorologischen Bedingungen, kann es auch 
zur Streuung an unbewegten bodennahen Objekten, so genanntem „ground clutter", sowie an 
bewegten Objekten, beispielsweise an Vögeln, Insekten und Staubpartikeln, kommen. 

Ausgehend von den Maxwell-Gleichungen, die die Ausbreitung von elektromagnetischen Wel- 
len allgemein beschreiben, hat Mie (1908) eine vollständige Theorie der Streuung von ebe- 
nen Wellen an einer homogenen dielektrischen Kugel in einem nicht absorbierenden Medium 
entwickelt. Dabei lässt sich die von einem Streu-Element emittierte Strahlung als eine Reihe 
von Partialwellen berechnen. Sind die Streuer mit Durchmesser D sehr klein, verglichen mit 
der Wellenlänge A, lassen sich die Streuteilchen als elektrische Dipole behandeln (Rayleigh- 
Näherung). Mit zunehmender Größe der Streuteilchen bei konstanter Wellenlänge werden elek- 
trische und magnetische Multipolmomente angeregt, und es treten Resonanzen auf. Dadurch 
gewinnt die Vorwärtsstreuung deutlich an Dominanz. An den Rayleigh-Bereich schließt sich 
somit der Mie-Bereich an, in welchem der Rückstreuquerschnitt als Funktion des Streuer- 
Durchmessers oszilliert. 
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Abbildung 2.1: Logarithmische Darstellung des normierten Rückstreuquerschnitts als Funkti- 
on des Verhältnisses vom Streuerumfang zur Radarwellenlänge für eine homogene Kugel, aus 
Skolnik (1990). Hierbei ist D = 2a, mit dem Kugelradius a 

Für die Radarmessung ist nur der Anteil der Strahlung von Bedeutung, der zurückgestreut wird. 
Der Streuquerschnitt ist dabei abhängig vom komplexen Brechungsindex m und von der radio- 
elektrischen Größe aR = Dtt/X = 2ti a/\, mit dem Kugelradius a. Der differentielle Streu- 
querschnitt für einen Winkel von 180 °, also gerade entgegengesetzt zur Strahlrichtung, ist der 
Rückstreuquerschnitt, der für den Fall der Rayleigh-Näherung gegeben ist durch 



(2.1) 



mit lifP = Dielektrizitätsfaktor. 



Die vollständige Mie-Theorie ist anhand der Abbildung 2.1 veranschaulicht: Der mit der Ku- 
gelschnittfläche normierte Rückstreuquerschnitt ist als Funktion der radioelektrischen Größe 
dargestellt. Werte von 2 tt a < 0.3 kennzeichnen den Rayleigh-Bereich (Steigung der Kur- 
ve ist ä; 4 bei der doppelt logarithmischen Auftragung). Ebenso erkennt man für > 0.3 den 
Resonanzbereich der Mie-Streuung. 
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Der Dielektrizitätsfaktor ist definiert als 

m-^ — 1 

(2.2) 



Der komplexe Brechungsindex m und somit \KY hängen bei einer gegebenen Substanz und 
deren Aggregatzustand von der Wellenlänge und der Temperatur des Streuers ab. So hat 
der Dielektrizitätsfaktor beispielsweise für die Wellenlänge des Wolkenradars MIRA36-S von 
A = 8.44 mm folgende Werte für Eis, für unterkühltes Wasser bei -8°C und für Wasser bei 
20 °C (Gunn und East, 1954): 



I^lice,,=0.92g/cm3 « O-l». 1^ Itater, T=-8 ° C ^ 0.84, | | tater, T=20 ° C ^ 0.91 (2.3) 

Wegen der Differenz der Dielektrizitätsfaktoren zwischen Wasser und Eis ist der Rückstreu- 
querschnitt <j für homogene, kugelförmige Hydrometeore des gleichen Durchmessers bis zu 5 
mal stärker für Wasserpartikel als für Eispartikel. Der Durchmesser eines Eishydrometeors muss 
gegenüber dem eines Wasserhydrometeors um ca. 30 % größer sein, um diesen Unterschied zu 
kompensieren. 



2.2 Dämpfung 

Die Dämpfung elektromagnetischer Strahlung in der Atmosphäre wird wesentlich durch Gase 
und Hydrometeore bestimmt. Signifikant dämpfende Gase sind dabei lediglich Sauerstoff und 
Wasserdampf. Die Dämpfung dieser beiden Gase nimmt grob mit der Frequenz der Strahlung 
zu (Abb. 2.2). Für Radar-Strahlung geringer Frequenz (Niederschlagsradar) ist die Dämpfung 
pro Kilometer klein (< O.Ol dB/km). Oberhalb etwa 10 GHz befinden sich erste ausgeprägte 
Absorptionsbanden des Wasserdampfs und des Sauerstoffs. Die Wahl der typischen Radar- 
frequenzen für Wolkenradare wird in Abhängigkeit dieser Absorptionsbanden getroffen und 
orientiert sich an den so genannten Absorptionsfenstem zwischen den Absorptionsbanden, 
bei 35 GHz und 95 GHz . Mit der Annahme einer homogenen Wasserdampf- Verteilung von 
7.5 g/m^ (entspricht Wasserdampfsättigung bei etwa 6.5 °C) ergibt sich nach Sauvageot (1992) 
für Messungen mit einem 35 GHz-Radar (MIRA36-S) eine maximale Dämpfung durch atmo- 
sphärische Gase von weniger als 0.1 dB/km. Bei MIRA36-S wird der Radarstrahl daher, bei 
maximaler Reichweite von 15 km, auf dem Weg zu den Streuern in 15 km Entfernung und wie- 
der zurück um bis zu « 3 dB gedämpft. Hydrometeore absorbieren ebenfalls elektromagneti- 
sche Strahlung und tragen zur Dämpfung bei. Betrachtet man zunächst die Absorption kleiner 
Hydrometeore, bzw. Wolkenpartikel, ergibt sich der Absorptionsquerschnitt aa in Rayleigh- 
Näherung für kleine Tropfen zu (Doviac und Zmic, 1984): 

a,=^^lmi-K)J2Dl (2.4) 

n=l 
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Abbildung 2.2: Dämpfung in dB/km (Einweg) durcli Sauerstoff (durchgezogene Linie) und 
Wasserdampf (7.5 g/m^, gestrichelte Linie) für Normaldrucl< (1013 hPa) als Funktion der Fre- 
quenz in GHz, aus Sauvageot (1992). 



Dabei bezeichnet Im den Imaginärteil und N ist die Anzahl der Tropfen, die zur Dämpfung bei- 
tragen. Aus Gl. 2.4 geht hervor, dass Verluste durch Absorption in Rayleigh-Näherung propor- 
tional zum Flüssigwassergehalt (M = Tr/6p [g/m^], P ist die Dichte in g/m^) entlang 
des Strahlwegs sind. Im Vergleich zum Streuquerschnitt ist der Parameter D/A um 3 Potenzen 
geringer vertreten. Dadurch ist die absorbierte Energie besonders bei kleinen Tropfen deutlich 
größer als die gestreute. Bei großen Tropfen (D/X > 1/10) kehrt sich diese Eigenschaft um, 
und die gestreute Energie übersteigt die absorbierte. Während die Dämpfung von Wolkeneis 
gering ist (< O.Ol (dB/km)/(g/m^), Sauvageot (1992)), kann die Dämpfung von Wolkentropfen 
0.1 (dB/km)/(g/m^) überschreiten, abhängig von der Radar- Wellenlänge. Dabei gilt: Je kürzer 
die Wellenlänge, desto höher die Dämpfung. Für MIRA36-S liegt die Dämpfung an Wolken- 
tropfen abhängig von der Temperatur unter der Annahme von Rayleigh-Näherung für einen 
Luftdruck von 1013 hPa zwischen 0.5 (dB/km)/(g/m^) und 2 (dB/km)/(g/m^) (Meneghini und 
Kozu, 1990; Sauvageot, 1992). Zu beachten ist allerdings, dass die Dämpfung in der Rayleigh- 
Näherung nur für Durchmesser von D < A/ 100 für Wassertropfen und D < X/2 für Eis- 
kugeln gültig ist und für größere Durchmesser unterschätzt wird (Blahak, 2005). Die stärkste 
Dämpfung elektromagnetischer Strahlung ist im Niederschlag zu beobachten. Die großen Hy- 
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drometeore können dabei, abhängig von der Radar- Wellenlänge, eine Dämpfiang von mehreren 
(dB/km)/(g/m^) verursachen. 



2.3 Radargleichung 

Der Zusammenhang zwischen der gesendeten Leistung Pt und der empfangenen Leistung Pr 
für Einzelstreuer wird mit der Radargleichung für meteorologische Radare hergestellt. Verschie- 
den ausführliche Formen der Gleichung sind in der einschlägigen Literatur beschrieben (z.B. 
Battan, 1973; Sauvageot, 1992; Rinehart, 1981; Doviac und Zmic, 1984). Die Radargleichung 
ist in Kurzform gegeben zu (Probert- Jones, 1962): 

Dabei ist C eine Konstante, die von den Radareigenschaften abhängt, wie beispielsweise die 
Wellenlänge A, der Antennengewinn G, der Antennenöffnungswinkel da, Verluste am Sender 
und Empfänger. Der Antennengewinn ist ein Maß für die Richtwirkung und den Wirkungs- 
grad einer Antenne. Hierfür wird die abgestrahlte Leistung der Antenne mit einer verlustlosen 
Bezugsantenne (homogene Abstrahlung in alle Richtungen) gleicher Sendeleistung verglichen, 
er ist der Rückstreuquerschnitt eines einzelnen Streuers in der Rayleigh-Näherung (Gl. 2.1). 
Der Index i bezieht sich dabei auf einen einzelnen Streuer innerhalb eines Volumenausschnitts. 
Dieser Volumenausschnitt wird auch Range Gate genannt. Der Begriff Range Gate findet in 
dieser Arbeit mehrfach Verwendung zur Beschreibung der Volumenausschnitte innerhalb ei- 
nes Radarstrahls. Der Volumenausschnitt ergibt sich aus der radialen Auflösung von Afi und 
der Fläche, die abhängig von der Entfernung R zum Radar mit dem rotationssymmetrischen 
Antennenöffnungswinkel 6a aufgespannt wird: 

7,- 5)2 

Ay = Afi— (2.6) 
8 In 2 

Die Summe über alle im Range Gate befindlichen Streu-Elemente ist die Reflektivität r]: 



PCO 

T]= a{D)Np{D)dD 
Jo 

= --\K^ \ Np{D)DUD (2.7) 



Np{D) ist die Partikelgrößenverteilung und das Integral darüber ist der Reflektivitätsfaktor 
Z. Dieser wird in den Einheit mm^/m^ angegeben. Da sich die Werte von Z über mehrere 
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Größenordnungen erstrecken können, wird in der Regel der Logarithmus des Verhältnisses von 
Z zu Zo = 1 mm^/m^ verwendet: 

C = 10 logio ^ (2.8) 

Zu beachten ist, dass der Herleitung der Radargleichung und der Anwendung folgende verein- 
fachende Annahmen zugrunde liegen: 



• Der Wert für \K'^ \ ist fest. Dadurch ist der Aggregatzustand und die Temperatur der Streu- 
er im Mess Volumen festgelegt. 

• Die Streuer sind klein im Vergleich zur Wellenlänge (Rayleigh-Näherung). 

• Mehrfachstreuung wird vernachlässigt. 

• Konstante Hardwarekomponenten, die Wellenlänge und die Sendeleistung, werden vor- 
ausgesetzt. 

• Der Reflektivitätsfaktor gilt im gesamten Streuvolumen gleichermaßen. Es kann keine 
Feinstruktur des Streuvolumens gemessen werden. Gemessen wird der Mittelwert der 
Reflektivität in einem Streuvolumen. 

Dies sind ebenfalls Fehlerquellen, die bei der Interpretation des Gemessenen berücksichtigt 
werden sollten. 



2.4 Geschwindigkeit von Streuelementen 

Neben der empfangenen Leistung und der daraus folgenden Reflektivität, kann mit den meisten 
Radaren heutzutage auch die Geschwindigkeit von Streuelementen bestimmt werden. Die mit 
Doppler-Radaren gemessene Radialkomponente der Geschwindigkeit der Streuer Vr,s setzt sich 
dabei zusammen aus der dreidimensionalen Windgeschwindigkeit vw und der Sedimentations- 
geschwindigkeit wp der Streuer, wobei letztere von der Form und Dichte der Streuer abhängt. 



(2.9) 

vs = ß vw — Wp k 



Hierbei ist eV der Einheitsvektor in radialer Richtung vom Radar zum Streu-Ort. Größere Par- 
tikel folgen nicht trägheitslos dem Wind, dies ist mit dem Korrekturfaktor ß berücksichtigt. 
Kann die Antenne nicht bewegt werden (trifft weltweit auf die meisten Wolkenradare zu), sind 
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vertikale Doppler-Messungen möglich und die Vertikalgeschwindigkeit der Streuer wird ge- 
messen. In diesem Fall spielt die horizontale Windkomponente keine Rolle und Gleichung 2.9 
vereinfacht sich zu Wr.s = vs = ß ww — wp (Abb.2.3). Ist die Antenne beweglich (MIRA36-S) 
werden radiale Geschwindigkeiten abhängig vom Azimut (et) und von der Elevation (e) gemes- 
sen (Abb.2.3). Dabei wird die Horizontalkomponente der Streuergeschwindigkeit hauptsächlich 
durch den vorherrschenden Horizontalwind bestimmt. 




Abbildung 2.3: Schematische Darstellung der Doppler-Messung der Radialgeschwindigkeit 
der Streuer und deren Abhängigkelten am Beispiel von ^ = i, i% = (7 m/s, 0 m/s, 1.5 m/s) und 
wp = 1 m/s für Doppler-Messungen im Zenit (grau) und in e = 45° (schwarz) 

Die Messung der Geschwindigkeit der Streuer v^^s erfolgt über die Differenz der Phasenän- 
derung zum Zeitpunkt ti des empfangenen Pulses im Vergleich zum gesendeten Puls und der 
Phasenänderung des nachfolgenden Pulses zum Zeitpunkt ii+i. 

2.4.1 Entfaltung 

Für Radargeräte, wie MIRA36-S, die eine gepulste elektromagnetische Welle abstrahlen, ist, 
abhängig von der Pulswiederholrate /p = 1/ (t^+i — ti) und der Wellenlänge A, die eindeutige 
maximal detektierbare Geschwindigkeit (s. Gl. 2.10) begrenzt und eine Faltung der gemesse- 
nen Radialgeschwindigkeit tritt auf. Gefaltete Geschwindigkeiten erschweren die Interpretation 
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der gemessenen Geschwindigkeiten und deren weitere Analyse. Eine Auswertung der Doppier- 
Messungen wird erst möglich, wenn diese zuvor entfaltet werden. 

Faltung tritt auf, wenn sich die Phase zwischen zwei Pulsen um mehr als ±180 ° ändert und die 
maximal eindeutig messbare Geschwindigkeit überschritten wird. Diese Geschwindigkeit, die 
so genannte Nyquist-Geschwindigkeit 

X fp c fp 

-N = ^ = ^, (2.10) 

wird durch die maximale eindeutige Phasenverschiebung von ±tt begrenzt. Hierbei ist c die 
Lichtgeschwindigkeit, fp die Pulswiederholrate, A die Wellenlänge und / die Frequenz des 
gesendeten Signals. Für MIRA36-S (/ = 35.5 GHz und /p = 5 kHz) ergibt sich eine Nyquist- 
Geschwindigkeit von vn = 10.56m/s (siehe auch Kap. 3.4.3). Abhängig vom Scanverfahren 
können verschiedene Entfaltungsmethoden angewendet werden. In dieser Arbeit werden zwei 
Methoden verwendet. Einerseits wird die Entfaltung von Geschwindigkeiten bei konstanter Ele- 
vation e und variablem Azimut a (so genannte PPI-Scans, siehe auch Kap. 3.2) und andererseits 
die Entfaltung von Dopplergeschwindigkeiten bei variabler Elevation e und konstantem Azimut 
a (so genannte RHI-Scans, siehe auch Kap. 3.2) angewandt. Die Entfaltung der Geschwin- 
digkeitsdaten von PPI-Scans wird mit Hilfe eines erweiterten VAD-Verfahrens (Kap. 2.4.2) 
ermöglicht und ist im anschließenden Kapitel dokumentiert. 

Die Entfaltung der Geschwindigkeitsdaten von RHI-Scans basiert im Wesentlichen auf einem 
Verfahren nach Hennington (1981). Dabei werden gefaltete Geschwindigkeiten mit Hilfe eines 
radialen Hintergrundwindfeldes Vb{R, e) (b für background) identifiziert und entfaltet. Die Be- 
rechnung des radialen Hintergrundwindfeldes basiert auf einem Profil des Horizontalwinds aus 
Modelldaten oder Messungen (z.B. von Radiosondenaufstiegen, Radar-, Lidarprofilern). 

Zur Detektion der Orte der Faltung wird eine Faltungszahl n{R, e), nach Hennington (1981) 
bestimmt. Sie ist definiert zu 

f,^j^^,^^MRA^^!iARA (2.11) 

Mit der Annahme, dass die Differenz zwischen den Werten des Hintergrundwindfeldes und dem 
wahren Windfeld klein sind, verglichen mit 2 (21.12 m/s für MIRA36-S), sind die Faltungs- 
zahlen dicht um ganze Zahlen gruppiert. Daher können die Faltungszahlen n durch Rundung 

ganzen Faltungszahlen n (so genannte Nyquistzahlen) zugeordnet werden. Die Richtung der 
Faltung und die Anzahl der Faltungen wird somit durch ganze Zahlen n wiedergegeben. Die 
wahre, entfaltete Geschwindigkeit kann dann berechnet werden zu 

Vr{R, e) = Vr,s{R, e) + 2 n{R, e) vn (2.12) 
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2.4.2 VAD-Verfahren 



Messungen der radialen Streuergeschwindigkeit Vr in einem festen radialen Abstand R vom Ra- 
dar unter einem Elevationswinkel e als Funktion des Azimuts a (PPI-Scan, siehe auch Kap. 3.2) 
bilden die Basis des VAD- Verfahrens. Mit der radialen Streuergeschwindigkeit ist es möglich 
die Windrichtung, sowie Winddrehungen aus einem so genannten PPI(V)-Diagramm (V steht 
für Velocity) per Auge abzuleiten. Dies lässt sich auch automatisieren und mit dem VAD- 
Verfahren detaillierter darstellen. Das VAD-Verfahren wurde von Browning und Wexler (1968) 
für meteorologische Radare großer Wellenlänge (Niederschlagsradare) entwickelt und ist auch 
auf Messungen mit Wolkenradaren übertragbar. Ziel ist es, mit diesem Verfahren einzelne Ge- 
schwindigkeitskomponenten abzuleiten. Dabei wird angenommen, dass sich die horizontale 
Geschwindigkeit linear entwickeln lässt. Dann ergibt sich für die, vom Azimut abhängige, Ge- 
schwindigkeit vs{ol): 



vsipt) 



^ + sin a + 1^ i? cos q; ^ 



?;o + 1^ Ä sin a + 1^ ß cos a 



V 



W — Wr, 



(2.13) 



Dabei ist das Koordinatensystem so angelegt, dass die Abszisse nach Osten und die Ordina- 
te nach Norden zeigt (a = 0 ° ist Norden). Die Sedimentationsgeschwindigkeit Wp, sowie 
die Änderungen der horizontalen Geschwindigkeitskomponenten in x- und y-Richtung du/dx, 
du/dy, dv/dx, dv/dy werden als konstant angenommen. 



Die radiale Geschwindigkeit, die im VAD dargestellt ist, ergibt sich zu: 



Vr,s{(y) = vs{a) er{a) 



(2.14) 



mit 



/ . \ 

sm a cos e 



cos CK cose 



v 



sm e 



(2.15) 



was der Strahlrichtung des Radars entspricht. Das Resultat lässt sich als Fourierreihenentwick- 
lung von Vr,si(x) schreiben (Abbruch nach zweitem Reihenglied): 



Vr,s{oi) = y + cos(nQ;) + 6„ sin(nQ;)) 



(2.16) 



n=l 
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mit 



R cos(e) [— + — ] ^-2{w -Wp) 




sin e 



Vq cos e 



uo cos e 



(2.17) 




Nur wenn auf dem gesamten Kreisring Daten vorhanden sind, ist es sinnvoll die horizontale 
Divergenz (V^ • = du/dx + dv/dy), die in Konkurrenz zur vertikalen Geschwindigkeits- 
komponente {w — Wp) steht, die Streckungsdeformation {dv/dy — du/dx) und die Scherungs- 
deformation (du/dy + dv/dx) zu berechnen, da die Fehler sonst zu groß werden. Zuverlässiger 
gestaltet sich die Berechnung des Horizontalwindes (uq und Vq). Sind wenig Messdaten auf ei- 
nem Kreisring vorhanden, wird statt der Fourieranalyse eine Minimierung zwischen den gemes- 
senen Werten und der Fitkurve durchgeführt, indem man uq und vq variiert. Der Horizontalwind 
ist dann repräsentativ für jene Bereiche auf einem Kreisring, in denen Streuer delektiert werden 
konnten. 

Entfaltung Tritt für manche Range Gates auf einem Kreisring Faltung auf, verfälschen die- 
se Werte die Berechnung der Horizontalgeschwindigkeiten. Um dem zu entgehen, wird das 
Verfahren durch Tabary et al. (2001) erweitert und die Fourierentwicklung erster Ordnung der 
Geschwindigkeit w^.s nach dem Azimut differenziert: 



Faltungssprünge lassen sich dadurch leicht detektieren, da sie große, unplausible differentiel- 
le Radialgeschwindigkeiten erzeugen. Zur Entfaltung der restlichen differentiellen Radialge- 
schwindigkeiten berechnen Tabary et al. (2001) zuerst den besten Fit an die gemessene Kurve 
und erhalten damit Werte für «i und bi. Diese werden in die modifizierte Gleichung 2.16 (Ab- 
bruch nach dem ersten Reihenglied und üq = 0) eingesetzt und die daraus berechnete Kurve 
der Radialgeschwindigkeiten als Funktion des Azimuts mit den gemessenen Radialgeschwin- 
digkeiten verglichen. Werte der gemessenen Radialgeschwindigkeiten, die mehr als | i^at | (Gl. 
2.10) von der berechneten Kurve abweichen, werden um 2 | ujv | in die richtige Richtung ver- 
schoben. Damit liegen die gemessenen Radialgeschwindigkeiten entfaltet vor und können mit 
dem VAD-Algorithmus nach Browning und Wexler (1968) ausgewertet werden. 



dVr,S 

da 



Ol cos a + bi sin a 



(2.18) 
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3 Wolkenradar MIRA36-S 

3.1 Technische Eigenschaften 

Die in diesem Kapitel beschriebenen Eigenschaften basieren im Wesentlichen auf dem Hand- 
buch des Herstellers des Wolkenradars MIRA36-S. MIRA36-S (35.5 GHz) ist ein scannen- 
des, Dopplerfähiges, polarimetrisches Wolkenradar, das zwei symmetrische Empfänger für den 
simultanen Empfang des Co- und Cross-polarisierten Signals hat. Die Cassegrain-Antenne 
ermöglicht Antennenbewegungen mit Azimutwinkeln a von 0 ° bis 360 ° und Zenitwinkeln 
ez bis zu 45 °. In Tab. 3.1 sind die wichtigsten technischen Details zusammengestellt. Manche 
Radarparameter, wie Pulsdauer r und Pulswiederholrate fp sind innerhalb bestimmten Grenzen 
variierbar. Die in Tab. 3.1 fett gekennzeichneten Werte wurden für die in dieser Arbeit analy- 
sierten Daten gewählt. Weitere Radargrößen werden im Folgenden angesprochen. 

Die Wellenlänge A des gesendeten Signals ist 

A = y = 8.44 mm (3.1) 

mit der Vakuum-Lichtgeschwindigkeit c. Die Pulslänge ergibt sich zu 

= CT = 59.96 m « 60 m (3.2) 

Die Entfernungsauflösung wird aus der zeitlichen Differenz zwischen gesendetem und empfan- 
genem Puls berechnet: 

Aß= y = 29.98 m« 30 m (3.3) 

Für eine eindeutige Bestimmung des maximalen Messbereichs muss der Puls reflektiert, bzw. 
am Empfänger gemessen worden sein, bevor der nächste gesendet wird. Es gilt: 

ßmax = ^ = = 30 km (3.4) 
Durch die Radarsoftware, die 512 Range Gates zulässt, wird die Reichweite auf 15 km begrenzt. 
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3 Wolkenradar MIRA3 6-S 



Radartyp 


monostatisch, gepulst, Magnetron 


Sendefrequenz / 


35.5 GHz 


Mittlere Sendeleistung Pt 


30 kW 


Antennenbauweise 


Cassegrain 


Antennendurchmesser 


1.2 m 


Antennengewinn G 


50.4 dB 


Antennenöffnungswinkel 9a 


0.52° 


Antennenpositionierungsgenauigkeit 


0.1° 


Scanbereich azimutal a 


0 ° bis 360 ° 


Scanbereich zenital e z 


-45 ° bis 45 ° 


Scangeschwindigkeit w«, uj^^, beide Richtungen 


67s,0°/sbis 10°/s 


Maximale Beschleunigung, beide Richtungen 


10° /s^ 


Leistungsverluste Fr 


7 dB 


Pulsdauer r 


100 ns, 200 ns, 400 ns 


Pulswiederholfrequenz (PRF) fp 


2.5 kHz, 5 kHz, 10 kHz 


Empfänger-Bandbreite Br 


5 MHz 


FFT-Länge tifft 


128, 256, 512, 1024 Puls-Zyklen 



Tabelle 3.1: MIRA36-S: Technische Eigenschaften, fett hervorgehoben sind die in dieser Arbeit 
verwendeten Einstellungen. 
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3.2 Abtastverfahren 

Eine Besonderheit des Wolkenradars MIRA36-S ist seine Eigenschaft, in gewissen Grenzen 
Azimut und Elevation der Antenne zu variieren (Scan-Fähigkeit). Mit diesen Eigenschaften 
kann das Wolkenradar einerseits eingesetzt werden, um konventionell vertikal zu messen, ande- 
rerseits aber auch zweidimensionale oder dreidimensionale Atmosphärenbereiche abzutasten. 
Für mehrdimensionale Messungen unterscheiden wir in RHI'- und PPI-Scans. Die Scange- 
schwindigkeit wurde hierfür auf 6 °/s festgelegt. 

Während eines RHI-Scans wird der Azimut konstant gehalten (meist in Richtung der mittle- 
ren Windrichtung) und die Antennenposition nur im Zenitwinkel verändert. Dieser Vorgang 
benötigt 15 Sekunden bei einer Scangeschwindigkeit von 6°/s. Dadurch wird ein vertikaler 
zweidimensionaler Kegelschnitt durch die Atmosphäre ermöglicht, und die Entwicklung einer 
Wolke im Vertikalschnitt kann über mehrere Scans verfolgt werden. Die Auflösung der in dieser 
Arbeit verwendeten RHI-Scans in der Elevation beträgt 



e = = 1.23 (3.5) 



Dabei ist uj^^ = 6°/s die Antennengeschwindigkeit in der Elevation, npFT = 256 die Anzahl 
der Pulszyklen pro FFT, = 4 die Anzahl der Spektren (FFTs) über die gemittelt wird und 
/p = 5 kHz die Pulswiederholfrequenz (PRF). 

Während eines PPI-Scans wird der Zenitwinkel konstant gehalten und der Azimut verändert 
(Kegelmantel-Schnitt). Ein voller PPI-Scan dauert 60 Sekunden bei einer Scangeschwindigkeit 
von Ua = 6°/s. Dieses Scanverfahren wird meist angewendet, um die horizontale Windrich- 
tung und Windgeschwindigkeit abzuschätzen (siehe Kap. 2.4.2). Die Auflösung im Azimut ent- 
spricht der Auflösung in der Elevation während eines RHI-Scans und beträgt Aa = 1.23 ° bei 

= 6°/s. 



3.3 Mittelung der Pulse 

Die Amplitude und Phase der gestreuten elektromagnetischen Wellen von allen Partikeln in- 
nerhalb des Radarpulses unterliegen Interferenzen und sind zufäUig. Atmosphärische Partikel 
sind unterschiedlich groß, nicht stationär und nicht im konstanten Abstand zu ihren Nachbar- 
partikeln. So lässt sich mit einem einzelnen Puls nur ein zufälliger Aufenthaltsort und eine 
turbulenzabhängige zufällige Geschwindigkeit der Hydrometeore in einem Range Gate erfas- 
sen. Ziel ist es allerdings, den mittleren Zustand der Streuer im Range Gate zu bestimmen. 
Deshalb ist es notwendig, eine gewisse Anzahl an Pulsen zu mittein, um den mittleren Zustand 



'Range Height Indicator 
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im Hinblick auf die Rückstreueigenschaften und Partikelgeschwindigkeiten im Range Gate zu 
erhalten. Eine wichtige Größe hierbei ist die Dekorrelationszeit. Idealerweise sollten mehrere 
voneinander unabhängige Werte gemittelt werden. Wenn man sich beispielsweise vorstellt, dass 
zwei fiktive Streuer in einem Range Gate vorhanden wären, würde eine doppelte Amplitude ge- 
messen werden, wenn deren Abstand genau A.R = A/2 wäre (bzw. eine Amplitude von 0, wenn 
deren Abstand genau AR = A/4 wäre). Wenn nur wenige aufeinander folgende Pulse gemittelt 
würden, würde dieser Extremwert das Ergebnis bestimmen. Die Zeit, die die Partikel innerhalb 
eines Range Gates benötigen, um sich umzugruppieren, ist die Dekorrelationszeit 

TD = 0.2— (3.6) 

Dabei entspricht cry der Standardabweichung der Teilchengeschwindigkeiten. Diese liegt ab- 
hängig von der Turbulenz nach Sauvageot (1992) zwischen ay = 0.5 m/s (geringe Turbulenz) 
und (Ty — 5 m/s (starke Turbulenz). Das heißt, dass die Dekorrelationszeit für A = 8.4 mm 
Werte zwischen 0.3ms (cry = 5 m/s) und 3.4ms (ay = 0.5 m/s) annehmen kann. Weiterhin 
ist aus Sauvageot (1992) bekannt, dass mehr als 30 unabhängige Pulse, bzw. 30 Dekorrelati- 
onszeiten notwendig sind, um den mittleren Zustand der Partikel im Range Gate messen zu 
können. 



Die Mittelung der Pulse erfolgt für MIR A3 6-S mit Hilfe der diskreten FFT über eine gewis- 
se Anzahl an Pulsen (71fft = 256, siehe Tab. 3.1) und einer anschließenden Mittelung (hier 
über n5 = 4 FFTs). Dies entspricht für die in dieser Arbeit durchgeführten Messungen einer 
Mittelung über nppT x ™s = 1024 Pulse und einer Mittelungszeit von 

tM = = 205 ms (3.7) 

jp 

Mit tM = 205 ms werden 60 bis 680 Dekorrelationszeiten (abhängig vom Grad der Turbulenz), 
bzw. unabhängige Pulse zur Mittelung verwendet. Diese Mittelung reicht aus, um die Amplitude 
des mittleren Eingangssignals repräsentativ für alle Streuer innerhalb eines Range Gates zu 
bestimmen. 



Aus dem Mittel über = 4 FFTs werden im nächsten Schritt drei spektrale Momente, die 
mittlere Leistung Pr, die Doppler- Geschwindigkeit Vj. und die spektrale Breite cr„^ berechnet. 
Dabei ist die mitüere Leistung proportional zum mittieren Quadrat der Amplitude des gemit- 
telten Eingangssignals. Die Doppler-Geschwindigkeit entspricht der mittleren radialen Streuer- 
geschwindigkeit innerhalb eines Range-Gates und die Doppler-Streubreite gibt die Inhomoge- 
nität der einzelnen Streuergeschwindigkeiten innerhalb eines Range Gates an. Die Doppler-Ge- 
schwindigkeiten und die Doppler-Streubreite werden in Strahlrichtung gemessen, abhängig von 
Zenitwinkel und Azimut. 
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3.4 Messgrößen 
3.4.1 Reflektivität 

Die Reflektivität wird nicht direkt gemessen, sondern muss über das Signal-Rauscli-Verliältnis 
(SNR^ mit der Radargleichung (Kap.2.3) berechnet werden. Das Signal-Rausch- Verhältnis ist 

SNR = — (3.8) 
Pn 

mit der mittleren empfangenen Leistung Pr (Gl. 2.5) und der thermischen Rauschleistung des 
Empfängers 

Pn = kB To Fr Br = kß (To + Tn) Br (3.9) 

Dabei ist kß = 1.38 10~^^ J/K die Boltzmann-Konstante, Tq = 290 K die Referenztemperatur 
zur Bestimmung der Rauschzahl Fr = 3.4 dB. Die Rauschtemperatur ist Tjv = Tq{Fr - 1) 
und Br = 5 MHz die Bandbreite des Empfängers (s. Tab. 3.1). 



Kombination der Gin. 2.5, 3.8, 3.9 liefert die Berechnung des Reflektivitätsfaktors 

_ SNR Pn R"^ _ Pr R'^ 

^- c -~c~ 

mit der Radarkonstanten 



(3.10) 



^ = ri0 1 O T X = 2 X 10 - (3-11) 

512 In 2 L,. m 

Die verwendeten Größen sind in Tab. 3.1 aufgeführt. Es folgt die Umrechnung des Reflekti- 
vitätsfaktors Z in das logarithmische Maß ( nach Gl. 2.8. In den folgenden Kapiteln dieser Ar- 
beit wird ausschließlich der logarithmierte Reflektivitätsfaktor besprochen und die Reflektivität 
rj taucht nicht wieder auf. Deshalb wird nachfolgend für den logarithmierten Reflektivitätsfaktor 
die im meteorologischen Sprachgebrauch übliche Bezeichnung Reflektivität verwendet. 



3.4.2 Lineares Depolarisationsverhältnis 

Das Lineares Depolarisationsverhältnis (LDR) zeigt an, wieviel Leistung im Cross-Kanal (Po- 
larisationsebene orthogonal zum gesendeten Signal), relativ zur Leistung im Co-Kanal (gleiche 
Polarisationsebene wie das gesendete Signal) empfangen wird. Nicht rotationssymmetrische 
Streu-Elemente, deren Hauptachse nicht in Richtung des elektrischen Feldes des gesendeten Ra- 
darstrahls ausgerichtet ist, depolarisieren einen kleinen Teil der Energie des gesendeten Signals, 
d.h. die Polarisationseben ändert sich. Dadurch wird im Co-Kanal und im Cross-Kanal jeweils 
ein Anteil des rückgestreuten Signals empfangen. Das lineare Depolarisationsverhältnis L [dB] 
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ist definiert als das Verhältnis der zurückgestreuten Leistung im Cross-Kanal zur rückgestreuten 
Leistung im Co-Kanal und wird bestimmt aus dem logarithmierten Verhältnis der beiden Re- 
flektivitätsfaktoren auf dem Co-Kanal (CO) und Cross-Kanal (CX): 



Das LDR ist abhängig vom Orientierungswinkel der Streu-Elemente. Es ist klein für ku- 
gelförmige Streu-Elemente, die unabhängig vom Orientierungswinkel sind. Das LDR erreicht 
maximale Werte für achsensymmetrische Streu-Elemente, die in einem Orientierungswinkel 
von 45 ° vorliegen. Betrachtet man ein Streuensemble, ist das LDR ein Maß für die Abweichung 
der Form der Streuer von der Kugelform. Hierbei spielt nicht nur die Form, sondern auch das 
Verhalten der Streuer während des Fallens eine Rolle. Je komplexer die Streu-Elemente inner- 
halb eines Range Gates sind, und je mehr fallende Streuer taumeln, desto größer wird das LDR. 
Typischerweise sind die Werte des LDR für Range Gates, in denen nur Regentropfen vorlie- 
gen, niedriger gegenüber den Werten für Range Gates, in denen sich ausschließlich Eispartikel 
befinden. Die höchsten Werte des LDR werden für nicht-meteorologische Streu-Elemente in 
der Grenzschicht (Insekten, Staub, usw.) und für nasse Eispartikel (Schmelzschicht) erreicht 
(Abb.3.1). 

3.4.3 Geschwindigkeitsgrößen 

Wie zuvor bereits beschrieben, strahlt das Wolkenradar Pulse elektromagnetischer Energie ab 
und empfängt ein Signal, das hauptsächlich durch meteorologische Ziele (z.B. Eiskristalle, - 
Partikel, Wassertröpfchen) rückgestreut wü-d. Das rückgestreute Signal kann als komplexes Si- 
gnal betrachtet werden (IQ-Signal). Aus der Veränderung der Phase des Signals lässt sich auf 
die (radiale) Geschwindigkeit schließen. 

Die radiale Dopplergeschwindigkeit läßt sich technisch nicht direkt aus der Dopplerverschie- 
bung eines Pulses bestimmen, denn selbst für sehr hohe Radialgeschwindigkeiten von z.B. 
30 m/s wäre die Dopplerverschiebung klein (fo ~ 7 kHz). Mit einer Sendepulsdauer von 
T = 200 ns ergibt sich eine Frequenzunbestimmtheit von 1/r = 5 MHz. Änderungen der Emp- 
fangsfrequenz von weniger als 5 MHz können demnach nicht direkt erfasst werden. 

Die radiale Dopplergeschwindigkeit läßt sich aber abhängig vom Radartyp über folgende 
Methoden abschätzen. Eine häufig verwendete Methode ist das so genannte „Pulse-Pair- 
Verfahren". Hierbei werden die Phasen zweier aufeinanderfolgender Pulse bestimmt. Die Pha- 
senänderung des komplexen Empfangssignals (IQ-Signal) von Puls zu Puls enthält dann die 
Information über die Bewegung der Ziele. Aus der Änderung der empfangenen Phase (fr mit 
der Zeit kann damit die Dopplergeschwindigkeit Vr berechnet werden: 




(3.12) 



dt 



= 2 Tt/d = —Vr,S ~ 



"KT 



(3.13) 
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Dabei ist fp (hier: 5 kHz) die Pulswiederholfrequenz, A = 8.44 mm die Radarwellenlänge und 
fo = —'^ t^r,s/^ die Dopplerfrequenz. 

Die Phase lässt sich eindeutig aber nur auf ganze Vielfache von 2 tt bestimmen, bzw. für die 
Phasendifferenz gilt der Eindeutigkeitsbereich — tt < Aip < tt. Nach Gl. 3.13 ergibt sich die 
maximale eindeutige Geschwindigkeit dann zu 

^^r,5,max = ^ = 10.56- (3.14) 

4 s 

Diese Geschwindigkeit ist die bereits in Kap. 2.4.1 (Gl. 2.10) eingeführte Nyquistgeschwindig- 
keit ■yjv- Die Beschränkung des Eindeutigkeitsintervalls führt bei hohen radialen Geschwindig- 
keiten zu einer Faltung. Auf dieses Problem wurde bereits in Kap. 2.4.1 eingegangen. 

Eine weitere, häufig verwendete Methode ist die Abschätzung der radialen Dopplergeschwin- 
digkeit aus der Dopplerfrequenz, die mit Hilfe einer schnellen Fouriertransformation (FFT) 
über aUe Pulse innerhalb eines eingestellten Zeitintervalls ermittelt wird. Dieses Verfahren wird 
im Wolkenradar MIRA36-S angewendet. Zur Bestimmung der Dopplerfrequenz wird das kom- 
plexe Signal Sk{l) im kten Range-Gate zur Zeit to + l At als Zeitreihe interpretiert. Dabei ist 
At = 1 / fp und / der Zählindex der Pulse. Dieser Zählindex wird begrenzt durch die eingestellte 
maximale Anzahl an Pulsen, die zur FFT beitiagen (hier: npFT = 256). Aus dem Leistungs- 
spektrum 



Pk{Sk,fm) OC 



"FFT 



(3.15) 



1=1 

kann dann die Dopplerfrequenz fo als Schwerpunkt des Hauptpeaks berechnet werden (Brin- 
gi und Chandrasekar, 2005). Dabei sind die Frequenzen des Leistungspektrums fm — m Af 
("^max = N)- Bei Kenntniss der Dopplerfrequenz ist es dann möglich die radiale Dopplerge- 
schwindigkeit nach Gl. 3.13 zu bestimmen. Auch hier ist durch die endliche zeitliche Auflösung 
die maximale, eindeutig zu bestimmende Geschwindigkeit Vr,s durch die Nyquistgeschwindig- 
keit (Gl. 2.10) beschränkt. 

Die Auflösung Av^.s, mit der die radialen Geschwindigkeiten aufgezeichnet werden können, ist 
abhängig von der Anzahl der Pulse, die zur FFT beitragen und ergibt sich zu: 

Avr,s = ^^^li^ = 0.08 - (3.16) 

J^FFT S 

Die spektrale Breite ist die empirische Varianz der Dopplergeschwindigkeit bzw. der FFT- 
Frequenz. Sie wird durch die Breite des Hauptpeaks im Leistungsspektrum bestimmt. Diese 
kann groß werden durch Turbulenz, wenn die Größenverteilung der Streuer breiter wird, ein 
Doppelpeak oder ein Faltungsübergang (siehe Kap. 2.4.1) existiert. Ein Doppelpeak kann ge- 
messen werden, wenn Streuer unterschiedlicher Eigenschaften innerhalb eines Range Gates 
vorhanden sind, z.B. Wassertropfen und Eispartikel. 
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3.4.4 Überblick 

Einen ersten Einblick in die Messgrößen liefern die Abbildungen von Messungen in vertikaler 
Richtung mit M1RA36-S vom 1. Juli 2007 (Abb. 3.1). Eine Schmelzschicht ist in der Höhe von 
ca. 2 km (zw. 12:00-16:00 UTC) anhand folgender Auffälligkeiten zu erkennen: erhöhte Reflek- 
tivitäten ( und erhöhte LDR-Werte L, sowie eine Grenzlinie zwischen langsam fallenden Par- 
tikeln und schneller fallenden Partikeln in der Vertikalgeschwindigkeit Vr,s und einer erhöhten 
spektralen Breite cr„^^ unterhalb 2 km aufgrund der erhöhten Sedimentationsgesch windigkei- 
ten der Streuer. Auf das Thema Schmelzschicht wird im nachfolgenden Kapitel 4 detaillierter 
eingegangen. Unterhalb der Schmelzschicht regnet es, das LDR nimmt die niedrigsten Werte 
an (runde Streuer) und die Sedimentationsgeschwindigkeit ist am größten. Streuer, die keine 
Hydrometeore sind, werden durch räumlich stark variable und zugleich hohe LDR-Werte un- 
terhalb einer Höhe von 2 km angezeigt. Diese Streuer sind meist Insekten und Staubpartikel. 
Weiterhin sind in der Abbildung Vertikalmessung der Radialgeschwindigkeit Vr,s Aufwindbe- 
reiche (gelb) und eine Mischung aus Abwindbereich und Sedimentationsgeschwindigkeit (grün, 
blau) in Wolken zu erkennen. Die Entwicklung eines einzelnen Aufwindbereichs kann alleine 
mit Vertikalmessungen nicht beobachtet werden. Dies wird durch zweidimensionale Messun- 
gen (RHI-Scans) ermöglicht und ist Gegenstand des Kapitels 5. 

3.5 Datengrundlage 

Seit 2007 wurden mit dem Wolkenradar MIRA36-S neben vertikalen Messungen auch RHI- 
Scans und PPI-Scans, sowie 3D-Messungen durchgeführt. Vom 25. Mai bis 10. September 2007 
wurde es während der „Convective-and-Orographically-induced-Precipitation-Study" (COPS, 
siehe Wulfmeyer et al., 2008; Kottmeier et al., 2008) eingesetzt, einer Messkampagne auf 
dem höchsten Berg des Nordschwarzwaldes, der Hornisgrinde (1164 m). Dabei wurden un- 
ter anderem langsame RHI-Scans (Geschwindigkeit der Antennenbewegung im Zenitwinkel 
dez/dt < 4°/s) in Windrichtung, orthogonal dazu und in Richtung weiterer Messsystem- 
Standorte im Westen (Rheintal) und Osten (Murgtal) durchgeführt. Der kürzeste zeitliche Ab- 
stand zwischen den jeweiligen RHI-Messungen betrug 4 Minuten. Vom 10. September 2007 
bis 31. Januar 2010 wurde das Wolkenradar auf der Messwiese am KIT Campus Nord einge- 
setzt und hauptsächlich schnelle RHI-Scans (dez/dt = 6°/s) in Windrichtung in zeitlichen 
Abständen von 30 Sekunden (RHI-Scans in die gleiche Richtung) gefahren. Damit konnte 
erstmals die zeitliche Entwicklung von Wolkenbereichen aufgezeichnet werden. Zwischen Juli 
2007 und Februar 2010 wurden insgesamt 37016 RHI-Scans durchgeführt. Darin enthalten sind 
auch Messungen ohne Wolken. Die Anzahl der Messungen, bei denen in 30% der Messfläche 
eines RHI-Scans bis zu einer Höhe von 10 km (meist befinden sich Wolken unterhalb dieser 
Höhe) Streuer delektiert werden konnten, liegt bei 8943. Detailliertere Angaben sind der Tabel- 
le 3.2 zu entnehmen. 
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Abbildung 3.1: Zeit-Höhen-Schnitte der Reflektivität C[dBz], des linearen Depolarisations- 
verhältnises L [dB] (dunkelrot eingefärbt sind Bereiche, in denen Werte im Co-Kanal, aber nicht 
im Cross-Kanal vorliegen, d.h. kein LDR berechnet werden konnte), der radialen Dopplerge- 
schwindigkeit Vr,s [m/s] und der spektrale Breite a^^g [m/s]. Details im Text. 
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Jahr 


Anzahl 


DJF 


MAM 


JJA 


SON 


2007 


RHI 


- 


- 


3274 


898 


RHI30 






400 


83 


2008 


RHI 




837 


159 


428 


RHI30 






105 


235 


2009 


RHI 


131 


147 




14343 


RHI30 


128 






4184 


2010 


RHI 


16799 








RHI30 


3808 








E 


RHI 


16930 


984 


3433 


15669 


RHI30 


3936 




505 


4502 



Tabelle 3.2: Anzahl gemessener RHI-Scans und Anzahl solcher RHI-Scans (RHI30), bei denen 
in mindestens 30% einer definierten Messfläche (200 km^), aufgespannt durch einen Zenitwin- 
kel (.z von 45° in beiden Richtungen und einer Höhe von 10 km, Streuer gemessen werden. 
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4 Schmelzschicht in Wolken 

Fallen Eispartikel durch warme Luftschichten (T > 0 ° C), beginnen sie zu schmelzen. Die 
Änderung der Form und des Aggregatzustandes der Hydrometeore führt zu Modifikationen ih- 
rer Streueigenschaften und ihrer Sedimentationsgeschwindigkeit. Dies äußert sich in einer in 
dieser Schicht erhöhten Reflektivität, für die mehrere Prozesse verantwortlich sind. Untersu- 
chungen der Schmelzschicht werden hauptsächlich aus zwei Gründen durchgeführt: Erstens sol- 
len die Prozesse, die zu den charakteristischen Radarsignaturen in der Schmelzschicht führen, 
ermittelt und ihr Zusammenspiel verstanden werden. Dies ist auch die Motivation, die diesem 
Kapitel zugrunde liegt. Zweitens sollen Algorithmen entwickelt werden, um diese für die Nie- 
derschlagsmessung am Boden störenden Artefakte zu eliminieren. 

Radarbeobachtungen der Schmelzschicht gibt es schon seit den Anfängen der Radarmeteoro- 
logie (Austin und Bemis, 1950). Detektionsalgorithmen wurden bereits für Radargeräte mit 
vertikal ausgerichteten Antennen als auch für Scan-fähige polarimetrische Radargeräte ent- 
wickelt. Detaillierte Untersuchungen der Schmelzschicht selbst stammen von Fabry und Za- 
wadzki (1995). Später wurden auch Detektionsalgorithmen entwickelt, die entweder Bestand- 
teile von Klassifikationsalgorithmen sind, z.B. in der jüngsten Vergangenheit von Giangran- 
de et al. (2008), Hogan und O'Connor (2004), White et al. (2002), oder die zur besseren 
Abschätzung der Regenrate dienen, z.B. Gysi et al. (1997), Matrosov et al. (2007). 

Mit einem X-Band-Radar (A = 3.2 cm) mit rein vertikal ausgerichteter Antenne und einem 
Windprofiler-Radar (A = 32.8 cm) arbeiteten beispielsweise Fabry und Zawadzki (1995), de- 
ren Ziel die Untersuchung der Schmelzschicht anhand von gemessenen Radarmessgrößen im 
Bereich der Schmelzschicht und in den Schichten direkt darüber und darunter war. Die Detek- 
tion der Schmelzschichtgrenzen und der Höhe über Grund, in der die Reflektivitäten innerhalb 
der Schmelzschicht maximal werden (Schmelzschichthöhe), erfolgte dabei über die maximalen 
Krümmungen der Kurve der minutengemittelten Vertikalprofile der Reflektivität. Dafür wurden 
die Profile allerdings mit den Daten des Windprofilers korrigiert, wodurch die Profile Trajekto- 
rien der Hydrometeore auf ihrem horizontalwindabhängigen Weg nach unten entsprachen. 

Vertikal zeigende Radargeräte (35 GHz, 94 GHz) wurden ebenfalls von Hogan und O'Connor 
(2004) eingesetzt. Sie erarbeiteten einen Algorithmus zur Schmelzschichtdetektion, der einge- 
bettet ist in einen umfangreichen Klassifikationsalgorithmus, der Radardaten, Lidardaten, Mi- 
krowellenradiometerdaten und ModeUdaten berücksichtigt. Die Schmelz schicht wird dabei mit 
Hilfe der vertikalen Dopplergeschwindigkeit und der Feuchttemperatur aus Modelldaten be- 
stimmt. Die Suche nach der Schmelzschicht wird eingeschränkt auf eine Schicht, in der die 
Feuchttemperatur Tf zwischen — 5°C < Tf < 5°C liegt. In dieser Schicht wird dann der 
stärkste Gradient der Vertikalgeschwindigkeit gesucht und als Schmelzschichthöhe bezeichnet. 
Ist die Vertikalgeschwindigkeit in der Schmelzschichthöhe kleiner 0.5 m/s, wird keine Schmelz- 
schicht angenommen. Ausgehend von der Schmelzschichthöhe wird in beiden Richtungen nach 
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Höhen gesucht, in denen die Vertikalgeschwindigkeit und der Gradient der Vertikalgeschwin- 
digkeit definierte Grenzwerte unterschreiten. Diese Höhen entsprechen dann der Ober- und Un- 
tergrenze der Schmelzschicht. Am Ende werden zeitlich aufeinander folgende Profile mitein- 
ander verglichen. Profile, deren Schmelzschichthöhe um mehr als 150 m von den benachbarten 
Schmelzschichthöhen abweicht, werden von der weiteren Betrachtung ausgeschlossen. 

Die Schmelzschicht wurde auch mit einer Kombination von Messgrößen aus PPI-Messungen 
polarimetrischer Radargeräte, beispielsweise von einem WSR-88D-Radargerät (S-Band, Ame- 
rikanischer Radarverbund) durchgeführt (z.B. Giangrande et al., 2008; Vivekanandan et al., 
1999; Park et al., 2008). Sie entwickelten über viele Jahre einen aufwendigen Klassifika- 
tionsalgorithmus. Dieser berücksichtigt Messfehler und baut auf Prinzipien der Fuzzylogik auf. 
Hierfür werden sich überschneidende Wertebereiche fürpolarimetrische Messgrößen festgelegt, 
die auf Regen, Schnee, Graupel oder nassen Schnee hinweisen. Die Schmelzschichtgrenzen 
werden dadurch festgelegt, dass für jede Elevation und Höhe ein Histogramm der polarimetri- 
schen Größen erstellt wird. Wenn der Wertebereich für schmelzenden Schnee der dominierende 
Wertebereich ist, werden das 20. und 95. Perzentil verwendet, um die Ober- und Untergrenze 
der Schmelzschicht zu bestimmen. 

Eigenschaften der Schmelzschicht sind auch in den gemessenen Größen des Wolkenradars, spe- 
ziell der Reflektivität, dem linearen Depolarisationsverhältnis und der Dopplergeschwindigkeit, 
zu erkennen. Dadurch kann die Höhe und Mächtigkeit der Schmelzschicht mit diesen Messun- 
gen abgeschätzt werden. Die Kenntnis der genauen Lage der Schmelzschicht ist wichtig für 
eine realistische Trennung von Bereichen, die durch Eispartikel, Tropfen oder der Mischphase 
bestimmt werden. Anhand der Schmelzschichtoberkante kann zudem die Höhe der Nullgrad- 
grenze abgeschätzt werden. 



4.1 Grundlagen 

Da in dieser Arbeit die räumliche und zeitliche Änderung dynamischer und wolkenphysi- 
kalischer Prozesse in Wolken im Vordergrund steht, wurden zweidimensionale Messungen 
durchgeführt (RHI-Scans). Mit diesen zweidimensionalen Messungen soll es möglich sein, 
die Schmelzschicht innerhalb eines Scans detektieren zu können. In diesem Kapitel wird ein 
Algorithmus vorgestellt, der ausschließlich auf den Messungen eines Wolkenradars, hier des 
MIRA36-S, beruht. Dafür werden räumlich gemittelte Messgrößen (Reflektivität, LDR und 
Vertikalgeschwindigkeit) kombiniert. Ziel ist es, in einer zeitlich hohen Auflösung (für jeden 
RHl-Scan, der 15 Sekunden dauert) eine Schmelzschicht zu detektieren. Der Algorithmus soll 
auch in Situationen, in denen nur zwei der drei genannten Größen zur Verfügung stehen, ei- 
ne Schmelzschicht detektieren. Auf solch einen Fall wird im weiteren Verlauf noch gesondert 



4.1 Grundlagen 



33 



eingegangen. Die Ober- und Untergrenzen der Schmelzschicht sind, wie auch in den zuvor be- 
schriebenen Arbeiten, Abschätzungen der Schmelzschichtgrenzen. Mittelungen der Daten sind 
notwendig, damit die Höhe der Schmelzschicht möglichst solide ermittelt werden kann. 

Um die Schmelzschicht verlässlich detektieren zu können, ist es zudem wichtig, die Eigenschaf- 
ten der Schmelzschicht und ihren Einfluss auf die Messgrößen zu verstehen. Deshalb wird im 
Folgenden beschrieben, was das Wolkenradar im idealen Fall während der unterschiedlichen 
Phasen des Schmelzvorgangs messen würde. Stellvertretend für einen idealen Fall wird der 
Schmelzvorgang anhand einer RHI-Messung vom 06. 10.2009 beschrieben und dargestellt (Abb. 
4.1), wobei auf vertikale Profile der Reflektivität CC-^), des linearen Depolarisationsverhältnisses 
L{z) und der Vertikalgeschwindigkeit der Streuer Wm{z) (siehe Kap. 5.1) eingegangen wird. 

Während des Schmelzens von Eispartikeln wirken sich drei Vorgänge auf die gemessene Re- 
flektivität aus, auf die im weiteren Verlauf eingegangen wird: Der Dielektrizitätsfaktor der Hy- 
drometeore nimmt zu, der Durchmesser der Hydrometeore ändert sich und die Anzahl der Hy- 
drometeore pro Volumeneinheit nimmt ab. 

Am Oberrand der Schmelzschicht, der mit der 0 °C-Isotherme in Zusammenhang steht, begin- 
nen fallende Eispartikel zu schmelzen. Dieser Prozess beginnt an der Oberfläche der Hydro- 
meteore. Bei größeren Eisaggregaten, wie Schneeflocken, entstehen Lufteinschlüsse, da sich 
eine Hülle aus Flüssigwasser um die Eispartikel bildet. Wegen der dann geringeren mittleren 
Materialdichte der anschmelzenden Eishydrometeore ist ihr Durchmesser größer als der eines 
Tropfens der gleichen Wassermenge. Sind die Oberflächen der Eispartikel flüssig, ist die Aggre- 
gation, d.h. die Büdung von größeren Eispartikeln durch Stoßprozesse, effektiver. Die Partikel 
können sich nicht nur verhaken, wie während der Bildung von Eispartikeln durch Dendrite, 
sondern können auch aneinanderhaften. Die Hydrometeore nehmen dadurch komplexere For- 
men an, wodurch das LDR deutlich ansteigt und der Durchmesser der Hydrometeore zunimmt. 
Gleichzeitig nehmen der Brechungsindex und somit der Dielektrizitätsfaktor höhere Werte an 
(für den vollständigen Phasenübergang der Hydrometeore um den Faktor |-R^lwater/l-^lice ~ 
Durch die Vergrößerung des Durchmessers und Erhöhung des Dielektrizitätsfaktors steigt der 
Rückstreuquerschnitt und das Radar misst eine, im Vergleich zu Eispartikeln, höhere Reflekti- 
vität (. 

Im Weiteren wird davon ausgegangen, dass die angeschmolzenen großen Eispartikel wei- 
ter sedimentieren und dementsprechend der Schmelzvorgang weiter fortschreitet. Sind noch 
Verästelungen der Eispartikel vorhanden, kollabieren diese. Der Durchmesser der so entstehen- 
den Hydrometeore nimmt ab, sie werden runder und die LDR-Werte sinken. Infolgedessen sinkt 
der Reibungswiderstand der schmelzenden und immer runder werdenden Partikel, wodurch sich 
ihre Sedimentationsgeschwindigkeit erhöht. Die Folge ist eine Zunahme der Distanz zwischen 
den Hydrometeoren während die Anzahl der Hydrometeore pro Volumeneinheit abnimmt. Da 
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der Reflektivitätsfaktor Z proportional zur Anzahl der Streuer ist, folgt dann eine Abnahme der 
Reflektivität C- 

Ist der Schmelzprozess abgeschlossen, sind die Eispartikel geschmolzen und im Unterschied 
zum Beginn des Schmelzprozesses näherungsweise kugelförmig (d.h. Tropfen). Gleichzeitig 
befinden sich keine Lufteinschlüsse mehr in den Hydrometeoren und die Materialdichte ist 
höher. Dadurch ist die Geschwindigkeit der Partikel am Ende des Schmelzprozesses höher und 
das LDR niedriger. Da der Effekt der Erhöhung des Dielektrizitätsfaktors gegenüber der Ver- 
ringerung der Anzahldichte und der Durchmessers der Hydrometeore meist überwiegt, sind die 
Reflektivitäten am Unterrand der Schmelzschicht meist höher als am Oberrand. 

Es ist allerdings schwierig, den Anteil des jeweiligen Vorgangs an der Gesamtänderung der 
Messwerte der Reflektivität, des LDRs oder der vertikalen Geschwindigkeit genau zu spezifi- 
zieren. Zusammenfassend kann festgehalten werden, dass die Form (z.B. kleine Eispartikel (Na- 
deln), zusammenhängende komplexe Eiskristalle (Dendrite), runde zusammenhängende Eispar- 
tikel (Graupe))), Anzahl und Größe der Eispartikel oberhalb der Schmelzschicht und der verti- 
kale Temperaturgradient einen großen Einfluss auf die Variation der Messgrößen haben dürften. 
Die erwähnten Einflussfaktoren wirken sich auf die Schmelzgeschwindigkeit der Eispartikel, 
die Änderung des Durchmessers und die Beschleunigung der Hydrometeore, aber nur mittelbar 
auf den Dielektrizitätsfaktor aus. Dadurch sind die Merkmale der Schmelzschicht nicht immer 
klar in den Profilen der Reflektivität, des LDRs oder der vertikalen Geschwindigkeit zu erken- 
nen. 

Ausgangspunkt für den Detektionsalgorithmus ist allerdings der oben beschriebene konzep- 
tionelle Verlauf. Der ideale Verlauf des Schmelzvorgangs, dem die Messung vom 06.10.2009 
sehr nahe kommt, ist in Abb. 4.1 dargestellt. Bevor diese Abbildung diskutiert wird, noch ein 
Wort zur Datenauswertung. Die gemessenen Werte für ^z) und L(i2, e^) wurden hierfür 
zunächst auf ein zweidimensionales kartesischen Gitter (x, z) interpoliert, mit Gitterabständen 
von Aa; = lS.z = 30 m in beiden Richtungen. Für jede Höhe z wurde dann ein Mittelwert aller 
finiten Werte von C,{x,z) und L{x,z) in x-Richtung berechnet, wenn mindestens zwei Drittel 
der möglichen Messwerte in einer Höhe besetzt waren. Im letzten Schritt wurde die Datenmen- 
ge zur besseren Detektion der Schmelzschichtgrenzen erhöht und die Profildaten auf ein Gitter 
mit einem Gitterabstand von Az = 10 m interpoliert. 

In der schon erwähnten Abb. 4.1 sind die Profile der Reflektivität C{z), des linearen Depo- 
larisationsverhältnisses L{z) und der Vertikalgeschwindigkeit der Streuer Wm{z) gezeigt. Die 
Beschreibung der relevanten Bereiche zur Schmelzschichtdetektion erfolgt für alle folgenden 
Abbildungen in Bewegungsrichtung der Partikel, also von oben nach unten. Die Werte für 
Wm{z) sind negativ, da sich die Partikel auf das Radar zu bewegen. Im stratiformen Fall (z.B am 
06.10.2009) ist davon auszugehen, dass die Vertikalgeschwindigkeit der Streuer hauptsächlich 
durch deren Sedimentationsgeschwindigkeit bestimmt wird. Deshalb wird eine Abnahme der 
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Abbildung 4.1: Vertikalprofile folgender Größen aus RHI-Messungen des Wolkenradars 
MIRA36-S vom 06.10.2009, 13:14 UTC: Reflektivität C(^) in dBz, lineares Depolarisations- 
verhältnis L{z) in dB und die Vertikalgeschwindigkeit Wm{z) in m/s. Grau unterlegt ist die de- 
lektierte Schmelzschicht. 

Vertikalgeschwindigkeit von oben nach unten einer zunehmenden Sedimentationsgeschwindig- 
keit zugeordnet. In Abb. 4.1 ist die delektierte Schmelzschicht grau unterlegt . 



Nun werden die maßgeblichen Charakteristika für die Detektion der Schmelzschicht im Umfeld 
des gesamten Messbereichs betrachtet. Der Messbereich in Abb. 4.1 wird von oben nach unten 
in drei Höhenabschnitte unterteilt, die sich an der grau unterlegten delektierten Schmelzschicht 
orientieren. Im Bereich oberhalb der Schmelzschicht liegen hauptsächlich Eiskristalle und - 
Partikel vor, die unterhalb 5 km wachsen. Ein Indiz für das Wachstum ist in der gleichzeitigen 
Zunahme der Reflektivität und der Sedimentationsgeschwindigkeiten zu sehen. Die Sedimen- 
tationsgeschwindigkeit der Hydrometeore ist in diesem Bereich gering, was darauf hindeutet, 
dass es sich um Eispartikel oder große Eiskristalle, etwa Dendrite, handelt. Einer deutlicheren 
Geschwindigkeitszunahme aufgrund der Zunahme an Masse steht, durch das meist flächenhafte 
Wachstum der Eiskristalle und Eispartikel, eine deutliche Erhöhung der Fläche und somit des 
Reibungswiders tands entgegen. 



Direkt darunter ist die Schmelzschicht zu erkennen, in der sich die oben beschriebenen 
Vorgänge abspielen. Merkmale, die deutlich hervorstechen, sind in allen Profilen zu erkennen 
und für die Detektion der Schmelzschicht nutzbar: Ein deutlich sichtbares lokales Maximum 
(Peak) im Profil der Reflektivität und im LDR und eine starke Änderung im Profil der Verti- 
kalgeschwindigkeit. In vielen Fällen ist die Spitze des Peaks gleichzeitig das absolute Maxi- 
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mum der Reflektivität und des LDR. In dieser Schicht liegen die Hydrometeore als zunehmend 
flüssiger werdende Eispartikel vor, deren Größenverteilung sich deutlich ändert. 

Im dritten Abschnitt, unterhalb der Schmelzschicht, sind die Eispartikel vollständig geschmol- 
zen, d.h. es liegen Tropfen vor, die näherungsweise eine sphärische Form besitzen und schnell 
fallen. Die Sedimentationsgeschwindigkeit liegt zwischen 4 m/s und 5 m/s und ist somit deut- 
lich höher als oberhalb der Schmelzschicht (1m/s bis 1.5 m/s). Die Werte der Reflektivitäten 
sind unterhalb der Schmelzschicht ebenfalls höher (4 dBz bis 7 dBz) als oberhalb der Schmelz- 
schicht (-7 dBz bis 3 dBz). Im Gegensatz dazu sind die LDR-Werte oberhalb der Schmelz- 
schicht höher (etwa -20 dB) als unterhalb (etwa -22 dB). 

Aufgrund der oben genannten Einflussfaktoren ist absehbar, dass nicht immer der ideale Fall 
vorliegt. Der Anspruch an den Detektionsalgorithmus ist allerdings, auch in solchen Fällen die 
Schmelzschicht zu detektieren. Ein solcher, vom „Ideal" abweichender Fall ist in Abbildung 
4.2 dargestellt, inklusive der detektierten Schmelzschicht. Auf die zum idealen Fall verschie- 
denen, meteorologischen Bedingungen soll an dieser SteUe nicht eingegangen werden, sondern 
lediglich die für den Algorithmus wichtigen Unterschiede der Profilmerkmale im Schmelz- 
schichtbereich aufgezeigt werden. Im Vergleich zum idealen FaU ist die Reflektivität oberhalb 
der Schmelzschicht höher als unterhalb, und der Peak innerhalb der Schmelzschicht ist zwar 
sichtbar, aber die Spitze des Peaks ist nur ein relatives und kein absolutes Maximum. Werte 
des LDR sind im Bereich der Schmelzschicht nur in Teilbereichen gemessen worden. Dennoch 
kann die Schmelzschicht mit dem nachfolgend präsentierten Algorithmus detektiert werden. 
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Abbildung 4.2: Profile aus RHI-Messungen des Wolkenradars MIRA36-S vom 04.09.2009, 

10:35 UTC folgender Größen: Reflektivität {Ciz)), lineares Depolarisationsverhältnis {L{z)) und 
Vertikalgeschwindigkeit {wm{z)). Grau unterlegt ist die delektierte Schmelzschicht. 
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4.2 Detektionsalgorithmus 

In jedem der drei Profile, die für einen idealen Fall in Abb. 4.1 gezeigt werden, sind deutlich 
detektierbare Merkmale zu erkennen. Diese Merkmale bilden die Gnmdlage des Detektionsal- 
gorithmus und werden im folgenden Abschnitt einzeln diskutiert, bevor auf die Kombination 
der Merkmale der drei Größen und die Bestimmung der Schmelzschichthöhe und -grenzen ein- 
gegangen wird. 



4.2.1 Eigenschaften der Profile als Ganzes 

Profil der Reflektivität Typisches Merkmal im gesamten Profil der Reflektivität ist ein loka- 
les Maximum (Peak) im Bereich der Schmelzschicht (Abb. 4.1). Die Analyse weiterer Profile 
der Reflektivität ergab auch: 

• Die stärkste konkave Krümmung ist oft im lokalen Maximum in der Schmelzschicht zu 
finden. Diese konkave Krümmung geht in den meisten Fällen zu jeder Seite in eine aus- 
geprägte konvexe Krümmung über. Dabei sind die stärksten konvexen Krümmungen in 
vielen Fällen am Anfang und Ende des Schmelzvorgangs zu finden. 

• Im Idealfall entspricht das lokale Maximum in der Schmelzschicht dem absoluten Maxi- 
mum des gesamten Profils. 

• Das lokale Maximum in der Schmelzschicht ist oft eingebettet in einen Trend der Werte 
in der Umgebung der Schmelzschicht (Abb. 4.2). 

• Die Anzahl der lokalen Maxima im gesamten Profil kann groß werden. 

Das augenscheinlich sicherste Merkmal, die Schmelzschicht zu identifizieren, ist die Abfolge 
starker Krümmungen (konvex, konkav, konvex, siehe Abb. 4.1 und 4.2). Deshalb ist ein geeig- 
netes Verfahren die Schmelzschicht zu detektieren, mit Hilfe der zweiten Ableitung der Profile 
der Reflektivitätswerte, ausgedrückt durch diskrete Messwerte 

C"(^) ^ (4.1) 

nach Maxima zu suchen. Die stärkste Krümmung im Profil der Reflektivität befindet sich meist 
innerhalb der Schmelzschicht und wird der Höhe der Schmelzschicht zugeschrieben. Ableitun- 
gen erhöhen das Rauschen der Daten und die Detektion der Schmelzschicht wird erschwert. 
Um dem entgegen zu wirken und möglichst glatte Profile zu erhalten, werden die Daten vor 
jeder Ableitung mit einem gleitenden Mittelwertsfilter über jeweils sieben Höhenstufen (70 m) 
geglättet. 
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Profil des LDR Ähnlich wie in den Profilen der Reflektivität ist auch im Profil des LDR ein 
lokales Maximum (Peak) im Bereich der Schmelzschicht (Abb. 4.1) zu erkennen. Die Untersu- 
chung vieler Profile ergab: 

• Die stärkste konkave Krümmung ist oft im lokalen Maximum in der Schmelzschicht zu 
finden. 

• Die konkave Krümmung des lokalen Maximum geht zu jeder Seite in eine ausgeprägte 
konvexe Krümmung über. Dabei sind die stärksten konvexen Krümmungen am Anfang 
und Ende des Schmelzvorgangs zu finden. 

• Es gibt Fälle, in denen keine oder nur wenige Werte in der Schmelzschicht vorhanden 
sind (Abb. 4.2). 

• Sind Werte an jedem Gitterpunkt innerhalb der Schmelzschicht vorhanden, ist das lokale 
Maximum in der Schmelzschicht meist das absolute Maximum im gesamten Profil. 

• Die Anzahl der lokalen Maxima im gesamten Profil kann groß werden. 

Auch für die Profile des LDR eignet sich die Abfolge starker Krümmungen als beständigstes 
Merkmal zur Detektion der Schmelzschicht. Wie auch für die Reflektivitäten wird die zweite 
Ableitung der Profile verwendet, um Maxima zu suchen: 

^"W « ^ (4.2) 

Die stärkste Krümmung befindet sich ebenfalls meist innerhalb der Schmelzschicht und wird 
der Höhe der Schmelzschicht zugeschrieben. Wie zuvor auf die Profile der Reflektivität ange- 
wendet, werden die Werte der LDR-Profile vor jeder Ableitung mit einem gleitenden Mittel- 
wertsfilter über jeweils sieben Höhenstufen (70 m) geglättet. 



Profil der Vertikalgeschwindigkeit Typisches Merkmal im Profil der Vertikalgeschwindig- 
keit ist eine deutliche Zunahme der Sedimentationsgeschwindigkeit der Streuer im Bereich der 
Schmelzschicht (Abb. 4.1). Die Sichtung einiger Profile der Vertikalgeschwindigkeit ergab: 

• Die stärkste Zunahme der Sedimentationsgeschwindigkeit ist meist in der Schmelzschicht 
zu finden. 

• Ein starker Gradient der Werte innerhalb der Schmelzschicht erstreckt sich über mehrere 
hundert Meter. 

• Starke Gradienten treten in manchen Fällen auch außerhalb der Schmelzschicht auf, er- 
strecken sich aber meist nur über wenige zehn Meter. 
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Zur Detektion der Schmelzschicht im Profil der Vertikalgeschwindigkeit ist der Gradient und 
somit die erste Ableitung das geeignete Mittel, um nach Maxima zu suchen: 

'i^M ~ (4.3) 

Der maximale Gradient von n',„{z) tritt meist in der Schmelzschicht auf. Wie zuvor werden 
auch die Werte der Profile der Vertikalgeschwindigkeit vor der Ableitung mit einem gleitenden 
Mittelwertsfilter über jeweils sieben Höhenstufen (70 m) geglättet. 



4.2.2 Kombination der Profile 

Mit der kombinierten Betrachtung der drei Größen (, L, Wm wird die Anzahl der Fehlde- 
tektionen reduziert und die Häufigkeit der richtigen Detektionen erhöht. Die Merkmale der 
Krümmung der Kurven der Reflektivitätsgrößen und die Merkmale der Steigung der Vertikalge- 
schwindigkeitskurve können zusammengeführt werden, wenn sie zuvor normiert wurden. Dies 

geschieht durch die Division mit definierten Grenzwerten (Gl. 4.5). Mit diesen Grenzwerten 
können die Messgrößen wie folgt kombiniert werden, wobei jeder Term das Gewicht 1/3 be- 
kommt^: 

S{z) = Usdz) + S,iz) + SUz)) = \ + ^ + ^) (4.4) 

Dabei ist S{z) als ein dimensionsloses Maß für das Vorhandensein einer Schmelzschicht, als 
Schmelzschichtindex, zu verstehen und soll für die Höhe der Schmelzschicht maximal werden. 
Die Krümmungen im Profil der Reflektivität und des LDR sind in der Höhe der Schmelzschicht 
konkav (negative Werte) und der Gradient im Profil der Vertikalgeschwindigkeit ist in der Höhe 
der Schmelzschicht positiv. Dies wurde bei der Wahl der Grenzwerte berücksichtigt, sodass 
S{z) in der Höhe der Schmelzschicht maximal wird. Diese Grenzwerte sollen festlegen, ab wel- 
chem jeweiligen Wert L'l und w'^ ^ die Schmelzschicht mit Hilfe von S{z) zuverlässig de- 
tektiert werden kann. Die Grenzwerte, die für den gesamten Datensatz verwendet werden, wur- 
den empirisch bestimmt: Anhand von Stichproben, für die die Schmelzschicht augenscheinlich 
war, wurden die Grenzwerte angepasst. Die beste Übereinstimmung der detektierten Schmelz- 
schichthöhen und -grenzen mit den per Auge definierten lieferten folgende Grenzwerte: 



Ct = -330 



dB 
kn? 



dB 
kn? 



m 



s km 



(4.5) 



Die hervorstechenden Merkmale im Vertikalprofil von S{z) sollten im idealen Fall im Bereich 
der Schmelzschicht erwartet werden: Ein Maximum innerhalb der Schmelzschicht und lokale 



'Der Index t für die Grenzwerte steht für threshold. 
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Minima an der Ober- und Untergrenze der Schmelzschicht (siehe Abb. 4.3, gelbe Linie zwi- 
schen 2.5 km und 3.5 km). Große Werte für S{z) treten allerdings mit Ausnahme des idealen 
Falls nicht nur im Bereich der Schmelzschicht auf. Dadurch sind mehrere Maxima innerhalb des 
Profils von S{z) möglich, die ähnliche Werte wie das Maximum innerhalb der Schmelzschicht 
annehmen können. Um das Maximum der Schmelzschicht zu detektieren, wird im ersten Schritt 
Eins als Schwellenwert für S{z) festgelegt. Somit ist eine Schmelzschicht vorhanden, wenn 
S{z) > 1. 




-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 
Schmelzschichtindex 



Abbildung 4.3: Vertikalprofile von S(^{z), Sl{z), S.u,{z) und deren gewichtetes Mittel S{z) als 
Funktion der Höhe z auf Basis der Messung vom 06.10.2009, 13:14 UTC. Grau unterlegt ist die 
Schmelzschicht. 

Im Idealfall ist damit die Detektion abgeschlossen. Es kann aber vorkommen, dass mehrere 
Maxima im 1. Schritt gefunden werden. Eine Reduktion der Anzahl der Maxima wird erreicht, 
indem die Form der Kurve von S{z) um das jeweilige Maximum analysiert wird. Dafür wird in 
jenen Höhen zq in denen S{z) > 1 ist, untersucht ob die Form der Kurve um 5(2:0), der Form 
entspricht, die im Bereich der Schmelzschicht zu erwarten wäre: Ein ausgeprägtes Maximum 
in der Mitte der Schmelzschicht, zwischen zwei ähnlich ausgeprägten Minima, die die Grenzen 
der Schmelzschicht bedeuten (Abb. 4.3). Diese Kurvenform lässt sich durch das sogenannte 
mexican hat-Profil beschreiben (Bronstein et al., 2008): 

F{z) = (1 - z^) ß-^ (4.6) 

mit Extremwerten bei z = { — 0, Zur ersten Abschätzung der Höhe und der Gren- 
zen der Schmelzschicht muss eine Höhe zq und eine Schmelzschichtdicke B in dem Sinne 
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gewählt werden, dass das Maximum in der Höhe zq liegt und die Oberkante und Unterkante der 
Schmelzschicht bei Zo±b (mit 6=1/2 B). Die Fitfunktion F{z) an die Werte von S{z) um die 
Höhe zq ergibt sich somit zu 



F{z,b)= 1 



3 (z - zoY 



_ 3 (»-^o) 

e 2 i.2 



(4.7) 



Für jede Höhe zq wird nun 6 in 10-Meter-Schritten von 30 m bis 500 m variiert. Für jene Höhe 
Zq und jene halbe Dicke b, für die der Korrelationskoeffizient nach Pearson 



QF,s{zo,b) = 



Goy{Fiz,b),S{z)) 
^V^T{F{z,b) VVar(5(z)) 



(4.8) 



maximal wird und größer 0.9 ist, ergibt sich die Abschätzung der Schmelzschichthöhe zu Zp = 
Zq und der Schmelzschichtdicke zu i? = 2b. Ist die Kurve von S{z) in der Schmelzschicht 
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Abbildung 4.4: Ausschnitt aus Abb. 4.3, nur Bereich der Schmelzschicht. Die rote Linie ist die 
Fitfunl<tion mit zq = 2.91 km und 6 = 170 m für die RHI-IVIessungen vom 06.10.2009, 

13:1 4 UTC (linl<s) und mit = 2.96 l<m und 6 = 230 m für die RHI-IVIessungen vom 06.10.2009, 
14:28 UTC (rechts). Grau unterlegt ist die delektierte Schmelzschicht. 

achsensymmetrisch zum Maximum in der Höhe zq, so kann die Oberkante (zt = zo + b) und 
Unterkante (Zb = Zq — b), sowie die Dicke (B = 2 b) der Schmelzschicht mit Hilfe des mexican 
hat-Fits direkt berechnet werden (Abb. 4.4, links). Die Indizes p, t und b stehen jeweils für 
peak, top und bottom. In einigen Fällen ist die Achsensymmetrie nicht gewährleistet und die 
Orte (Höhen) des Maximums und der Minima der Fitfunktion in der Schmelzschicht stimmen 
nicht genau mit denen von S{z) überein (Abb. 4.4, rechts). Eine Verbesserung wird mit einem 
zweiten Fit erreicht. Nachdem die Höhen der Minima (zt, z^) und des Maximums {Zp) ungefähr 
bekannt sind, wird der Fitbereich in drei Segmente um den jeweiligen Extremwert eingeteilt 
und mit dem passenden Ausschnitt der mexican hat-Funktion gefittet. Für die Bestimmung der 
Höhe der Oberkante der Schmelzschicht wird für alle z E — 100 m, + 100 m] mit Hilfe der 
Gleichungen 4.7 und 4.8 (^o wird dabei durch Zt ersetzt) der beste Fit für den Höhenbereich um 
die Oberkante der Schmelzschicht und somit ein neues Zt und b berechnet. Gleichermaßen wird 
für den Bereich um die Unterkante (zq wird durch Z}, ersetzt und z G [zh— 100 m, z^ + 100 m]) 
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Abbildung 4.5: Beste Fits der mexican hat-Funktion F(z,6) (rot) an die Kurve S{z) (gelb) 
und der beste Fit F{zt,b), F{zb,b), F{zt,b) (schwarz gepunl<tet) jeweils im Bereich der Minima 
und des Maximums an die Kurve S{z) für RHI-Messungen vom 06.10.2009, 13:14 UTC. Grau 
unterlegt ist die delektierte Schmelzschicht. 




Abbildung 4.6: wie Abb. 4.5, für RHI-Messungen vom 06.10.2009, 14:28 UTC. 
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und den Bereich um das Maximum (2:0 wird durch Zp ersetzt und z E[zp— 100 m, Zp + 100 m]) 
verfahren. 

Der zweite Fit an drei unabhängige Höhenbereiche von S{z) bestätigt im symmetrischen Fall 
den ersten Fit (Abb. 4.5), während im asymmetrischen Fall die Detektion von Zu Zb und Zp 
verbessert wird (Abb. 4.6). Diese Zt und Zb sind die Grundlage für die Berechnung der Dicke 
der in allen Abbildungen grau unterlegten Schmelzschicht. 

4.3 Eigenschaften der Schmelzschicht 

In diesem Unterkapitel wird gezeigt, was das Wolkenradar in der Schmelzschicht und in der 
Nähe der Schmelzschicht messen kann, wie diese Messungen aus meteorologischer Sicht inter- 
pretiert werden können und wie sich Messungen und Annahmen mit Literaturwerten decken, 
bzw. ergänzende Informationen liefern. 

Zunächst wird auf die detektierten Schmelzschichtdicken eingegangen und darauf, wodurch 
diese bestimmt werden. Dabei werden die zur Detektion verwendeten Messgrößen direkt ober- 
halb der Schmelzschicht mit der Dicke der Schmelzschicht in Relation gesetzt. Insgesamt ste- 
hen 3781 RHl-Scans zur Verfügung, für die eine Schmelzschicht detektiert wurde. Die Höhe 
direkt oberhalb der Schmelzschicht wurde auf Zt + 100 m festgelegt. Diese Distanz ist einer- 
seits weit genug von der Schmelzschicht entfernt, um deren Einfiuss zu minimieren und liegt 
andererseits ausreichend nahe an der Schmelzschicht, um Aussagen über die Eigenschaften der 
Eispartikel, bevor sie schmelzen, treffen zu können. Aus der Literatur ist bekannt, dass die 
Dicke der Schmelzschicht durch die Schmelzeigenschaften der jeweiligen Eispartikel bestimmt 
wird, hauptsächlich durch die Größe, Form (Verhältnis von Volumen zu Oberfläche) und Dichte 
(Lufteinschlüsse) der Eispartikel an der Obergrenze der Schmelzschicht. 

Die Reflektivität ist stark abhängig von der Größe der schmelzenden Eispartikel (D^- 
Abhängigkeit in der Rayleigh-Approximation, siehe dazu Gl. 2.7). Es ist zu erwarten, dass 
größere Eispartikel eine dickere Schmelzschicht verursachen als kleinere. Die Reflektivität soll- 
te deutlich mit der Schmelzschichtdicke korrelieren. Dies zeigt sich anhand Abb. 4.7. In dieser 
Abbildung, wie in allen Streudiagrammen in diesem Kapitel, ist die Anzahl der Werte loga- 
rithmisch aufgetragen, wobei die Grauskalen die Zahlenwerte angeben. Dadurch wurde jede 
einzelne Messung oder Detektion in der Darstellung berücksichtigt. Um Zusammenhänge bes- 
ser erkennen zu können, ist in Abb. 4.7 das 20., 50. und das 80. Perzentil der Reflektivität direkt 
oberhalb der Schmelzschicht als Funktion der Schmelzschichtdicke aufgetragen (rote Linien). 
In den nachfolgenden Abbildungen 4.8, 4.10 und 4.11 sind die Perzentile der jeweiligen Größe 
ebenfalls als Funktion der Schmelzschichtdicke rot eingezeichnet. Anhand des 50. Perzentils 
in Abb. 4.7 ist für Reflektivitäten geringer als OdBz zu erkennen, dass die Schmelzschicht- 
dicke im Mittel um ca. 30 m zunimmt, wenn sich die Reflektivität um 5 dB erhöht. Oberhalb 
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Abbildung 4.7: Streudiagramme der Schmelzschichtdicke ß in m als Funktion der Reflektivität 
C in dBz direkt oberhalb der Schmelzschicht {zt + 100 m). Logarithmische Grautondarstellung 
der Anzahl von Detektionen. Die roten Linien kennzeichnen das 20., 50. und 80. Perzentil. 

von 0 dBz ist mit Zunahme der Schmelzschichtdicke eine abnehmende Abhängigkeit mit der 
Reflektivität festzustellen. Das 50. Perzentil der Reflektivität läuft bei zunehmender Schmelz- 
schichtdicke gegen einen Grenzwert von etwa 5 dBz. Die Vermutung liegt nahe, dass die Partikel 
oberhalb der Schmelzschicht mit zunehmender Schmelzschichtdicke größer sind, deren Durch- 
messer aber vermehrt außerhalb der Rayleigh-Annahme (A ^ D) liegen und dadurch vom 
Radar unterschätzt werden. Diese Vermutung wird durch Literaturdaten bestätigt: Anhand flug- 
zeuggestützter /n-^ito-Messungen konnte gezeigt werden, dass einzelne Eispartikel an der Ober- 
kante der Schmelzschicht Größen von bis zu 12 mm (Stewart et al., 1984) bzw. 14.5 mm (Evans 
et al., 2005) erreichen. 

Abhängig von der Radarwellenlänge erreicht der Rückstreuquerschnitt und damit zusam- 
menhängend die Reflektivität mit zunehmendem Partikeldurchmesser beim Übergang vom 
Rayleigh- zum Miebereich ein Maximum und sinkt für weiter wachsende Partikel im Miebe- 
reich auf ein Minimum, bevor der Rückstreuquerschnitt wieder ansteigt (Resonanzeffekt, vgl. 
Abb. 2.1). Das erste Minimum, auch Dark Band genannt, wurde bereits bei Messungen mit Ra- 
daren größerer Wellenlängen direkt oberhalb des Bright Band beobachtet (Fabry und Zawadzki, 
1995; Sassen et al., 2005). Für die Wellenlänge des Wolkenradars ergibt sich näherungsweise, 
dass ab einem Durchmesser einer vergleichbaren homogenen, sphärischen Eiskugel von et- 
wa 2.7 mm (vgl. Abb. 2.1) die Werte, die durch Rayleigh-Näherung ermittelt wurden, deutlich 
von denen abweichen, die mittels der Mie-Streuung berechnet wurden und ein Maximum errei- 
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Abbildung 4.8: Wie Abb. 4.7, für das Streudiagramm der Schmelzschichtdicke ß in m als 
Funktion der Vertikalgeschwindigkeit in m/s direkt oberhalb der Schmelzschicht (24+ 100 m). 
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Abbildung 4.9: Streudiagramm der Reflektivität ({zt + 100 m) in 6Bz als Funktion der Verti- 
kalgeschwindigkeit Wmizt + 100 m) in m/s. Logarithmische Grautondarstellung der Anzahl von 
Detektionen. 
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chen. Dies führt zu der Hypothese, dass oberhalb der Schmelzschicht bei Schmelzschichtdicken 
größer als 300 m zunehmend Partikel mit Durchmessern größer als D ~ 3 mm vorliegen. Die 
Sedimentationsgeschwindigkeit der Eispartikel, die für Reflektivitäten im Bereich von etwa 
5dBz verantwortlich sind, liegen zwischen 1m/s und 2.5 m/s (Abb. 4.9). Es ist davon auszu- 
gehen, dass die Eispartikel, die für Schmelzschichtdicken größer als 300 m verantwortlich sind, 
eher von bereiften Eispartikeln, wie z.B. Graupel, verursacht werden und weniger von regulären 
Schneeflocken. Während unbereifte, mehrere Millimeter große Eispartikel (Schneeflocken) ty- 
pische Sedimentationsgeschwindigkeiten von etwa 1 m/s annehmen, kann Graupel auch Sedi- 
mentationsgeschwindigkeiten von deutlich über 2 m/s erreichen (siehe auch Pruppacher und 
Klett, 1980). 

Die Schmelzschichtdicke wächst mit zunehmender Sedimentationsgeschwindigkeit der Eispar- 
tikel oberhalb der Schmelzschicht (Abb. 4.8). Für Schmelzschichtdicken < 200 m führt eine 
Zunahme der Sedimentationsgeschwindigkeit der Eispartikel um 10 cm/s zu einer Zunahme der 
Schmelzschichtdicke von ca. 13 m. Schmelzschichtdicken > 200 m sind bei einer Zunahme der 
Sedimentationsgeschwindigkeit von 10 cm/s um ca. 44 m dicker. Eine mögliche Erklärung für 
die verringerte Zunahme der Sedimentationsgeschwindigkeit ist, dass mit zunehmender Größe 
zwar die Masse der Partikel zunimmt, aber die Oberfläche der Partikel deutlicher ansteigt und 
damit auch deren Reibungswiderstand. Wenn man annimmt, dass die Eispartikel in diesem 
Zusammenhang zunehmend komplexere Strukturen aufweisen und dabei von der sphärischen 
Form abweichen, müssten auch die Werte des LDR mit zunehmender Dicke der Schmelzschicht 
ansteigen. Der Zusammenhang ist schwach sichtbar, wenn nur das 50. Perzentil betrachtet wird 
(Abb. 4.10). Die große Streuung der Werte lässt allerdings keinen klar sichtbaren Zusammen- 
hang erkennen. 

Der näherungsweise lineare Zusammenhang der Sedimentationsgeschwindigkeit der Eispartikel 
mit der Dicke der Schmelzschicht (Abb. 4.8) legt die Vermutung nahe, dass die Schmelzdauer 
der Eispartikel unabhängig von ihrer Größe sein könnte. Um dies zu überprüfen, wurde die 
Schmelzdauer (Der Index m steht für melting) berechnet: 



Der Berechnung liegt die Annahme zugrunde, dass das Profil der Sedimentationsgeschwindig- 
keit der Partikel durch die Schmelzschicht dem Profil der Vertikalgeschwindigkeit (Wm) des 
jeweiligen RHI-Scans entspricht, der auch zur Abschätzung der Schmelzschichtdicke verwen- 
det wurde. 

Eine Darstellung der Schmelzschichtdicke als Funktion der Schmelzdauer findet sich in Abb. 
4.11. Für Schmelzschichtdicken zwischen 200 m und 500 m trifft die oben angestellte Ver- 
mutung zu. Die Schmelzdauer der Eispartikel ist annähernd konstant und liegt bei ca. 150 s 
(50. Perzentil). Der Effekt der Abkühlung der Luft aufgrund des Phasenübergangs, wodurch. 




(4.9) 
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Abbildung 4.10: wie Abb. 4.7 für die Streudiagramme der Schmelzschichtdicke als Funktion 
des LDR L in dBdirekt oberhalb der Schmelzschicht. 

abhängig vom Temperaturgradienten bevor Niederschlag einsetzt, Isothermie innerhalb der 
Schmelzschicht erzeugt werden könnte, ist nur für Schmelzschichtdicken unterhalb 200 m zu 
beobachten. Eine konstante Schmelzdauer würde dadurch möglich, dass die Temperaturzunah- 
me der Umgebung der Eispartikel während des Fallens durch die Schmelzschicht den Schmelz- 
vorgang beschleunigt und somit große Eispartikel näherungsweise genau so schnell vollständig 
geschmolzen sind wie kleinere. Gestützt wird diese Theorie dadurch, dass isotherme Luft- 
schichten bisher überwiegend im Bereich des Oberrands der Schmelzschicht beobachtet wur- 
den, während die Temperatur innerhalb der Schmelzschicht meist von oben nach unten zunimmt 
(Stewart et al., 1984). 

Die Größe der Eispartikel an der Oberkante der Schmelzschicht wird bestimmt durch Wachs- 
tumsprozesse oberhalb der Schmelzschicht. Vorwiegend im Bereich zwischen -10 °C und 
0°C findet das stärkste Wachstum der Eispartikel statt (Evans et al, 2005), dabei nimmt die 
Reflektivität zu (Fabry und Zawadzki, 1995). Der überwiegende Wachstumsprozess wird dabei 
der Aggregation der Eispartikel zugeschrieben. In dieser Schicht, die in den häufigsten Fällen 
zwischen -8 °C und -3 °C liegt, ist meist auch das Maximum der Anzahl der Eispartikel zu 
finden, deren Zahl in Richtung der Schmelzschicht wieder abnimmt, während die Größe der 
Eispartikel zunimmt (Stewart et al., 1984; Evans et al., 2005; Woods et al, 2008). Dabei wird 
vermutet, dass die maximale Anzahl der Eispartikel erreicht wird, indem neue Eispartikel durch 
sekundäre Eisproduktion in Form von Nadeln entstehen, die durch das Absplittern von Eispar- 
tikeln während des Bereifens entstanden sind (Hallett und Mossop, 1974; Woods et al., 2008). 
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Abbildung 4.11: wie Abb. 4.7 für die Streudiagramme der Schmelzschichtdicke als Funktion 
der Schmelzdauer tm in s. 

Diese Wachstumsprozesse sind andeutungsweise im Median der Profile in Abb. 4.12 zu er- 
kennen. Gezeigt sind Abweichungen der Werte der Reflektivität, des LDR und der Vertikalge- 
schwindigkeit vom Wert der jeweiligen Größe in der Schmelzschichthöhe Zp. Die Höhen sind 
ebenfalls als Abweichung von der Höhe Zp dargestellt. Das 50. Perzentil, der Median, spiegelt 
im Wesentlichen den in Kap. 4 beschriebenen idealen Fall zur Detektion einer Schmelzschicht 
wieder. Dafür wird angenommen, dass das Medianprofil der Einzelprofile die Eigenschaften 
der Einzelprofile wiedergibt. Oberhalb von 2.5 km nimmt die Reflektivität von oben nach un- 
ten zu, bleibt zwischen 2.5 km und 2 km fast konstant und wächst dann mit ca. 3 dB/km bis 
zum Oberrand der Schmelzschicht. Gleichzeitig weist die Sedimentationsgeschwindigkeit der 
Eispartikel ab etwa 2 km oberhalb der Schmelzschicht eine Zunahme auf. Der zuvor beschrie- 
bene Anstieg der Anzahl, gefolgt von der Größenzunahme der Eispartikel durch Aggregation 
unterhalb der — 10 °C-Isotherme wird dadurch angedeutet. Mit der Annahme eines feuchta- 
diabatischen Temperaturgradienten von 0.5 K/100 m ist eine Temperatur von — 10°C in einer 
Höhe von 2 km oberhalb der Schmelzschichtoberkante zu finden. Interessant ist zudem, dass 
die Zunahme der Reflektivität unterhalb 2 km oberhalb der Schmelzschichthöhe für niedrige 
Reflektivitäten (20. Perzentil) deutlicher wird. Für höhere Reflektivitäten (80. Perzentil) ist dies 
nicht mehr zu erkennen. Dort spielen zunehmend konvektive Prozesse eine Rolle, die die Ein- 
zelprofile beeinflussen. Ein solches Profil der Reflektivität ist in Abb. 4.2 gezeigt. 

Nachdem die Messgrößen in Bezug zu Prozessen oberhalb der Schmelzschicht und deren Ein- 
fluss auf die Schmelzschicht gesetzt wurden, folgt nun die Diskussion der Messgrößen in Bezug 
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auf die Vorgänge innerhalb der Schmelzschicht. Aus der Literatur ist bekannt, dass zu Beginn 
des Schmelzprozesses ebenfalls Aggregation stattfindet, die zu einer weiteren Massezunahme 
der einzelnen Partikel führt. Weitere Prozesse, wie die Bereifung und die Kondensation, bzw. 
Deposition, tragen ebenfalls bei zum Massenwachstum und zur Erhöhung der Niederschlagsra- 
te, teilweise Verdopplung (Evans et al., 2005). Dass zudem kleine Wolkenpartikel direkt ober- 
halb und innerhalb der Schmelzschicht vorhanden sein können, wurde durch Melchionna et al. 
(2008) gezeigt. 



Das Verhalten der drei Messgrößen innerhalb der detektierten Schmelzschichten wird anhand 
der Ableitungen der 50. Perzentile der Messgrößen gezeigt. Dafür werden, wie in Kap. 4.2 
beschrieben, aus den Profilen der 50. Perzentile die Profile S(;{z), Sl{z), Sw{z) und S{z) be- 
rechnet (Abb. 4.13) und die Schmelzschichthöhe, -Oberkante und -unterkante detektiert. Das 
Ergebnis ist in Abb. 4.13 dargestellt. Die Höhe ist als Differenz zur Schmelzschichthöhe Zp auf- 
getragen. Interessant ist, dass die Höhe der stärksten konvexen Krümmungen der Reflektivität 
(bzw. der Minima von 5^), die die Grenzen des BrightBand andeuten, von denen des LDR (bzw. 
der Minima von Sl) und von den Grenzen des stärksten Anstiegs der Vertikalgeschwindigkeit 
abweichen. Diese Erkenntnis war nicht in der gesichteten Literatur zu finden. 
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Abbildung 4.12: 20.(rot), 50.(schwarz) und 80. (grün) Perzentil der Profile der Reflekti- 
vitätsdifferenz ({z) - ({zp) in dB, der Differenz des LDRs L{z) - L{zp) in dB und der Differenz 
der Vertil<algesclnwindigkeit w,„(z) - Wm{zp) in m/s als Funktion der Höhe z - Zp\r\ m. Grau 
unterlegt ist die delektierte Schmelzschicht. 
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Die automatische Detektion einer Schmelzschicht, bzw. deren Grenzen Zf,, Zt und deren Höhe Zp 
ist erfolgreich, wenn alle drei Größen verwendet werden. Deshalb wird für die mittleren Eigen- 
schaften der drei Größen in Bezug auf die Schmelzschichtgrenzen das Profil von S{z) auf Basis 
der Medianprofile berechnet. Auf Basis des Profils von S{z) wird die Oberkante der Schmelz- 
schicht durch die Reflektivität und die Vertikalgeschwindigkeit und die Unterkante vorwiegend 
durch das LDR und die Vertikalgeschwindigkeit bestimmt (Abb. 4.13). In der zitierten Literatur 
lag das Hauptaugenmerk, aufgrund der meist deutlich gröberen Auflösung der Messung, auf der 
Detektion der Schmelzschicht. Die unterschiedliche Reaktion der Messgrößen in Bezug auf die 
Schmelzschichtgrenzen wurde nicht analysiert. Eine detaillierte Betrachtung darüber, welche 
Messgröße in welcher Höhe ein Signal zur Detektion der Schmelzschichtgrenzen liefert, wurde 
nicht durchgeführt. Das in dieser Arbeit gezeigte Ergebnis (Abb. 4.13) ermöglicht erstmals die 
Angabe von Zahlenwerten für den Unterschied der delektierten Schmelzschichtoberkanten und 
-unterkanten verschiedener Messgrößen. 

In Abb. 4.13 ist zu erkennen, dass oberhalb der Schmelzschichthöhe Zp das Minimum des 
Schmelzschichtindex auf Basis der Reflektivität (S(;{z)) und der Anstieg des Schmelzschichtin- 
dex auf Basis der Vertikalgeschwindigkeit (Sw{z)) knapp oberhalb der Oberkante der Schmelz- 
schicht zu finden ist. Etwa 40 m darunter, bereits innerhalb der detektierten Schmelzschicht ist 
das Minimum des Schmelzschichtindex auf Basis des LDR (Sl{z)) zu sehen. Der Vergleich der 
Maxima ergibt, dass das Maximum des Profils von S(;{z) ca. 20 m oberhalb und die Maxima 
von SfX^)^ Si,,{z) ca. 20m unterhalb des Maximums des Schmelzschichtindex S{z) liegen. 
Unterhalb der Schmelzschichthöhe Zp ist das Minimum des Profils von Sl{z) an der Schmelz- 
schichtuntcrkantc, der Übergang zu näherungsweisc konstanten Werten für Su,{z) direkt darun- 
ter und das Minimum des Profils von 5*^(2:) ca. 50 m darüber zu erkennen. Die Oberkante der 
Schmelzschicht wird im Wesentlichen durch Änderungen in der Reflektivität und in den Verti- 
kalgeschwindigkeiten bestimmt, während die Unterkante durch Änderungen in den Werten des 
LDR und der Vertikalgeschwindigkeit bestimmt werden. 

Zur Veranschaulichung der Vorgänge innerhalb der Schmelzschicht, in Bezug auf die detek- 
tierte Schmelzschicht aus den Medianprofilen, werden aus der Literatur bekannte Vorgänge in 
den Zusammenhang der gezeigten Profile (Abb. 4.13) zur Detektion der Schmelzschichtgrenzen 
gebracht. Eine Hypothese wird aufgestellt, welche die Reihenfolge und Relevanz der wesentli- 
chen Prozesse und Auswirkungen der Prozesse, die in der detektierten Schmelzschicht während 
des Sedimentierens der Partikel zum Tragen konmien, spezifizieren und die gezeigten Profile 
erklären soll: 

1. Die Eispartikel (Aggregate, hauptsächfich Schneeflocken und Graupel) beginnen an den 
Oberflächen zu schmelzen. Dadurch wird die Aggregation effizient, da die Hydrometeore 
nun effektiver zusammenhaften können. Die mittlere Größe der Hydrometeore nimmt 
dabei weiter zu, die Reflektivität und die Sedimentationsgeschwindigkeit steigen an. 
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Abbildung 4.13: Profile relativ zur Schmelzschichthöhe Zp von Sc,{z), Sl{z), S^{z) auf Basis 
des 50. Perzentils aller Profile, für die eine Schmelzschicht delektiert wurde. Grau unterlegt ist 
die daraus delektierte Schmelzschicht. Beschriftung 1. bis 6. siehe Text. 



2. Durch die Aggregation und Bereifung von Hydrometeoren wachsen diese weiter an und 
die Form der Hydrometeore ändert sich hin zu komplexeren, asphärischen Partikeln und 
eine flüssige Oberfläche bildet sich um die Eispartikel. Das LDR steigt an. Gleichzeitig 
nimmt der Dielektrizitätsfaktor mit zunehmendem Wassergehalt der Hydrometeore ab 
und die Reflektivität steigt weiter an (Größeneffekt). 

3. Da die Aggregation anhält, nimmt die Sedimentationsgeschwindigkeit der Hydrometeo- 
re weiter zu und die Anzahl der Hydrometeore pro Volumeneinheit ab. Dies wirkt der 
Zunahme der Reflektivität entgegen. Ist die Reflektivität maximal, wirkt dieser Prozess 
genau der Zunahme des Dielektrizitätsfaktors und dem Wachstum der Hydrometeore ent- 
gegen. 

4. Im weiteren Schmelzen der Hydrometeore fließen zunehmend flüssige Partikel zusammen 
(Koagulation), und ehemals nichtsphärische Partikel beginnen in sphärische überzugehen. 
Das Maximum des LDR ist dabei erreicht, wenn dieser Prozess genau der Aggregati- 
on von noch vorhandenen Eispartikeln zu komplexeren Partikelformen entgegen wirkt. 
Durch die Annäherung der Form der Hydrometeore an eine Kugel sinkt auch der 
Strömungswiderstand der Hydrometeore und die höchste Geschwindigkeitszunahme der 
Hydrometeore wird erreicht. 

5. Die Geschwindigkeitszunahme wird dadurch reduziert, dass immer weniger langsam fal- 
lende Eispartikel vorhanden sind. Dadurch steigen die Abstände der Partikel zueinander 
kaum noch und die Anzahl der Partikel pro Volumeneinheit nimmt noch minimal ab. 
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Dieser Abnahme steht die Zunahme des Dielektrizitätsfaktors entgegen und führt dazu, 
dass die Reflektivität kaum noch sinkt. Die Unterkante der Schmelzschicht, die durch ei- 
ne starke konkave Krümmung der Reflektivitätswerte (als Funktion der Höhe) bestimmt 
wird, ist erreicht, der Schmelzprozess aber noch nicht abgeschlossen. 

6. Sind alle Eispartikel vollständig geschmolzen, liegen die Hydrometeore am Ende des 
Schmelzprozesses näherungsweise sphärisch vor, das LDR erreicht niedrige Werte und 
die Geschwindigkeitszunahme der Hydrometeore endet. Es liegen nur noch Wassertrop- 
fen vor. Diese können koagulieren und größere Regentropfen bilden, aber auch durch 
Stoßprozesse zerfallen (Breakup). Anhand Abb. 4.12 ist zu erkennen, dass sich die Wer- 
te aller Größen unterhalb der Schmelzschicht nur noch geringfügig ändern. Bei ge- 
nauerer Betrachtung ist allerdings unterhalb der Schmelzschicht, bis zu einer Höhe von 
1 km unterhalb der Unterkante, eine leichte Zunahme der Sedimentationsgeschwindig- 
keit zusammen mit einer geringen Zunahme der Reflektivität erkennbar. Daraus lässt sich 
schließen, dass die Tropfen weiter wachsen und die Zunahme der Tropfengröße durch 
Koagulation etwas effektiver abläuft als der Breakup-Prozess. 



4.4 Zusammenfassung 

Die Schmelzschicht wurde bisher meist mit einer Kombination von Messgrößen aus Volumen- 
scans mit Niederschlagsradaren oder mit einem zeitlichen Mittel einzelner Messgrößen vertikal 
zeigender Wolkenradare detektiert. In dieser Arbeit wurde erstmals ein Algorithmus vorge- 
stellt, der mit der Kombination dreier Größen arbeitet, die aus Messungen eines Scan-fähigen 
Wolkenradars abgeleitet wurden. Dieser Algorithmus ermöglicht die Berechnung der Höhe der 
Schmelzschicht, sowie die Ober- und Unterkante der Schmelzschicht. Die zur Detektion ver- 
wendeten Größen sind die Reflektivität, das lineare Depolarisationsverhältnis und die Vertikal- 
geschwindigkeit der Streuer. Der Vorteil dieses Detektionsverfahrens ist, dass für Messungen, 
die der Untersuchung der Entwicklung der Wolke dienen, auch die Schmelzschicht mit deren 
Dicke detektiert werden kann, ohne das Scan- Verfahren ändern zu müssen. Messungen einer 
vertikal ausgerichteten Antenne zwischen den RHI-Scans sind dadurch nicht notwendig, um 
die Schmelzschicht zu detektieren. Ein weiterer Vorteil gegenüber bisher angewandten Detek- 
tionsalgorithmen ist die große Höhenauflösung der Schmelzschichtgrenzen. Die Resultate der 
statistische Auswertung der detektierten Schmelzschichtgrenzen sind nachfolgend aufgeführt. 

Anhand der ermittelten Schmelzschichtgrenzen konnten im Messzeitraum Schmelzschicht- 
dicken von 150 m bis 500 m beobachtet werden. Die detektierten Schmelzschichten wurden 
auch auf Zusammenhänge der drei Größen mit der Schmelzschichtdicke untersucht. Dabei 
wurde gezeigt, dass für Reflektivitäten kleiner OdBz (entspricht einer Schmelzschichtdicke 
von etwa 300 m) bei einer Zunahme der Reflektivität direkt oberhalb der Schmelzschicht um 
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ca. 5 dB die Schmelzschichtdicke um 30 m größer ist. Für Reflektivitäten größer 0 dBz ist je- 
doch eine Abweichung dieses Zusammenhangs zu erkennen, der bereits in der Literatur be- 
schrieben wurde. Betrachtet man den Median der Schmelzschichtdicke als Funktion der Re- 
flektivitäten, nimmt die Schmelzschichtdicke für Reflektivitäten größer 0 dBz weiter zu, von 
300 m auf knapp über 500 m , während die Reflektivitäten gegen einen Grenzwert von etwa 
5 dBz laufen. Dies Hegt in der Größe der Partikel begründet. Die Rayleighannahme ist nicht 
mehr gültig, da sich der Durchmesser der Eispartikel der Radarwellenlänge annähert. 

Anhand der gezeigten Messwerte und theoretischen Werte für den Rückstreuquerschnitt wur- 
de gezeigt, dass der mittlere Durchmesser der Partikel direkt oberhalb der Schmelzschicht für 
Schmelzschichtdicken von mehr als 300 m mehr als « 3 mm beträgt. 

Nimmt die Sedimentationsgeschwindigkeit der Streuer direkt oberhalb der Schmelzschicht um 
ca. 9 cm/s zu, erhöht sich die Schmelzschichtdicke ebenfalls um 30 m. Dies gilt für Schmelz- 
schichtdicken größer als etwa 200 m. Daraus geht hervor, dass die Eispartikel mit zunehmender 
Größe eine dickere Schmelzschicht verursachen. 

Bei Betrachtung der Mediane fällt auf, dass für Partikel mit Reflektivitäten > —10 dBz und 
Sedimentationsgeschwindigkeiten > 1 m/s die Schmelzdauer konstant ist. Der Median der 
Schmelzschichtdicke als Funktion der Schmelzdauer zeigt dabei konstante Werte von 150 s für 
Schmelzschichtdicken zwischen 200 m und 500 m. Das bedeutet, dass ein Ensemble von Eispar- 
tikeln unabhängig ihrer mittleren Größe und Form die gleiche Zeit zum Schmelzen benötigen 
würde. Daraus lässt sich schließen, dass ein Ensemble von Hydrometeoren, die aufgrund ih- 
rer größeren Masse im Laufe des Schmelzvorgangs größere Sedimentationsgeschwindigkeiten 
annehmen als ein Ensemble von Hydrometeoren geringerer Masse, in wärmere Luftschichten 
fallen und der Schmelzprozess in dem Maße beschleunigt wird, dass die Schmelzdauer kon- 
stant bleibt. Für geringere Reflektivitäten (< —10 dBz ) und Sedimentationsgeschwindigkeiten 
(< 1 m/s ) nimmt die Schmelzdauer mit der Dicke der Schmelzschicht zu. Vermutlich ist der 
Effekt der Abkühlung der Luft aufgrund des Phasenübergangs von Eis-Partikeln zu Wasser- 
Partikeln dominierend und innerhalb der Schmelzschicht herrscht weitgehend Isothermie vor. 
Hierfür müssten, im Unterschied zu Schmelzschichtdicken innerhalb der die Schmelzdauer kon- 
stant ist, eine große Anzahl kleinerer Hydrometeore vorliegen, die einen großen Beitrag zur 
Abkühlung leisten. Durch den geringeren Durchmesser ist zudem der Weg während des Schmel- 
zens geringer und die Schmelzdauer ist abhängig vom Durchmesser des Partikel-Ensembles. 

Zum Schluss wurden die Median-Profile, sowie Profile des 20. und 80. Perzentils der drei 
Messgrößen diskutiert. Die Mediane zeigen andeutungsweise ein Wachstum von Eiskristallen 
oberhalb der Schmelzschicht. Für niedrigere Reflektivitäten (unterhalb des Medians) ist dies 
deutlicher zu erkennen. Zur Schmelzschicht hin nehmen die Reflektivität und die Sedimenta- 
tionsgeschwindigkeit der Streuer zu. Die Größe der Eispartikel nimmt zu, Schneeflocken und 
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Graupelp artUcel werden gebildet. In der Schmelzschicht sind typische lokale Maxima der Re- 
flektivität und des LDR, sowie eine Zunahme der Sedimentationsgeschwindigkeit zu erkennen. 
Die hohe räumliche Auflösung der Messung erlaubt eine differenziertere Betrachtung der Ei- 
genschaften der drei Messgrößen innerhalb der Schmelzschicht. Die Schmelzschichtgrenzen 
der jeweiligen Größe weichen dabei ab von der klassischen Vorstellung (z.B. Houze, 1993) und 
Unterschiede in der Höhe der detektierten Grenzen sind erkennbar und können abgeschätzt wer- 
den. So wird die Oberkante der Schmelzschicht bestimmt durch ausgeprägte Änderungen der 
Sedimentationsgeschwindigkeit und der Reflektivität, die beide ansteigen. Die schmelzschicht- 
bedingten Änderungen der Werte des LDR erfolgen ca. 40 m unterhalb. Die Unterkante der 
Schmelzschicht wird dagegen bestimmt durch Änderungen der LDR- Werte und der Sedimen- 
tationsgeschwindigkeit. Die Änderungen der Reflektivität ist ca. 50 m oberhalb zu erkennen. 

Das LDR zeigt die deutlichsten Signale innerhalb der Schmelzschicht und ist der beste Indikator 
für das Vorhandensein einer Schmelzschicht. Da das LDR allerdings nicht immer zur Verfügung 
steht, ist die Detektion mit kombinierten Messgrößen erfolgreicher. Für die möglichst genaue 
Bestimmung der Schmelzschichtgrenzen ist die kombinierte Betrachtung der drei Größen wie- 
derum von Vorteil, da diese abhängig vom Stadium des Schmelzprozesses zu verschiedenen 
Zeitpunkten des Schmelzvorgangs reagieren. 

Der in diesem Kapitel vorgestellte Algorithmus in Verbindung mit räumlich hoch aufgelösten 
Radar-Messungen eines Scan-fähigen Radars ermöglicht eine detailliertere Betrachtung der 
Schmelzschicht als dies bisher mit Hilfe von Radardaten geschehen ist. 
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Die Vielfältigkeit der Größe und Form von Wolken spiegelt die Variationen der dynamischen 
Prozesse wieder, die für die Entstellung von Wolken und deren weiterer Entwicklung verant- 
wortlich sind. In diesem Kapitel werden Verfahren gezeigt, aus Radardaten Größen abzuleiten, 
die dynamische Vorgänge in Wolken aufzeigen. Im Zuge dessen wird eine neue Größe ein- 
geführt, die residuale Geschwindigkeit, die bisher in der Radarmeteorologie noch nicht ver- 
wendet wurde. Am Ende des Kapitels schließt sich die Interpretation der durch die Verfahren 
sichtbar gemachten dynamischen Prozesse in Wolken an. 

Um dynamische Prozesse in Wolken erkennen und verfolgen zu können, sind zeitlich hoch 
aufgelöste zweidimensionale Messungen der Streuergeschwindigkeiten von Vorteil. Zu diesem 
Zweck wurden Dopplermessungen während RHI-Scans in die gleiche Scanrichtung im zeitli- 
chen Abstand von 30 Sekunden durchgeführt (Kap. 3.4.3). Die dabei aufgezeichneten radialen 
Geschwindigkeiten können, bei Überschreiten der Nyquistgeschwindigkeit, gefaltet vorliegen 
und müssen zur weiteren Analyse zuvor entfaltet werden (Kap. 2.4.1). Die Faltungen der Ra- 
dialgeschwindigkeiten, die größer als die Nyquistgeschwindigkeit sind, werden am Beispiel 
einer RHI-Messung mit MIRA36-S gezeigt (Abb. 5.1). Am Farbumschlag von maximalen po- 
sitiven Geschwindigkeiten (dunkelrot) auf maximale negative Geschwindigkeiten (dunkelblau) 
ist dieses Artefakt erkennbar. Deutlich sieht man auch, dass die Obergrenze der Wolkenschicht 
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Abbildung 5.1: Gemessene radiale Hydrometeor-Geschwindigkeit Vr,s rnit Faltung vom 
24.1 1 .2009, 10:53 UTC (MIRA36-S, RHI, Azimut a = 288°). 
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Abbildung 5.2: Links: Profil der Inorizontalen Hintergrundgesclnwindigl<eit Ub- IVlitte: Profil der 
mittleren Vertikalgeschwindigkeit wb aus der zenitalen RHI-Messung. Rechts: Aus den Profilen 
(links) berechnete Radialgeschwindigkeiten in Scan-Richtung {a = 288°) 

bei 6.5 km liegt. Bei genauer Betrachtung erkennt man in einer Höhe von gut 2 km einen Ge- 
schwindigkeitssprung. Dies ist die Schmelzschicht (siehe Kap. 4). 

Es ist möglich, die gemessenen Radialgeschwindigkeiten zu entfalten, indem ein radiales Hin- 
tergrundwindfeld abgeschätzt wird (Kap. 2.4.1). Dieses radiale Hintergrundwindfeld wird auf 
Basis von Messungen des Horizontalwindes weiterer Messinstrumente (Abschätzung eines Ho- 
rizontalwindprofils) und ergänzend mit Messungen der Vertikalgeschwindigkeit der Streuer mit 
dem Wolkenradar (Abschätzung eines Vertikalwindprofils) bestimmt. Nach der Entfaltung nach 
Hennington (1981) (Kap. 2.4.1) werden im letzten Schritt des Entfaltungsvorgangs fehlerbehaf- 
te Werte, die im Übergangsbereich der Faltung liegen, für die weitere Bearbeitung ausgeschlos- 
sen. Schrittweise wird für die Entfaltung folgendermaßen vorgegangen: 

1. Berechnung des Horizontalwind-Profils Ub{z): Mit einem in der Nähe stehenden C-Band 
Radar werden im operationeilen Betrieb jede zehn Minuten Volumenscans durchgeführt 
und mittels des VVP' -Algorithmus von Waldteufel und Corbin (1979) die Richtung 
und Stärke des Horizontalwinds bestimmt. Diese Werte sind repräsentativ für ein Volu- 
men, das aufgespannt wird mit Elevationen zwischen 0.4° und 30° und einer maximalen 
Reichweite von 120 km. Nun kann der Fall eintreten, dass das C-Band Radar in einem 
Höhenintervall keine genügend großen Streuer detektiert, das Wolkenradar aber doch. 
Dann werden Windprofile aus operationeilen Radiosondenmessungen verwendet. Dabei 
sollten möglichst viele Stationen verwendet werden, die nahe am Radarstandort liegen. 
Der Grund dafür ist, dass standortbezogene Abweichungen der Windrichtung und Wind- 
geschwindigkeit durch sich schnell ändernde meteorologische Bedingungen, beispiels- 
weise bei Frontdurchgängen, minimiert werden sollen. Im vorliegenden Fall wurden Sta- 

' Velocity Volume Processing 
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Abbildung 5.3: Wie Abb. 5.1, links: Nyquistzahl n mit Werten von -1 (rot), 0 (grün), 1 (blau), 
rechts: entfaltete Radialgeschwindigkeiten 

tionen gewählt, die möglichst nahe am Radarstandort liegen und an denen mindestens je- 
de 12 Stunden Radiosonden gestartet werden. Zum einen ist dies Stuttgart/Schnarrenberg, 
die ca. 63 km in süd-östlicher Richtung vom Radarstandort entfernt liegt, zum anderen 
die Station Idar-Oberstein, die ca. 106 km in nord-westlicher Richtung vom Radarstand- 
ort entfernt ist. Die Protile beider Stationen werden für jeden Aufstieg (Startzeiten sind 
OUTC und 12UTC) gemittelt und zeitlich interpoliert auf den Zeitpunkt des jeweiligen 
RHI-Scans. Um aus den mit Radar und Radiosonde gemessenen Geschwindigkeiten ein 
einheitliches Windprofil zu erstellen, werden die Profile zunächst auf ein Gitter interpo- 
liert, das den Höhenstufen des Wolkenradars entspricht. Die Abschätzung des Horizon- 
talwinds in Richtung des RHI-Scans Ub{z) ist dann gegeben durch 

Ub{z) = uc,r{z) cos[ac,R{z) - a] (5.1) 

Dabei ist a der Azimut des RHI-Scans, ucfi{z) die kombinierte horizontale Windge- 
schwindigkeit, ac,R{z) die kombinierte horizontale Windrichtung (Indizes C, R und b 
stehen jeweils für C-Band, Radiosonde und background). Ein Profil-Beispiel für U}){z) 
ist für den 24.1 1.2010, 10:53 UTC in Abb. 5.2, links, gezeigt. Die Komponente des Hori- 
zontalwindes in Richtung des RHI-Scans (a = 288 °) nimmt im untersten Kilometer stark 
mit der Höhe zu und überschreitet 20 m/s, stagniert bis zu einer Höhe von ca. 4.5 km und 
steigt dann mit der Höhe weiter an. 

2. Berechnung des Profils der Vertikalgeschwindigkeit Wh{z): Da keine Messungen der Ver- 
tikalgeschwindigkeit existieren, wird das Profil aus Messungen des Wolkenradars ab- 
geleitet. Hierfür werden die radialen Geschwindigkeiten der RHI-Messungen des Wol- 
kenradars für Zenitwinkel von ez = 0 ° zt 4 ° gemittelt. Der Zenitwinkelbereich wur- 
de so gewählt, dass einerseits möglichst viele Messungen berücksichtigt werden können 
und andererseits der Betrag der Horizontalgeschwindigkeit des Windes zur radialen Ge- 
schwindigkeit gering bleibt (bei einer Auflösung von Aez = 1.23° (Kap. 3.2) stehen 
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Abbildung 5.4: Links: Doppler-Streubreite, rechts: Faltungsübergangsbereinigte radiale Ge- 
schwindigkeiten 

6 bis 7 Vertikal-Messungen eines RHI-Scans zur Verfügung). Die Variabilität der Werte 
ist zwar noch immer hoch (Abb. 5.2, mitte), dennoch sind klar Bereiche unterschiedli- 
cher Vertikalgeschwindigkeit sichtbar. Ein Bereich langsam fallender Partikel befindet 
sich oberhalb etwa 2.2 km, hier fallen Eispartikel mit ca. 1 m/s. Ein Bereich schnell fal- 
lender Partikel liegt unterhalb 2.0 km, hier fallen Tropfen mit 2 — 4 m/s. Dazwischen 
befindet sich die Schmelzschicht, die mit einer starken Zunahme der Fallgeschwindigkeit 
einhergeht (siehe auch Kap. 4). 

3. Berechnung des radialen Hintergrundwindfeldes vi,{R, ez): Die abgeschätzte mittlere Ho- 
rizontal- und Vertikalgeschwindigkeit des Windes müssen nun noch auf eine Ebene ra- 
dialer Geschwindigkeiten, vergleichbar mit einer RHI-Messung, projiziert und addiert 
werden zu 

Wb(i?, €z) = Ub{z) sm{ez) + Wb{z) cos(ez), (5.2) 

mit z = R cos(e2^). Ein Beispiel für das berechnete radiale Hintergrundwindfeld auf 
Basis der Hintergrundprofile der mittleren Horizontal- und Vertikalgeschwindigkeit ist, 
konsistent zu Abb. 5.1, in Abb. 5.2, rechts, für den 24.1 1.2009, 10:53 UTC gezeigt. 

4. Mit Hilfe des radialen Hintergrundwindfeldes folgt die Berechnung der Nyquistzahl (Gl. 
2.1 1). Damit sind die Faltungsbereiche und die Richtung der Faltung detektiert (Abb. 5.3, 
links) und die gemessene Radialgeschwindigkeit kann entfaltet werden (Gl. 2.12). Das 
Ergebnis der Entfaltung ist in Abb. 5.3, rechts dargestellt. Einzig am Übergangsbereich 
der Faltung sind Artefakte zu erkennen, die sich von den umgebenden Werten deutlich 
abheben. 

5. Bereinigung des Faltungsübergangs: Der Faltungsübergang zeichnet sich dadurch aus, 
dass die Geschwindigkeiten der Streuer in jedem Range Gate oberhalb und unterhalb der 
Nyquistgeschwindigkeit (Kap. 2.4.1) liegen und so die spektrale Breite deutlich erhöhen 
(Abb. 5.4, links, brauner vertikaler, linienartiger Bereich). Mit der Wahl eines empirischen 
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Grenzwertes für die spektrale Breite von 2.2 m/s werden nahezu alle artifiziellen Radi- 
algeschwindigkeiten am Faltungsübergang detektiert und zur weiteren Verarbeitung der 
Geschwindigkeitsdaten nicht mehr berücksichtigt (Abb. 5.4, rechts). Eine Interpolation 
der dadurch entstandenen Lücken wurde bewusst vermieden. Nach der Berechnung der 
residualen Geschwindigkeiten (nachfolgendes Kapitel) und der anschließenden Betrach- 
tung einzelner Strukturen würden gemittelte Bereiche die Deutung der Werte zusätzlich 
erschweren. 



5.1 Berechnung der residualen Geschwindigkeiten 

Allein mit den entfalteten radialen Geschwindigkeiten lässt sich noch nicht sehr viel über die 
Eigenschaften und die Entwicklung der gemessenen Wolke, bzw. deren Niederschlag aussagen. 
Direkt ableiten lässt sich das Profil der Vertikalgeschwindigkeit im Zenit (Abb. 5.2, mitte) und, 
mit der Annahme, dass der horizontale Anteil der Radialgeschwindigkeit der bestimmende ist, 
bei einem Zenit- Winkel von 45 Grad, die mittlere Horizontalgeschwindigkeit der Streuer. Bei- 
spielsweise anhand Abb. 5.4, rechts, läßt sich die Horizontalgeschwindigkeit in einer Höhe von 
6 km zu ca. (18 m/s)/ sin(45 °) = 25.5 m/s abschätzen. 

Die Vertikalgeschwindigkeit im Zenit ist weder die mittlere Vertikalgeschwindigkeit im gesam- 
ten Messbereich, noch charakterisiert sie vertikale Umlagerungen im gesamten Messbereich, 
sondern lediglich die vertikalen Bewegungen in einer, im Vergleich zum gesamten Messbe- 
reich, sehr kleinen Fläche (bzw. eine Linie). Wie bereits im vorigen Kapitel beschrieben, lässt 
sich mit dieser Vertikalgeschwindigkeit zumindest der Niederschlagsbereich flüssiger Hydro- 
meteore (Regen), der mit höheren Vertikalgeschwindigkeiten einhergeht (Abb. 5.2, mitte), von 
dem Bereich trennen, in dem vorwiegend Eispartikel vorliegen. Ob vertikale Umlagerungen in- 
nerhalb der Wolke stattfinden und in welchem Umfang, lässt sich mit den entfalteten radialen 
Daten nicht direkt bestimmen. Deshalb wird nun eine neue Größe eingeführt, mit deren Hilfe 
es möglich sein wird, diese Umlagerungen und deren Entwicklung abzuschätzen. Diese Größe 
wird nachfolgend als residuale Geschwindigkeit bezeichnet. 

Ausgangspunkt zur Berechnung der residualen Geschwindigkeit ist die Beobachtung, dass 
die entfalteten radialen Geschwindigkeiten im Wesentlichen von großräumigen, horizontalen 
Windfeldem und dem mittleren Profil der Vertikalgeschwindigkeit der Streuer abhängen. Diese 
großräumigen Strukturen sollen von den radialen Doppler-Geschwindigkeiten separiert wer- 
den, um Strukturen auf kleineren Skalen, gekennzeichnet durch die residuale Geschwindigkeit, 
sichtbar zu machen. 

Zur Separation der mittleren horizontalen und vertikalen Komponenten der Radialgeschwin- 
digkeiten ist es sinnvoll, die Daten auf ein gleichmäßiges kartesisches Gitter zu interpolieren. 
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Für Daten des Wolkenradars ist die Gitterlänge in x- und ^-Richtung jeweils 30 m und ent- 
spricht somit der Länge eines Range Gates. Damit wird eine schichtweise Analyse der Daten 
gewährleistet. D.h. für jedes Höhenintervall werden die Windfelder separat abgeleitet mit nach- 
folgend beschriebener Methode. Hierzu wird der Ansatz gemacht, dass sich die in kartesischen 
Koordinaten vorliegenden entfalteten radialen Geschwindigkeiten Vr{x, z) formal aufgespalten 
lassen in ein großskaliges Windfeld v*{x, z) und ein residuales Windfeld v^iesi^, z) gemäß 



Vr{x, Z) = V*{X, Z) + ÜResCiC, z) (5.3) 

Das großskalige Windfeld wird konzeptionell zerlegt in einen Anteil, der von der lokalen Ho- 
rizontalgcschwindigkeit «(.x, 2) abhängt und in einen Anteil, der von der lokalen Vertikalge- 
schwindigkeit w{x, z) abhängt. Für die lokalen Geschwindigkeiten wird eine lineare Näherung 
angenommen: 



du 

Ü{X, Z) = Um{z) + -^{z) X 

? (5.4) 

W{X, Z) = Wm{z) + ^{Z) X 
OX 



Dabei sind Ujn{z) und Wjn{z) jeweils die mittleren Horizontal- und Vertikalgeschwindigkeiten, 
du/dx{z) und dw/dx{z) sind die jeweiligen mittleren Gradienten der Horizontal- und Verti- 
kalgeschwindigkeiten in x-Richtung. Einzelne Terme der rechten Seite (Gl. 5.4), bzw. Kom- 
binationen davon werden aus den entfalteten Radialgeschwindigkeiten durch die Minimierung 
einer Kostenfunktion, wie im Folgenden ausführlich beschrieben, abgeleitet. 

Terme, die Ableitungen nach z enthalten, werden nicht berücksichtigt, da während der Berech- 
nung schichtweise vorgegangen wird und jedes Höhenintervall separat berücksichtigt wird. Ein 
separiertes Höhenintervall enthält keine Information über die vertikale Änderung der lokalen 
Geschwindigkeiten. Weitere Terme der lokalen Geschwindigkeit werden vernachlässigt, mit 
der Annahme, dass diese klein sind, verglichen mit den gewählten Termen. Das großskalige 
Windfeld ist durch die Approximation definiert zu 



v*{x, z) = ü{x, z) sm{ez{x, z)) + w{x, z) cos{ez{x, z)) 

du \ ( dw 

Um{z) + —{z)x\ siB.{ez{x,z)) + iwm{z) + —{z)x ) cos{ez{x,z)) 



(5.5) 
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mit dem Zenitwinkel ez{x, z), der horizontalen Koordinate x und der vertikalen Koordinate z. 
Mit X = z tan(e^(a;, z)) kann Gleichung 5.5 geschrieben werden als 



*f X f \ ■ f ^ ^^ , ^ s^^'^i^z{x, z)) 
V [x, z) = Um{z) sm[ez{x, z)) + —{z) z 



dx cos(e^(a;, z)) 

dw 

+w^{z) cos{ez)+ -Q^i^)^ su\{ez{x,z)) 



Um{z) + ^i^) ^ sin(e^(a:, z)) 
+ I Wm{z) - ^(z) cos{ez{x, z)) 



dx ) Q.o^{ez{x,z)) 

Aus Messungen können dann mit der Methode der kleinsten Quadrate drei linear unabhängige 
Terme 



«m(^) + l^(^) ^ J • \Wyy^{z) - ^{z) z) , (^{z) Z ) . (5.7) 

ox J \ ox J \ox 



berechnet werden. Dies wird realisiert mit der Minimierung der Kostenfunktion 



Ji^) = Y^Mxi, z) - v*{xi, z)y (5.8) 

i=l 

Dabei entspricht i dem jeweiligen Gitterpunkt auf der Abszisse. 

Es stehen nicht immer in jedem Höhenintervall ausreichend Daten zur Verfügung. Dadurch 
können die Abweichungen der berechneten drei Terme von einem Höhenintervall zum nächsten 
groß werden. Um dem entgegen zu wirken, wird vor der Berechnung der drei Terme für jedes 
Höhenintervall ein gleitendes Mittel für Vr{x, z) über jeweils zwei benachbarte Höhenintervalle 
darüber und darunter berechnet. Das heißt, die Werte der radialen Geschwindigkeit Vr(x, z) 
entsprechen dem Mittel der Werte zwischen Zk-2 bis Zk+2- Wenn dennoch für diese Vr{x, z) 
weniger als ein Drittel der möglichen Datenmenge pro Höhenintervall vorhanden ist, wird keine 
Auswertung durchgeführt. 

Durch Addition des zweiten zum dritten Term (Gl. 5.7) kann die mittlere Vertikalgeschwin- 
digkeit der Streuer Wyn{z) für jedes Höhenintervall direkt abgeschätzt werden. Mit der Annah- 
me, dass die mittlere Horizontalgeschwindigkeit der Streuer deutlich größer ist, als die mittlere 
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Änderung der Vertikalgeschwindigkeit, d.h. Umiz) » z ^(z), kann diese ebenfalls für jede 
Höhenstufe abgeschätzt werden. 

Sind die drei Terme (Gl. 5.7) bekannt, kann unter Verwendung von Gl. 5.6 folglich das großska- 
lige Windfeld v*{x, z) berechnet werden. Wird dann das großskalige Windfeld vom gemessenen 
Windfeld subtrahiert, erhält man die residualen Geschwindigkeiten (Gl. 5.3). 

Mit dieser Methode können für jeden RHI-Scan drei verschiedene Geschwindigkeitsgrößen 
abgeschätzt werden: 

• Ein Profil der mittleren horizontalen Streuergeschwindigkeit in Scan-Richtung, die 
näherungsweise dem Horizontalwind in Scan-Richtung entspricht (Abb. 5.5, oben links). 

• Ein Profil der mittleren Vertikalgeschwindigkeit der Streuer, die sich aus der Sedimen- 
tationsgeschwindigkeit der Streuer und aus der vertikalen Luftbewegung zusammensetzt 
(Abb. 5.5, oben rechts). 

• Ein zweidimensionales Feld der residualen Geschwindigkeit. Diese Geschwindigkeit re- 
sultiert aus der radialen Geschwindigkeit abzüglich der Anteile der mittleren Horizon- 
talgeschwindigkeit, der mittleren Vertikalgeschwindigkeit und des mittleren horizontalen 
Trends der horizontalen und vertikalen Geschwindigkeiten (Abb. 5.5, unten). 

Diese Größen werden im Folgenden zur Analyse von gemessenen Strukturen in den gemessenen 
Reflektivitäten und zur Deutung dynamischer Vorgänge in Wolken genutzt. 

Ein Beispiel der aus Vj.{x, z) abgeleiteten Größen u^{z), w^{z), und v-^si^, z) ist in Abbildung 
5.5 für eine Messung vom 24.11.2009, 10:53 UTC dargestellt. Die Horizontalgeschwindigkeit 
Ujn{z) zeigt dabei von unten nach oben eine starke Zunahme der Geschwindigkeit mit der Höhe 
bis zu einer Höhe von 1.5 km, bis 5.5 km sind die Änderungen der Horizontalgeschwindig- 
keit gering und steigen oberhalb, bis zur Oberkante der Wolke in 6.4 km wieder an. Qualitativ 
stimmt dies mit dem Profil des Hintergrundwinds überein (Abb. 5.2, links). Die mittlere Ver- 
tikalgeschwindigkeit liegt unterhalb 2 km um —3.5 m/s und oberhalb 2.2 km um —1m/s. Da- 
zwischen, im Bereich der Schmelzschicht, zeigt sich ein starker Gradient der Vertikalgeschwin- 
digkeit mit der Höhe. 

Die Werte der residualen Geschwindigkeit (Abb. 5.5, unten) zeigen ein deutlich differenzier- 
teres Bild als die Werte der entfalteten radialen Geschwindigkeit (Abb. 5.4) und einzelne 
teilweise abgeschlossene Strukturen sind zu erkennen. Auf Gebiete, in denen sich die Streu- 
er relativ zu v*{x,z) um bis zu 0.6m/s vom Radar weg bewegen (gelb-rot gekennzeichnet, 
z.B. in 4km < z < 5 km, 0 km < x < 2 km) folgen Gebiete, in denen sich die Streu- 
er mit bis zu —0.6m/s relativ auf das Radar zu bewegen (grün-blau gekennzeichnet, z.B. in 
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Abbildung 5.5: Oben: Mittlere horizontale Windgeschwindigkeit Um{z) in RHI-Scanrichtung, 
daneben: Mittlere Vertikalgeschwindigkeit w,n{z)- Unten: Residuale Geschwindigkeit VRes{x,z) 
in m/s, blau-grün bedeuten Bewegungen auf das Radar zu und gelb-rot Bewegungen vom 
Radar weg, grau sind Fehldaten (zu geringe Datenmenge, Faltungsübergang). 
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3 km < z < A km, 0 km < x < 2 km). Die größte Variation der Werte ist unterhalb der 
Schmelzschicht, im Regen zu erkennen. Grau gekennzeichnet sind Gebiete, in denen Werte 
der Reflektivität vorliegen, aber keine Werte der residualen Geschwindigkeiten. In diesen Ge- 
bieten ist entweder die Datengrundlage zu schwach, um die residualen Geschwindigkeiten zu 
berechnen (oberhalb 6.5 km), oder die Gebiete liegen im Faltungsübergang. Die Interpretation 
der Werte der residualen Geschwindigkeit für verschiedene meteorologische Zustände soll in 
den nächsten Kapiteln diskutiert werden. 

5.2 Identifikation von Strulituren mittels residualer Geschwindigkeit 

Wie im vorigen Abschnitt gezeigt, ist die residuale Geschwindigkeit eine neue Größe, die zur 
Visualisierung dynamischer Strukturen dient und die Interpretation von Strukturen in Wolken 
und deren Entwicklung unterstützen kann. 

Die Auswertung von Daten der residualen Geschwindigkeit werden in den folgenden Unterka- 
piteln anhand einzelner, ausgewählter Zeiträume gezeigt. Lokale Strukturen und wellenartige 
Strukturen in den residualen Geschwindigkeiten, die mit Strukturen in den Reflektivitäten ein- 
hergehen, und dynamische Strukturen, die zwar in den residualen Geschwindigkeiten, nicht 
aber in den Reflektivitäten augenscheinlich sind, werden diskutiert. Zu beachten ist dabei, dass 
es sich um zweidimensionale Messungen handelt, die Windrichtung allerdings nie im komplet- 
ten Messabschnitt der Orientierung des RHI-Scans entspricht. So kann es vorkommen, dass 
starke Änderungen in den gemessenen Größen ihre Ursache in der Advektion von Luftmassen 
in einem Winkel abweichend zum RHI- Azimut haben. 

5.2.1 Lokale Strukturen 

Lokale, kleinräumige Strukturen haben ihre Ursache oft in turbulenten Aufwinden und wur- 
den in der Vergangenheit schon mehrfach untersucht. In diesem Zusammenhang wurden be- 
reits räumlich eindimensionale Messungen bodengestützter, vertikal ausgerichteter Doppier- 
Radare durchgeführt, z.B. Frisch et al. (1995), mit einem 35 GHz-Radar in maritimen strati- 
formen Wolken, KoUias und Albrecht (2000) mit einem 95 GHz-Radar in kontinentalen Stra- 
tocumulus, sowie KoUias et al. (2001) in Schönwettercumuli. Außerdem kamen auch zwei- 
dimensionale Dopplermessungen mit zwei 95 GHz-Radaren zum Einsatz, die vom Flugzeug 
aus vorgenommen wurden (z.B. Damiani et al. (2006) in Cumulus-Wolken sowie Geerts et al. 
(2006) in einer Kaltfront). Es konnten bereits zweidimensionale Geschwindigkeitsfelder aus 
flugzeuggestützten Dual-Doppler-Messungen abgeleitet werden. Dabei war der Strahl eines Ra- 
dars senkrecht nach unten und der Strahl eines zweiten Radars in Flugrichtung schräg nach vor- 
ne ausgerichtet. Ein Nachweis von Wirbeln in Cumulus-Wolken konnte damit erbracht werden 
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(Damiani et al., 2006). Die räumliche und zeitliche Entwicklung einzelner Strukturen innerhalb 
einer Wolke konnte allerdings nicht zufriedenstellend verfolgt werden, da die Flugzeugmes- 
sungen nur sporadisch erfolgten. Demgegenüber ist diese Möglichkeit mit einem Wolkenradar 
mit schwenkbarer Antenne geschaffen. Mit Hilfe der zeitlich hoch aufgelösten RHI-Messungen 
können dynamische Prozesse, bzw. die Veränderung von Strukturen in Wolken, sichtbar durch 
die residualen Geschwindigkeiten, zeitlich verfolgt werden. 

Eine einzelne Struktur, die in den Reflektivitäten sichtbar wird, sowie die Erhöhung der Reflek- 
tivität einer vorhandenen Struktur kann ihre Ursache in vielfältiger Weise haben, beispielsweise 
in der horizontalen und/ oder vertikalen Advektion von Hydrometeoren großer Durchmesser 
oder von Hydrometeoren sehr großer Anzahl. Auch die horizontale Advektion feuchter Luft 
oder das Anheben einer Luftmasse kann durch Kondensations-, bzw. Resublimationswachstum 
die Reflektivität erhöhen. Turbulente Durchmischung, beispielsweise durch Konvektion, kann 
die Rate der Stoßprozesse der Hydrometeore erhöhen und so zum Wachstum der Partikel bei- 
tragen und die Reflektivität erhöhen. Genauso finden sich auch Prozesse, die zu einer Senkung 
der Reflektivität innerhalb einer Struktur bis zur Auflösung der Struktur führen können. Die Ur- 
sache und Entwicklung einer lokalen Struktur soll am Beispiel der Messungen vom 06.10.2009 
besprochen werden. Der Durchzug eines Bereichs hoher Fluktuationen der radialen Geschwin- 
digkeit, die man anhand der residualen Geschwindigkeiten sehen kann, steht in Zusammenhang 
mit einem Gebiet deutlich abgegrenzter, erhöhter Reflektivitäten und eine damit verbundene 
Verstärkung des Niederschlags unterhalb der Schmelzschicht. Die Eigenschaften dieser Struk- 
tur und deren Änderungen werden im Folgenden diskutiert. 

Die meteorologische Situation ist so, dass das Messgebiet des Wolkenradars vorderseitig eines 
Troges liegt und eine mit dem Trog in Verbindung stehende Wellenstörung den Radarstandort 
überquert. Es werden feuchtwarme Luftmassen aus Süd- West advehiert, die ein Niederschlags- 
gebiet über den Radarstandort führen. Während dieses Ereignisses wurde eine Serie von RHI- 
Scans in süd-westlicher, bzw. nord-östlicher Richtung bei einem Azimut von 255 °, bzw. 75° 
gefahren. Der Zeitraum, in dem die Struktur in den Messungen zu sehen ist, liegt zwischen 
14:07:46 UTC und 14:11:17 UTC. 

Einen ersten Überblick über die meteorologische Situation zeigen Messungen vom gesamten 
Messbereich zu den beiden Zeitpunkten 14:08:16 UTC und 14:09:16 UTC, die also nur 60 s 
auseinanderliegen. Folgende Größen sind hierfür in den nachfolgenden Abbildungen darge- 
stellt: 

1. In den Abbildungen 5.6-5.9 die Reflektivität, die residuale Geschwindigkeit (Kap. 5.1), 
das lineare Depolarisationsverhältnis (LDR) und die spektrale Breite. 

2. In den Abbildung 5.10,5.11 das Profil der Windgeschwindigkeit und Windrichtung, be- 
rechnet mit Hilfe des VAD-Algorithmus (Kap. 2.4.2) aus zwei PPI-Messungen vor und 
nach der RHI-Scan-Serie. 
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3. In Abbildung 5.12 die abgeleiteten Profile der mittleren Horizontalgeschwindigkeit und 
der mittleren Vertikalgeschwindigkeit (Kap. 5. 1 ) in Richtung der gemessenen RHI-Scans. 

Grau waagerecht liniert ist in den Abb. 5.6-5.9 die abgeleitete mittlere Begrenzung der Schmelz- 
schicht (Kap. 4.2) der jeweiligen Messung. Die jeweilige schwarze Linie zum Zeitpunkt 
14:08:16 UTC und 14:09:47 UTC steUt den Verlauf der im Folgenden diskutierten Struktur 
der Werte erhöhter Reflektivität dar und ist zur besseren Orientierung in allen Bildern einge- 
zeichnet. 

Zunächst werden die Daten besprochen, die in den Abb. 5.6-5.9 präsentiert sind. Im gesamten 
Scan-Bereich werden bis zu einer Höhe von 9 km Hydrometeore delektiert. Die Schmelzschicht 
befindet sich zwischen 2970 m und 3270 m und ist gekennzeichnet durch erhöhte Reflekti- 
vitätswerte (bis 10 dBz), deutlich erhöhte LDR-Werte (bis —7 dB) und einem deutlichen Gra- 
dienten der mittleren Vertikalgeschwindigkeit (Abb. 5.12). Auf die mittlere Vertikalgeschwin- 
digkeit wird im weiteren Verlauf noch gesondert eingegangen. 

Oberhalb der Schmelzschicht sind bis zu einer Höhe von 5.5 km Strukturen erhöhter Reflek- 
tivität zu erkennen. Die größte Änderung innerhalb einer Minute erfährt die mit schwarzen 
Linien hervorgehobene Struktur in Form einer Gabel. Während zum Zeitpunkt 14:08:16 UTC 
die linke „Gabelspitze" der Struktur die größeren Reflektivitäten zeigt, ist es zum Zeitpunkt 
14:09:47 UTC deren rechte „Gabelspitze". Im Bereich dieser „GabeP'-Struktur sind in den Wer- 
ten der residualen Geschwindigkeit große Schwankungen von bis zu ±0.8 m/s zu erkennen, die 
sich deutlich von den residualen Geschwindigkeiten im restlichen Messbereich (±0.2 m/s ober- 
halb der Schmelzschicht) abheben. Oberhalb 5.5 km ist die Verteilung aller Werte homogener. 
Weder in den Reflektivitäten noch in den residualen Geschwindigkeiten sind große Variationen 
der Werte erkennbar. 

Das LDR und die spektrale Breite zeigen deutliche Unterschiede hinsichtfich der Bereiche un- 
terhalb und oberhalb der Schmelzschicht. Die niedrigsten Werte des LDR befinden sich un- 
terhalb der Schmelzschicht (ca. -23 dB) und die höchsten Werte innerhalb der Schmelzschicht 
(bis -7 dB). Oberhalb der Schmelzschicht ist eine erhebliche, weitgehend unstrukturierte Va- 
riabilität der LDR-Werte zu sehen. Im Bereich der „GabeF'-Struktur zeichnet sich allerdings 
ein Bereich ab, in dem die LDR-Werte, verglichen mit der Umgebung, meist geringer ausfallen 
und Werte unterhalb —20 dB annehmen. Oberhalb 7 km sind die Werte der Reflektivität im 
Cross-Kanal unterhalb der Detektionsgrenze und es kann kein LDR berechnet werden (grauer 
Bereich). Für die spektrale Breite hingegen sind die größten Werte (bis zu 1 m/s) unterhalb der 
Schmelzschicht zu erkennen. Bereits innerhalb, als auch oberhalb der Schmelzschicht liegen 
die Werte meist unter 0.2 m/s. Höhere schichtartig erscheinende Werte oberhalb der Schmelz- 
schicht, die horizontal vom Zenit zu den Rändern des Messbereichs zunehmen, sind Variationen 
der horizontalen Windgeschwindigkeit. Im Bereich der "Gabel" -Struktur sind lediglich in der 
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Abbildung 5.6: RHI-Scan vom 06.10.2009, 14:08:1 6 UTC (oben) und 60s später, d.h. 
14:09:16 UTC (unten). Dargestellte Größe ist die Reflektivität in [dBz]. Graue Linien: Grenzen 
der Schmelzschicht (Kap. 4). Schwarze Linien: Struktur in der Reflektivität (siehe Text). 



Copyrighied matenal 



68 



5 Dynamik in Wolken 





Abbildung 5.7: wie Abb. 5.6, dargestellte Größe ist die residuale Geschwindigkeit in [m/s]. 
Zusätzliche graue Bereiche: Gebiete in denen Werte der Reflektivität, aber keine Werte der 
dargestellten Größe vorliegen 
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Abbildung 5.8: wie Abb. 5.7, dargestellte Größe ist das lineare Depolarisationsverhältnis in 
[dB]. 
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Horizontale Distanz [Wm] 
Abbildung 5.9: wie Abb. 5.7, dargestellte Größe ist die spektrale Breite in [m/s]. 
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unmittelbaren Nähe der rechten „Gabelspitze" deutlich höhere Werte (bis zu 0.6 m/s) der spek- 
tralen Breite zu erkennen. Da diese erhöhten Werte auch in der Nähe des Zenits auftreten, ist 
die Ursache des breiteren Spektrums in der Variation der Vertikalgeschwindigkeiten zu sehen. 
Wie auch für die residuale Geschwindigkeit ist der Faltungsübergang grau unterlegt. 



Jetzt wenden wir uns der Analyse der Windprofile (Abb. 5.10,5.1 1) zu. Zur Berechnung der 
mittleren Windrichtung, sowie mittleren Windgeschwindigkeit (vgl. Kap. 2.4.2) sind PPI-Mes- 
sungen nötig. Diese wurden 13:31 UTC und 14:34 UTC, jeweils zeitlich vor und nach dem 
Durchzug der „Gabel-Struktur" durchgeführt. Die Windprofile aus PPI-Messungen des Wol- 
kenradars, ermittelt mit Hilfe des VAD-Algorithmus, sind in den Abbildung 5.10,5.1 1 rot dar- 
gestellt. Zum Vergleich wurde das räumliche und zeitlich gewichtete Mittel der Winddaten der 
Radiosondenaufstiege von Stuttgart/Schnarrenberg und von Idar-Oberstein (siehe Anfang des 
Kapitel 5) eingezeichnet (schwarz). Sind zu wenig Daten auf einem Kreisring eines PPI-Scans 
vorhanden, d.h. auf einer Höhenstufe, werden keine Winddaten berechnet. Dies erklärt die Da- 
tenlücken in den Windprofilen unterhalb 3 km. Die Profile der Messungen mittels Wolkenradar 
lassen sich anhand der Windgeschwindigkeit in drei Höhenbereiche einteilen, von unten nach 
oben in einen Bereich zunehmender Werte mit der Höhe, gefolgt von einem Bereich konstanter, 
bzw. leicht fallender Werte und abschließend ein Höhenabschnitt wieder zunehmender Werte. 




Windgeschwindigkeit [m/s] Windriclitung [Grad] 



Abbildung 5.10: Windgeschwindigkeit und Windrichtung vom 06.10.2009, 13:31 UTC. Rot: Be- 
rechnete Werte aus PPI-IVIessungen des Wolkenradars iVIIRASS-S. Schwarz: Zeitliches und 
räumliches Mittel von je zwei Radiosondenaufstiegen. 
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Abbildung 5.11: Wie Abb. 5.10, hier für die Messung vom 06.10.2009, 14:34UTC. 



Von besonderem Interesse ist der zweite Höhenabschnitt, in dem auch die ,,Gabel"-Struktur zu 
finden ist. Dieser Höhenabschnitt erstreckt sich zum Zeitpunkt 13:3 1 UTC von 3 km bis 6.5 km. 
Die Windgeschwindigkeit ist näherungs weise konstant und liegt bei etwa 20 m/s. Eine leichte 
Rechtsdrehung des Windes von 240° auf 250° ist ebenfalls zu erkennen und tritt hauptsächlich 
in einem Höhenintervall zwischen 4 km und 5.5 km auf. Nach ca. einer Stunde, zum Zeitpunkt 
14:34 UTC, erstreckt sich der zweite Höhenabschnitt von etwa 3 km bis 7.5 km. Im Unterschied 
zum vorigen Zeitpunkt nimmt die Horizontalgeschwindigkeit höhere Werte an und sinkt leicht 
mit zunehmender Höhe von ca. 24 m/s knapp oberhalb 3 km auf ca. 21m/s knapp oberhalb 
7.5 km. Auch die Windrichtung hat sich zwischen den Zeitpunkten verändert von einer nahezu 
linearen Rechtsdrehung des Windes hin zu wechselnden Drehrichtungen. Eine Rechtsdrehung 
um ca. 10° findet statt zwischen etwa 3km-4km, gefolgt von einer Linksdrehung um ca. 10° 
zwischen etwa 4km-5.5km. Oberhalb, bis ca. 7.5 km bleibt die Windrichtung annähernd kon- 
stant. Aus der Zunahme der Windgeschwindigkeit lässt sich schließen, dass sich der horizonta- 
le Druckgradient im zweiten Höhenabschnitt erhöht hat. In Luftschichten oberhalb 3km kann 
der beobachtete Wind als geostrophischer Wind angenommen werden. Die Rechtsdrehung des 
als geostrophisch angenommenen Windes mit der Höhe deutet auf Warmluftadvektion und die 
Linksdrehung auf Kaltluftadvektion hin. Somit findet zum Zeitpunkt 13:31 UTC im zweiten 
Höhenabschnitt Warmluftadvektion statt, die sich auch in den dritten Abschnitt fortsetzt. Zum 
Zeitpunkt 14:34 UTC tritt Warmluftadvektion nur noch im unteren Bereich des zweiten Ab- 
schnitts auf, an den sich direkt darüber ein Höhenabschnitt anschließt in dem Kaltluftadvektion 
stattfindet. Dies spricht für eine Labilisierung der Schichtung im zweiten Höhenabschnitt. 
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Abbildung 5.12: Horizontalwind und Vertikalgescliwindigkeit der Streuer in Riclitung 
der RHI-Scans (Azimut: 255 Grad) vom 06.10.2009, 14:08:16UTC (oben) und 60s später, 
14:09:16 UTC (unten). Rot: Berechnete Werte aus RHI-Messungen des Wolkenradars 
MIRA36-S. Schwarz: Zeitliches und räumliches Mittel von je zwei Radiosondenaufstiegen. 
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Die Geschwindigkeiten der Streuer aus RHI-Messungen Um und Wm sind, wie für die bisher 
gezeigten RHI-Messungen (Abb. 5.6-5.9), zum Zeitpunkt 14:08:16 UTC und 14:09:47 UTC in 
Abb. 5.12 dargestellt. Mit horizontalen, grauen Linien eingezeichnet ist die detektierte Schmelz- 
schicht (Kap. 4.2) und schwarz markiert sind die Horizontalgeschwindigkeiten aus Radioson- 
denmessungen in Richtung des jeweiligen RHI-Scans (255 Grad). Negative Werte der horizon- 
talen Geschwindigkeiten bedeuten eine Bewegung der Partikel in den RHI-Scans von links 
nach rechts und positive den umgekehrten Fall. Negative Geschwindigkeiten der vertikalen 
Geschwindigkeiten bedeuten meist sedimentierende Partikel. Die Horizontalgeschwindigkeiten 
aus MIRA36-S schwanken dabei zu beiden Zeitpunkten minimal um die Horizontalgeschwin- 
digkeiten der Radiosondenmessung. 

Auffällig im Profil des Horizontalwindes aus MIRA36-S ist zu beiden Zeitpunkten ein Ma- 
ximum der Windgeschwindigkeit zwischen 2.5 km und 5.3 km. Die Windgeschwindigkeit er- 
reicht dabei Werte von bis zu 24 m/s. Diese Werte entsprechen damit jenen, die 14:34 UTC für 
die zweidimensionale horizontale Windgeschwindigkeit gemessen wurden (Abb. 5.11). Das be- 
deutet, dass die Änderung der Drehrichtung des Windes in einer Höhe von etwa 4 km bereits 
zum Zeitpunkt 14:08 UTC stattgefunden hatte. Damit steht die „Gabel"-Struktur im Zusam- 
menhang mit einer Labilisierung der Luftschichten durch eine höhenabhängige Veränderung 
der Advektion von Luftmassen. 

Die Sedimentationsgeschwindigkeit der Hydrometeore nimmt von oben nach unten von nahe 
0 m/s bis auf 1.3 m/s zu, von der Oberkante der Wolke bis zur Schmelzschichtobergrenze. Zum 
Zeitpunkt 14:09:16 UTC ist zudem eine starke Zunahme der Sedimentationsgeschwindigkeit in 
einer Höhe von etwa 5 km zu erkennen. Diese Zunahme steht vermutlich in Zusammenhang 
mit Änderungen der Zusammensetzung der Hydrometeore im Bereich der „Gabel" -Struktur 
(Abb. 5.6). An der Oberkante der Schmelzschicht erfolgt zu beiden Zeitpunkten eine deutliche 
Zunahme der Sedimentationsgeschwindigkeit und erreicht unterhalb der Schmelzschicht (im 
Regen) Werte von ca. 3.5 m/s. 

Zusammenfassung Die betrachtete Messfläche kann von oben nach unten in vier Schichten 
aufgeteilt werden: Die erste Schicht zwischen der Oberkante der Wolke in ca. 9 km und ei- 
ner Höhe von 5.5 km, die zweite zwischen 5.5 km und der Oberkante der Schmelz schicht in 
3.27 km, die dritte ist die Schmelzschicht selbst und die vierte Schicht befindet sich unterhalb 
der Schmelzschicht, unterhalb 2.97 km. 

Wenden wir uns zunächst den unteren beiden Schichten zu. In der dritten Schicht, der Schmelz- 
schicht, findet der Phasenübergang der fallenden Hydrometeore von Eispartikeln zu Wasser- 
tropfen statt (Kap. 4). Diese Schicht zeichnet sich aus durch die höchsten Werte der Reflekti- 
vität und des LDR und einem starken Gradienten der mittleren Vertikalgeschwindigkeit und der 
spektralen Breite. 
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In der vierten Schicht befinden sich hauptsächlich flüssige, schnell fallende Tropfen, es regnet. 
Dass die Hydrometeore unterhalb der Schmelzschicht kugelförmig sind, zeigt sich an niedrigen 
Werten des LDR. Ebenfalls charakteristisch für Regen sind hohe Werte der spektralen Breite: 
Liegen Hydrometeore als Tropfen vor, ist ihre Sedimentationsgeschwindigkeit hoch und ih- 
re Form nahezu kugelförmig. Mit zunehmender Sedimentationsgeschwindigkeit variiert diese 
stark und die spektrale Breite ist erhöht. Oberhalb der Schmelzzone dagegen sind die Eisparti- 
kel nicht kugelförmig, taumeln leicht und deren Reibungswiderstand an Luft ist deutlich erhöht. 
Dadurch ist die Variabilität der möglichen Sedimentationsgeschwindigkeiten der Hydrometeore 
und somit die spektrale Breite oberhalb der Schmelzschicht deutlich geringer als unterhalb. 

Oberhalb der Schmelzschicht müssen die später fallenden Hydrometeore zuerst gebildet wer- 
den. Dies passiert bereits in der ersten Schicht. In dieser Schicht sind niedrige Werte und nur 
geringe Variationen der Reflektivitäten zu erkennen. Ebenso nehmen die Werte der residualen 
Geschwindigkeit nur geringe Werte an. Daraus geht hervor, das keine vertikalen Austausch- 
prozesse im Messbereich stattfinden und diese Wolkenschicht als stratiform angenommen wer- 
den kann. Die von oben nach unten zunehmenden Werte der Reflektivität, die mit einer zu- 
nehmenden Sedimentationsgeschwindigkeiten einhergehen, lassen darauf schließen, dass die 
Größe der Eispartikel von oben nach unten zunimmt. Weiterhin findet in dieser Schicht eine 
Rechtsdrehung des Horizontalwindes mit der Höhe statt, die zum Zeitpunkt 13:31 UTC ausge- 
prägter ausfällt als zum Zeitpunkt 14:34 UTC. Dies lässt darauf schließen, dass in dieser Schicht 
Warmluftadvektion erfolgt, die sich im Laufe der Zeit abschwächt. 

Die interessanteste Schicht ist die zweite und befindet sich oberhalb der Schmelzschicht. Auch 
hier ist davon auszugehen, dass die meisten Hydrometeore als Eispartikel vorliegen. In dieser 
Schicht sind die höchsten positiven und niedrigsten negativen residualen Geschwindigkeiten zu 
identifizieren. Da diese auch im Zenit sichtbar sind, ist der vertikale Anteil der residualen Ge- 
schwindigkeiten in dieser Schicht der dominierende. Die Bereiche maximaler, bzw. minimaler 
residualer Geschwindigkeiten sind entweder auf Bereiche unterschiedlich großer Streuer oder 
auf vertikale Luftbewegungen, oder auf einer Mischung aus beidem zurückzuführen. So sind 
im Umfeld einer „Gabel" -Struktur Bereiche zu erkennen, in denen hohe Werte der Reflekti- 
vität mit negativen Werten der residualen Geschwindigkeiten einhergehen. Dies spricht dafür, 
dass diese Bereiche im Vergleich zur Umgebung durch größere Eispartikel bestimmt werden. 
Daneben gibt es innerhalb der „Gabel" -Struktur aber auch Bereiche hoher Reflektivität, die 
mit hohen positiven residualen Geschwindigkeiten in Zusammenhang stehen. Hier sind Auf- 
windbereiche wahrscheinlich, die die Sedimentation der Partikel verlangsamen. Finden verti- 
kale Austauschprozesse statt, ist der stratiforme Charakter der Wolke in dieser Schicht nicht 
mehr gegeben. Die Variation der Reflektivitäten ist ebenfalls größer, und kleinräumige zusam- 
menhängende Strukturen hoher Reflektivität zeichnen sich ab. Hervorstechend ist hierbei die 
Struktur in Form einer Gabel, die sich im Laufe der Zeit in ihren Reflektivitätswerten stark 
verändert. Dabei nimmt die Reflektivität im Bereich der „Gaber'-Struktur in Höhen zwischen 



76 



5 Dynamik in Wolken 



4 km und 5 km innerhalb einer Minute um mehrere dB zu. Gleichzeitig nimmt die Sedimenta- 
tionsgeschwindigkeit in dieser Schicht um mehrere 10 cm/s zu. Daraus läßt sich folgern, dass 
die Größe der Hydrometeore im oberen Bereich der „GabeF'-Struktur zunimmt. Die „Gabel"- 
Struktur steht zudem in direktem Zusammenhang mit erhöhten Werten und stärkeren Variatio- 
nen der residualen Geschwindigkeiten und in Teilbereichen mit der spektralen Breite. Am un- 
teren Ende dieser Struktur schließt sich ein Bereich intensiveren Niederschlags an, sichtbar in 
den Schichten drei und vier durch erhöhte Reflektivitäten, negative Residualgeschwindigkeiten 
und einer erhöhten spektralen Breite. Die Werte des LDR sind schwieriger interpretierbar, da 
immer eine Mischung aus Hydrometeoren verschiedener Formen gemessen wird, z.B. oberhalb 
der Schmelzschicht von Graupel, unterkühlte Wassertropfen, Dendrite, Plättchen und weitere 
Eiskristalle und -aggregate. Im Bereich der „Gabel" -Struktur fallen die LDR-Werte im Ver- 
gleich zur Umgebung meist niedriger aus. Die anschauliche Folgerung ist, dass die Eispartikel 
innerhalb der „Gabel"-Struktur tendenziell mehr oder weniger rund sind, oder die Eispartikel 
aufgrund turbulenter Umlagerungen stärker taumeln als außerhalb der Struktur. Ist die Form der 
Partikel im Bereich der „GabeF'-Struktur runder, befinden sich möglicherweise neben Eiskri- 
stallen, -Partikeln auch unterkühlte Wassertropfen im Messvolumen. 

Die schnelle und deutliche Zunahme der Partikelgrößen im oberen Bereich der „GabeF'- 
Struktur hat ihre Ursache zum einen in der seitlichen Advektion großer Partikel in das zwei- 
dimensionale Messgebiet und zum anderen möglicherweise im homogenen Wachstum der Par- 
tikel durch ein größeres Feuchteangebot aufgrund von kleinräumigen Hebungsprozessen vor- 
derseitig der „GabeF'-Struktur Das homogene Wachstum der Eispartikel kann durch vorhan- 
dene unterkühlte Wassertropfen noch beschleunigt werden. Die vordcrscitigc Hebung ist an- 
hand hoher positiver Residualgeschwindigkeiten vorderseitig der „GabeF'-Struktur angedeu- 
tet. Aufgrund der getroffenen Annahmen ist davon auszugehen, dass die ,, GabeF'-Struktur in 
Verbindung mit einem frontähnlichen Bereich konvektiver Einlagerungen in das großräumige 
Hebungsgebiet steht. 

5.2.2 Wellenartige Strukturen 

Neben lokalen Strukturen in Wolken werden zeitweise auch periodisch angeordnete, wellen- 
artige Strukturen beobachtet, die mit Hilfe der residualen Geschwindigkeiten untersucht wer- 
den können. Wellenartige Strukturen sind gekennzeichnet durch eine Abfolge starker Reflekti- 
vitäten gekoppelt mit Auf- und Abwindbereichen. Am 08.01.2010 wurden solche wellenartige 
Strukturen sowohl in den Reflektivitäten als auch in den residualen Geschwindigkeiten beob- 
achtet. Die großräumige Wetterlage, die diese Strukturen im Messbereich verursachte, wurde 
an diesem Tag bestimmt durch ein "Abtropfen" eines Sekundärtrogs von einem Haupttrog über 
Skandinavien. Dieser Vorgang löste vorderseitig des Trogsystems eine Zyklogenese aus. Das 
dadurch entstandene, sich verstärkende Tiefdruckgebiet advehierte wärmere Luftmassen in der 
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Abbildung 5.13: Temperaturprofil in Grad Celsius der Radiosondenaufstiege in Stutt- 
gart/Schnarrenberg (ST, schwarz) und Idar-Oberstein (10, grün) vom 08.01.2010, 12:00 UTG 

mittleren Troposphäre aus Süd-West. Die wärmere Luft glitt dabei über kältere Luft polaren 
Ursprungs, die in der unteren Troposphäre um ein Hoch über Skandinavien geführt und zum 
Radarstandort aus nordöstlicher Richtung advehiert wurde. Der dadurch initiierte, großräumige 
Hebungsprozess führte zu leichtem Schneefall. Die Temperatur lag ganztägig in allen Höhen 
der Troposphäre bei Werten deutlich unter 0°C (siehe Abb. 5.13). 

Anhand der Temperaturprofile (Abb. 5.13) ist zu erkennen, dass sich am Standort Idar-Oberstein 
(lO) die kältere Luftmasse bereits bis in Höhen von 4.5 km bemerkbar macht. Die Temperatu- 
ren liegen um bis zu 6°C unter denen am Standort Stuttgart/Schnarrenberg (ST). Oberhalb einer 
Höhe von etwa 4.5 km entspricht der Verlauf des Temperaturprofils dem einer Feuchtadiabate 
und der vertikale Temperaturgradient (Abb. 5.13) ist ca. 0.7 K/ 100 m. Übereinstimmend zeigen 
beide Profile eine Inversion knapp unterhalb einem Kilometer. Direkt darüber ist die Tempera- 
turabnahme näherungsweise feuchtadiabatisch, bevor die nächste Inversionsschicht folgt. Die- 
se befindet sich am Standort Stuttgart/Schnarrenberg in einer Schicht knapp unterhalb 2.5 km. 
Darüber ist eine nahezu isotherme Schicht zu erkennen und die Temperaturabnahme mit der 
Höhe liegt lediglich bei 0.2 K/100 m. Oberhalb von 3.5 km findet der Übergang zum weite- 
ren feuchtadiabatischen Temperaturverlauf statt. Im Profil der Temperatur am Standort Idar- 
Oberstein hingegen befindet sich die Inversionsschicht knapp unterhalb 3 km und ist schwächer 
ausgeprägt. Weitere Inversionsschichten folgen in ca. 3.3 km und ca. 4.5 km. 

Die verschiedenen Richtungen der großräumigen Advektion von Luftmassen unterschiedlicher 
Herkunft sind in Abbildung 5.14 zu erkennen. Dargestellt ist die horizontale Windgeschwin- 
digkeit (links) und die Windrichtung (rechts). Die roten Profile basieren auf einem PPI-Scan 
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Abbildung 5.14: Windgeschwindigkeit und Windrichtung vom 08.01.2010, 14:1 8 UTC. Rot: Be- 
rechnete Werte aus PPI-IVIessungen des Wolkenradars IVIIRASS-S. Schwarz: Zeitliches und 
räumliches Mittel von je zwei Radiosondenaufstiegen (Vgl. Abb. 6.3: Einzelne Sondenaufstie- 
ge zum Zeitpunkt 08.01.2010, 12:00 UTC). 

mit dem Wolkenradar MIRA36-S vom 08.01.2010, 14:18 UTC, bei einem Zenitwinkel von 
ez = 45°. Die schwarzen Profile basieren auf Messungen der beiden Radiosondenaufstiege 
in Stuttgart/Schnarrenberg und Idar-Oberstein. Dabei wurde ein zeitliches und räumliches Mit- 
tel von jeweils zwei Radiosondenaufstiegen für jeden Standort berechnet und gleichermaßen 
vorgegangen, wie auch schon zuvor für die Berechnung des Hintergrundwindfelds. Die Advek- 
tion von Luftmassen aus südwestlicher Richtung ist oberhalb einer Höhe von 3.5 km und einer 
Windrichtung (Abb. 5.14, rechts) von etwa 210 ° zu sehen. Die Advektion von Luftmassen aus 
nord-östlicher Richtung ist unterhalb 2 km und einer Windrichtung von 30 ° — 60 ° zu erkennen. 
Zwischen 2 km und 3.5 km findet eine Rechtsdrehung des Windes um 150 ° statt. 



Bei Betrachtung der horizontalen Windgeschwindigkeit aus Messungen des Wolkenradars 
(Abb. 5.14, links) fällt auf, dass unterhalb der Scherungsschicht maximale Geschwindigkei- 
ten von 17 m/s-2ü m/s in einer Schicht zwischen 1 km und 2 km erreicht werden und minimale 
Geschwindigkeiten innerhalb der Scherungsschicht in einer Höhe von etwa 3 km und Hori- 
zontalgeschwindigkeiten zwischen 2 m/s und 3 m/s . Der Höhenabschnitt, in dem die stärksten 
Änderungen der Windgeschwindigkeit und der Windrichtung auftreten erstreckt sich von 
2.6km bis 3.2 km. Dieser Höhenabschnitt wird im weiteren Verlauf als Scherungsschicht be- 
zeichnet. 
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Zeitlich vor der PPI-Messung wurden RHI-Scans in nordöstlicher, bzw. südwestlicher Richtung 
durchgeführt. Der Azimut betrug 58 °, bzw. 238 °. Ein RHI-Scan zum Zeitpunkt 13:44:31 UTC 
ist in den Abb. 5.15, 5.16 dargestellt. Die gezeigten Größen sind die Reflektivität, die residuale 
Geschwindigkeit, das lineare Depolarisationsverhältnis und die spektrale Breite. Die Obergren- 
ze der Wolke wird in einer Höhe von etwa 5.5 km detektiert. Der Messbereich kann in drei ver- 
schiedene Schichten eingeteilt werden: a) Eine Schicht geringer Gradienten der Reflektivitäten 
in einer stratiformen Region zwischen 3.1 km und der Obergrenze der Wolke, b) die Scherungs- 
schicht mit periodischen Fluktuationen der Reflektivität zwischen 2.6kmund 3.1 km und c) eine 
Schicht unterhalb 2 km, mit einem minimalen horizontalen und einem deutlich sichtbaren verti- 
kalen Gradienten der Reflektivität. Die mittlere Schicht (b) der periodischen Fluktuationen der 
Reflektivität wurde mit weißen Linien abgegrenzt und in die Darstellung weiterer Größen in 
Abbildung 5.15 übertragen. 

Die mittlere Schicht steht in Zusammenhang mit nahezu homogenen LDR-Werten von etwa 
— 25 dB. Diese Werte deuten auf abgerundete oder aufgrund der Variation der Geschwindigkei- 
ten stark taumelnde Eispartikel in dieser Schicht hin. Direkt darüber sind die Werte erheblich 
höher und liegen ungefähr bei -16 dB. Die Eispartikel in dieser Schicht sind mehrheitlich unre- 
gelmäßiger geformt und bestehen wahrscheinlich aus Dendriten oder Plättchen, die horizontal 
ausgerichtet sind. Die horizontale Ausrichtung ist daran zu erkennen, dass die LDR-Werte mit 
zunehmender Abweichung vom Zenit abnehmen. Ausgehend von diesen hohen Werten sinkt 
das LDR mit der Höhe auf -19 dB in einer Höhe von 4.5 km. In den grauen Bereichen liegen die 
Werte der Reflektivität im Cross-Kanal unterhalb der Detektionsgrenze und es kann kein LDR 
berechnet werden. Die LDR-Werte in der unteren Schicht, unterhalb 2 km, sind weitgehend 
homogen verteilt und liegen bei etwa -19 dB. 

Oberhalb der Scherungsschicht sind große Bereiche abwechselnd gering positiver sowie nega- 
tiver residualer Geschwindigkeit zu erkennen. Unterhalb dieser Schicht sind die Werte in der 
gleichen Größenordnung wie oberhalb der Schicht und liegen meist unter 0.1 m/s, lediglich in 
Höhen um 1 km sind minimal höhere Werte zu erkennen. Die räumliche Ausdehnung von re- 
sidualen Geschwindigkeiten gleicher Vorzeichen sind unterhalb der Scherungsschicht deutlich 
kleiner ausgeprägt. Innerhalb der Scherungsschicht sind die extremsten Werte der residualen 
Geschwindigkeit von bis zu ±0.8 m/s und periodische Strukturen der residualen Geschwindig- 
keit zu erkennen. Diese Strukturen ragen teilweise bis zu 400 m weit aus dieser Schicht heraus. 

Die Werte der spektralen Breite zeigen oberhalb der Scherungsschicht nur erhöhte Werte 
bei großen Zenitwinkeln, so dass sie von Schwankungen im Horizontalwind herrühren. Die- 
se Schwankungen sind in Höhenintervallen von 4.0km-4.4km, 3.4km-3.6km, sowie 3.1km- 
3.3 km erkennbar. In der Scherungsschicht sind auch im Bereich des Zenit über die gesamte 
Breite erhöhte Werte (bis zu 0.4 m/s) zu erkennen, was auf erhöhte Variationen der Vertikalge- 
schwindigkeiten der Eispartikel hindeutet. Unterhalb dieser Schicht ist die Verteilung der Werte 
wiederum homogener, die Werte sind gegenüber der Schicht oberhalb 3.1 km etwas erhöht. 
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Abbildung 5.15: RHI-Scans vom 08.01.2010, 13:44:31 UTC. Von oben nach unten sind die 
dargestellten Größen die Reflektivität in dBz und die residuale Geschwindigkeit in m/s. Grau 
unterlegt sind Bereiche, in denen Reflektivitäten vorhanden sind, aber keine Daten der jeweili- 
gen Größe vorliegen. Weiße Linien grenzen den Bereich ab, in denen wellenartige Strukturen 
in den Reflektivitäten vorliegen. 
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Abbildung 5.16: Wie Abb. 5.15. Von oben nach unten sind die dargestellten Größen das LDR 
in dB und die spektrale Breite in m/s. 
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Abbildung 5.17: Horizontalwind xim und Vertil<algescliwindigkeit der Streuer in Richtung 
der RHI-Scans (Azimut: 58 Grad) vom 08.01.2010, 13:44:31 UTC. Rot: Berechnete Werte aus 
RHI-IVIessungen des Wolkenradars IVIIRA36-S. Schwarz: Zeitliches und räumliches Mittel von 
je zwei Radiosondenaufstiegen. Graue Linien grenzen den Bereich ab, in denen wellenartige 
Strukturen in den Reflektivitäten vorliegen. 



In Abb. 5.17 sind die Profile der mittleren horizontalen Windgeschwindigkeit in Scanrichtung 
und der mittleren Vertikalgeschwindigkeit der Streuer auf Basis des gezeigten RHI-Scans darge- 
stellt. Die Trennung der verschiedenen Advektionsrichtungen der Luftmassen ist im Profil der, 
aus RHI-Messungen abgeleiteten, horizontalen Windgeschwindigkeit Um in Abbildung 5.17 
ebenfalls deutlich zu erkennen. Unterhalb einer Höhe von 3 km bewegen sich die Hydrometeo- 
re in Scanrichtung (58 °) mit bis zu 18 m/s (MIRA36-S) von links nach rechts und oberhalb mit 
bis zu 19 m/s von rechts nach links. 



Anhand des Profils der mittleren Vertikalgeschwindigkeit Wm ist festzustellen, dass die Sedi- 
mentationsgeschwindigkeit der Streuer von oben nach unten, ab einer Höhe von 5.2 km bis 
zu einer Höhe von 4.2 km um 0.4 m/s zunimmt. Es ist davon auszugehen, dass die Größe der 
Eispartikel in diesem Höhenintervall von oben nach unten zunimmt, da ebenfalls die Reflek- 
tivität höhere Werte annimmt. Der weitere Verlauf wird durch minimale Änderungen geprägt, 
bevor die Sedimentationsgeschwindigkeit in einer Schicht von 3.2 km bis 2.9 km, im oberen 
Bereich der Scherungsschicht, wieder um 0.4 m/s abnimmt. Obwohl die Werte der Reflekti- 
vität in diesem Höhenintervall im Vergleich zu den Werten knapp oberhalb wieder deutlich 
höher sind, nehmen die Sedimentationsgeschwindigkeiten geringere Werte an. In Verbindung 
mit den residualen Geschwindigkeiten ist zu erkennen, dass sich die Eispartikel in dieser Schicht 
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auch stellenweise nach oben bewegen. Nach einer erneuten Zu- und Abnahme in der gleichen 
Größenordnung, nimmt die Geschwindigkeit unterhalb 2 km wieder zu und erreicht Werte, die 
im Schneefall unterhalb 1 km, zwischen 0.9m/s und 1.1m/s liegen. 

Zusammenfassung Im gezeigten Fall liegt die Temperatur im gesamten Messbereich unter- 
halb von 0°C. Dadurch liegen die Hydrometeore ausschließlich in der Eisphase vor und eine 
Schmelzschicht ist nicht vorhanden. Die periodischen Fluktuationen der Reflektivitäten stehen 
in direktem Zusammenhang mit einer ausgeprägten Scherungsschicht. Dies sind ideale Bedin- 
gungen für interne Schwerewellen. Neben der Reflektivität ist die residuale Geschwindigkeit 
die einzige Größe, die deutlich sichtbare periodische Strukturen aufzeigt. Diese Fakten wer- 
den so interpretiert, als dass einerseits regelmäßige Auf- und Abwindbereiche initiiert werden 
und vertikale Austauschprozesse stattfinden. Andererseits wird eine Sperrschicht generiert, die 
wenige (niedrige Reflektivität) große (hohe Sedimentationsgeschwindigkeit) Eiskristalle und 
"Partikel innerhalb der Schmelzschicht von vielen (höhere Reflektivität) kleinen (niedrige Sedi- 
mentationsgeschwindigkeit) oberhalb der Schmelzschicht trennt. Die Eisteilchen innerhalb der 
Schmelzschicht besitzen zudem eine komplexere Form (hohes LDR), die einer Scheibe ähnelt, 
welche horizontal ausgerichtet ist (Abnahme des LDR mit zunehmender Abweichung vom Ze- 
nit). Die Eisteilchen darüber sind hingegen deutlich runder oder taumeln stark (niedriges LDR). 
Des weiteren deuten überlappende Bereiche der residualen Geschwindigkeiten auf ein gelegent- 
liches Eindringen der Partikel oberhalb der Sperrschicht in den darunter- und darüberliegenden 
Bereich hin (siehe Abb. 5.16). 

Im Gegensatz zu den dargestellten Messungen vom 06.10.2010 stehen die meisten Struktu- 
ren der Reflektivität in der Scherungsschicht in Zusammenhang mit Aufwinden (positive resi- 
duale Geschwindigkeiten). Messungen zeitlich vor und nach dem gezeigten RHI-Scan zeigen 
beständige räumlich periodische Strukturen in dieser Schicht. 

Die Erkenntnisse in diesem Kapitel bilden die Grundlage für das nächste Kapitel, in dem interne 
Schwerewellen in der Atmosphäre anhand von Messungen und der aus Messungen abgeleiteten 
residualen Geschwindigkeit diskutiert werden. 
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6 Detektion von Wellen in Wolken 

Wie bereits zum Schluss des Kapitels 5 gezeigt, können regelmäßige, wellenartige Strukturen 
in Wolken mit Hilfe der Reflektivität und der residualen Geschwindigkeit beobachtet werden, 
z.B. am Fall vom 08.01.2010. Mit der Annahme, dass es sich in diesem Fall um interne Schwe- 
rewellen handelt, werden in diesem Kapitel Möglichkeiten diskutiert, die Wellenlänge, die Pha- 
sengeschwindigkeit und die Ausbreitungsrichtung mit Hilfe der residualen Geschwindigkeit 
abzuschätzen. Der Begriff „Welle" bezieht sich in diesem Kapitel ausschließlich auf interne 
Schwerewellen. 



6.1 Wellendetektion 

Bisher gibt es kaum Arbeiten zur Messung von internen Schwerewellen in Wolken. Wenn le- 
diglich die direkt gewonnenen Größen, wie die Reflektivität, das LDR, die Dopplergeschwin- 
digkeit und die spektrale Breite berücksichtigt werden, ist es schwierig, Wellen zu detektieren. 
Dies wird erleichtert mit der Hinzunahme der residuale Geschwindigkeit. Damit können re- 
gelmäßige Strukturen in einem festen Höhenbereich Auf- und Abwindbereichen zugeordnet 
werden. In diesem Kapitel soU ein Verfahren gezeigt werden, in Wolken Wellen detektieren zu 
können. Hierfür werden folgende Annahmen getroffen: 

• Die Atmosphäre ist stabil geschichtet. 

• Die Ausbreitungsebene der Welle ist horizontal ausgerichtet. 

• Die Ausbreitungsrichtung der Welle aw weicht maximal 45° von der Richtung des Azi- 
mutwinkels a des RHI-Scans ab. 

Weiterhin sollen für die verwendeten Daten folgende Bedingungen gelten: 

• Es soUen nur Wellen detektiert werden, deren Wellenlängen Xw maximal der Länge des 
Messbereichs Mk in einer festen Höhe Zk entsprechen. Dabei ist k der laufende Index für 
jede Höhenstufe von unten nach oben, beginnend bei 1. 

• Der Messbereich in jeder Höhe soll mindestens 10 Radarstrahlen entsprechen, d.h. 
mindestens 10 Messwerte sollen für die Interpolation der Dopplergeschwindigkeit auf 
kartesischen Koordinaten, zur Berechnung der residualen Geschwindigkeit (Kap. 5.1), 
berücksichtigt werden. 
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Sind diese Voraussetzungen gegeben, wird für jede Höhenstufe mittels Fourier-Transformation 
die Frequenz der Abfolge von Auf- und Abwinden (schnell wechselnde Vorzeichen der residua- 
len Geschwindigkeit) berechnet. Dabei wird folgendermaßen vorgegangen: 

Zunächst werden in jeder Höhenschicht Fehlwerte (z.B. im Faltungsbereich) linear interpoliert. 
Vom Mittelwert und Trend innerhalb einer Höhenschicht müssen die Daten nicht mehr befreit 
werden, da dies schon während der Berechnung der residualen Geschwindigkeit geschehen ist. 
Außerhalb des Messbereichs werden die Daten auf eine feste Datenlänge von 4096 Werten 
mit dem Wert 0 ergänzt. Dabei nimmt die Länge des Messbereichs aufgrund des Messverfah- 
rens (RHI-Scan) mit zunehmender Höhe zu. Es fließen also mit zunehmender Höhe mehr be- 
rechnete Werte der residualen Geschwindigkeit in die Berechnung der Frequenzen ein. Durch 
die Ergänzung der Werte außerhalb des Messbereichs treten an den Rändern des Messbereichs 
Sprünge auf, die das Ergebnis der nachfolgenden Fouriertransformation verfälschen können. 
Diese Sprünge werden mit einem passenden Fensterfilter geglättet. Um ein möglichst deutliches 
Signal der Wellen zu bekommen und gleichzeitig Unterschwingungen der gefensterten Fou- 
riertransformation zu minimieren, wurde die nachfolgende Fensterfunktion (Hamming-Fenster) 
gewählt: 

i(;* = 0.54 + 0.46 cos (6.1) 

^'^ Mk 

Dabei ist rij = — Mfc/2, ...,Mk/2 — 1 der aktuelle Index des Eingangssignals (residuale Ge- 
schwindigkeit) in x-Richtung, in einer festen Höhe Zk, zentriert im Zenit des RHI-Scans. j ist 
der Zählindex in x-Richtung, beginnend bei 1. Mf, ist die Fensterbreite, d.h. die Anzahl der 
Messwerte, die abhängig von der Messung (RHI) einer Höhe Zk zugeordnet werden können 

(Mk^2k tan(45°) = 2 k). 

Die Fensterfunktion wird für jede Höhe Zk mit den Werten der residualen Geschwindigkeit der 
selben Höhe gefaltet und nachfolgend mittels dislcreter Fouriertransformation in den Frequenz- 
raum transformiert: 

N(k) 

E. 27r(3-l)(i-l) 
VRes,j,k-wl^e ' 4096 (6.2) 

i=l 

Mit der komplexen Konjugation wird dann die spektrale Signalstärke 



^l,k = ■^Res,/,fe • VRes,l,k (6.3) 



und die Amplitude 
sowie die Wellenlänge 



2 



Ak = -v^JSlk (6.4) 



4096 ^ 4096 
A^,-= Az = ^^30 m (6.5) 



berechnet. Dabei ist / = G N* A / < 4096/2}. 
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Die berechnete Wellenlänge basiert auf einer zweidimensionalen Messung (RHI-Scan). Da 
die Richtung der RHI-Scans nicht in allen Höhen der Ausbreitungsrichtung der detektierten 
Wellen entsprechen kann, muss die Wellenlänge noch korrigiert werden (Kapitel 6.2). Des- 
weiteren ist der Wellenlängenbereich der detektierbaren Wellen durch die zu Anfang des Ka- 
pitels gemachten Annahmen eingeschränkt. Der maximal mögliche Zenitwinkel (±45°) ei- 
nes RHI-Scans, die Abtastgenauigkeit der Messung und die Gitterauflösung des kartesischen 
Gitters limitieren den Wellenlängenbereich der detektierten Wellen. Deshalb wird zur Inter- 
pretation die minimale Wellenlänge, die in jeder Höhe detektiert werden kann, festgelegt auf 
^w,mm = 10 (2 Ax) = 10 [2 Zk tan (Ae/2)]. Dies entspricht zehn mal dem Messabstand zwi- 
schen zwei Zenitwinkeln (Ae), abhängig vom Abstand im Zenit. Die maximale Wellenlänge 
AvK,max wird begrenzt durch M^. 




Wellenlänge |km] Wellenlänge [km] 



Abbildung 6.1: Die mittlere Amplitude der fouriertransformierten residualen Geschwindigkeit 
abhängig von der Wellenlänge [km] und der Höhe [km] vom 08.01.2010 aus 96 RHI-Scans 
zwischen 13:30UTC und 14:17 UTC. Links: Der gesamte Höhenbereich in dem Hydrometeore 
detektiert wurden. Rechts: Ein höher aufgelöster Ausschnitt im Bereich der Scherungsschicht. 

Die Anwendung des gezeigten Verfahrens zur Abschätzung der Wellenlänge in Abhängigkeit 
der Höhe wird am Beispiel vom 08.01.2010 veranschaulicht. Hierfür werden RHI-Scans nach 
dem oben beschriebenen Verfahren ausgewertet und die Amplitude Ai^^ über eine Stunde 
gemittelt. Das arithmetische Mittel der Amplituden von 96 RHI-Scans des Zeitraums von 
13:30 UTC-14:17 UTC ist in Abb. 6.1 dargestellt. Auffällig sind drei Bereiche hoher Amplitu- 
den, unterhalb 0.5 km, in einem Höhenbereich um 1 km und zwischen 2.2 km und 3.4km (Abb. 
6.1, links). Für die Bereiche unterhalb 1.5 km ist die Zuordnung von Wellenlängen schwie- 
rig, da die maximalen Amplituden am Rand des detektierbaren Bereichs innerhalb des je- 
weiligen Höhenabschnitts liegen und somit wenig Daten in die Fourier-Transformation ein- 
gehen. Deshalb wird im Folgenden hauptsächlich der Höhenabschnitt zwischen 2.2 km und 
3.4km diskutiert, in dem sich eine ausgeprägte Scherungsschicht befindet (siehe Kap. 5.2.2). 
Dieser Höhenabschnitt ist in Abb. 6.1, rechts, detaillierter dargestellt. Es sind drei Berei- 
che erhöhter Amplituden zu erkennen, die sich von der Umgebung abgrenzen. Im ersten 
Höhenbereich von 2.5 km bis 2.8 km liegt das jeweilige Maximum der Amplituden (abhängig 
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von der Höhe) bei einer Wellenlänge von etwa 1.7 km und die höchste Amplitude ist in einer 
Höhe von ca. 2.7km zu finden. Der zweite Bereich erhöhter Amplituden befindet sich in ei- 
nem Höhenbereich zwischen 2.8 km und 3.1 km und die Maxima der Amplituden können einer 
Wellenlänge von ca. 1.6 km zugeordnet werden. Der dritte Höhenabschnitt erhöhter Amplitu- 
den ist zwischen 3.15 km und 3.35 km angedeutet. In diesem Höhenabschnitt sind die Maxima 
der Amplituden für kürzere Wellenlängen von ca. 0.9 km zu erkennen. 

Die gezeigten abgegrenzten Bereiche erhöhter Amplituden deuten auf das Vorhandensein von 
Wellen hin. Die berechneten Wellenlängen aus zweidimensionalen Daten sind dabei ein erstes 
Maß für die „wahre" Wellenlänge. Für die Abschätzung der „wahren" Wellenlänge der drei- 
dimensionalen Welle ist die Kenntnis der Ausbreitungsrichtung nötig. Eine erste Abschätzung 
der Ausbreitungsrichtung wird mit Hilfe des horizontalen Grundstroms, bzw. des mittleren Ho- 
rizontalwinds vorgenommen und in Kapitel 6.2 diskutiert. 

Mit den fouriertransformierten residualen Geschwindigkeiten konnten Periodizitäten, mögliche 
Wellen, in Höhen zwischen 2.5 km und 3.35 km detektiert werden. Die Eigenschaften der Peri- 
odizitäten werden im Folgenden untersucht und verglichen mit den Eigenschaften von internen 
Schwerewellen. 



6.2 Eigenschaften der detektierten Wellen 

Ziel ist es, Wellenparameter wie die Wellenlänge und die Phasengeschwindigkeit, aus Messun- 
gen von RHI-Scans mit einem Wolkenradar abzuschätzen. Für die anschließende Diskussion 
sollen Referenzwerte dieser Parameter auf der Basis theoretischer Grundlagen aus Radioson- 
denmessungen abgeleitet werden. 

6.2.1 Abschätzung der Phasengeschwindigkeit und Wellenlänge mittels Radiosondenda- 
ten 

Voraussetzung für das Auftreten interner Schwerewellen ist eine thermisch geschichtete At- 
mosphäre. Ist die Atmosphäre stabil geschichtet, erfahren Luftpakete bei einer Auslenkung in 
großen Höhen eine rücktreibende Kraft durch die Schwerkraft. Wenn man von Reibungsein- 
flüssen absieht, folgt aus einer vereinfachten Form der Vertikalbeschleunigung des Luftpakets 
die S ch wingung s diff erentialgleichung : 

^ + N'^w = 0 (6.6) 
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mit der Vertikalgeschwindigkeit w und der Brunt-Väisälä-Frequenz für trockene Luft 




(6.7) 



Hierbei ist g = 9.81 m/s die Schwerebeschleunigung. Die Brunt-Väisälä-Frequenz ist ein Sta- 
bilitätsmaß der Atmosphäre und abhängig von der potentiellen Temperatur trockener Luft: 

\p 

Diese Temperatur würde erreicht werden, wenn ein Luftpaket adiabatisch vom Druck p auf 
Normaldruck {po = 1000 hPa) gebracht würde. Ri = 287 J/(kg K) ist die spezifische Gaskon- 
stante für trockene Luft und Cp = 1004 J/(kg K) die spezifische Wärmekapazität bei konstantem 
Druck. 

Mit den Anfangsbedingungen für kleine, beschleunigungsfreie Vertikalgeschwindigkeiten 
w{to) = u'o und dw{to)/dt = 0 zum Zeitpunkt to = 0 ergeben sich Lösungen der Schwin- 
gungsdifferentialgleichung für folgende Fälle: 

{Wo cosh{\N\t) für < 0 
Wo für = 0 (6.9) 
Wocos{Nt) für N^>0 

Für < 0 entfernt sich ein Luftpartikel mit zunehmender Zeit von der Ausgangslage. Der in- 
stabile Fall liegt vor und das Luftpartikel ist bei einer virtuellen Auslenkung nach oben wärmer 
als ihre Umgebung und ein Auftrieb wird wkksam. Für = 0 liegt der neutrale Fall vor und 
das Luftpartikel verharrt in der Ausgangslage. Für > 0, im stabilen Fall, findet eine har- 
monische Schwingung um die Ausgangslage statt. D.h. bei einer virtuellen Auslenkung nach 
oben ist das Luftpartikel kälter und schwerer als ihre Umgebung und die rücktreibende Kraft 
(Schwerebeschleunigung) wird wirksam. Im umgekehrten Fall, bei einer Auslenkung nach un- 
ten, kehren sich die Verhältnisse um und der Auftrieb ist die bestimmende Kraft. Somit muss 
für interne Schwerewellen eine stabile Schichtung vorhanden sein und für die Brunt-Väisälä- 
Frequenzen gelten: iV^ > 0 (36 /dz > 0, Gl. 6.7). 

Die horizontale Phasengeschwindigkeit cw einer internen Schwerewelle kann unter Verwen- 
dungung der Brunt-Väisälä-Frequenz abgeschätzt werden zu (siehe u.a. Houze, 1993): 
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cw = uo± — ^ uo± — (6.10) 

k 2 TT 

Dabei ist uq der horizontale Grundstrom, k die horizontale Wellenzahl und Xw die horizontale 
Länge der Welle. Dies ist eine vereinfachte Form der Berechnung der Phasengeschwindigkeit 
unter der Annahme, dass die vertikale Ausbreitungskomponente der Welle viel kleiner ist als die 
horizontale. Dies trifft zu, wenn sich die Strömung im hydrostatischen Gleichgewicht befindet. 

Die Phasengeschwindigkeit cw einer internen Schwerewelle wird im Folgenden zunächst mit 
Daten von Radiosondenaufstiegen berechnet und im weiteren Verlauf in Relation gesetzt mit 
einer Phasengeschwindigkeit, die rein aus Wolkenradar-Messungen auf der Basis von RHI- 
Scans abgeschätzt wird. 

Um die Phasengeschwindigkeit mit Hilfe der Gl. 6.10 abschätzen zu können, muss das Profil 
der potentiellen Temperatur und das Profil des horizontalen Grundstroms bekannt sein. Da das 
Wolkenradar keine Temperaturen und Drücke misst, werden als Näherung die Werte operatio- 
neil messender Radiosonden in Stuttgart/Schnarrenberg (ST) und Idar-Oberstein (10) herange- 
zogen. Aus den dort gemessenen Profilen der Temperatur T und des Drucks p werden Profile 
der potentiellen Temperatur (Gl. 6.8) und der Brunt-Väisälä-Frequenz (Gl. 6.7) berechnet und 
für Messungen vom 08.01.2010, 12:00 UTC in Abb. 6.2 dargestellt. 




270 275 280 285 290 295 300 305 310 315 5 10 15 20 25 30 35 

Potentielle Temperatur [K] Brunt-Väisälä-Freauenz [1 0"^ s^l 



Abbildung 6.2: Berechnete potentielle Temperatur (links) und Brunt-Väisälä-Frequenz für 
trockene Luft (rechts) auf der Basis von Radiosondenmessungen vom 08.01.2010, 12:00 UTC 
am Standort ST (schwarz) und IG (grün). 

Die potentielle Temperatur (Abb. 6.2) nimmt für beide Profile (ST, 10) oberhalb 500 m mit der 
Höhe zu und die Atmosphäre ist trockenstabil geschichtet. Die stabile Schichtung zeigt sich 
auch in den daraus berechneten positiven Frequenzen der Brunt-Väisälä-Frequenz für trockene 
Luft (A^ > 0 1/s). In den Profilen am Standort ST sind zwei Bereiche auffällig: Die stabilste 
Schichtung befindet sich knapp unterhalb 1 km, erkennbar an der deutlichsten Zunahme der po- 
tentiellen Temperatur und den höchsten Brunt-Väisälä-Frequenzen. Ein weiterer Bereich hoher 
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Brunt-Väisälä-Frequenzen, in dem sich auch die ausgeprägte Scherungsschicht befindet, ist in 
einem Höhenintervall zwischen 2.2 km und 3.7 km zu erkennen und wird lediglich in einem 
Höhenbereich von etwa (2.7 ± 0.1 )km, direkt oberhalb einer Inversion (Abb. 5.13) unterbro- 
chen. Die niedrigen Brunt-Väisälä-Frequenzen im Höhenbereich von (2.7 ± 0.1)km stehen in 
Zusammenhang mit einer nahezu isothermen Schicht, die auch an der geringen Zunahme der 
potentiellen Temperatur mit der Höhe sichtbar wird. In den Profilen am Standort 10 ist ebenfalls 
eine stabile Schichtung zu erkennen, die unterhalb 1 km beginnt, sich aber bis knapp 1.5 km 
erstreckt. Die weiteren Höhenintervalle stabiler Schichtung sind verständlicherweise verbunden 
mit den Inversionsschichten (siehe Abb. 5.13). Sie liegen oberhalb der stabilen Bereiche, die im 
Profil der Brunt-Väisälä-Frequenz am Standort ST angezeigt werden und befinden sich knapp 
unterhalb 3 km, zwischen 3.1 km und 4 km, sowie knapp oberhalb 4.5 km. 

Um nun die horizontale Phasengeschwindigkeit cw (Gl. 6.10) abzuschätzen, ist die Kenntnis 
des mittleren Windfeldes am Radarstandort notwendig. Zur Verfügung stehen Messungen vom 
Wolkenradar selbst (PPI-Scan von 13:05 UTC und anschließende Berechnung des Windfeldes 
mittels VAD-Algorithmus (Kap. 2.4.2)) und von Radiosondenaufstiegen an den Standorten ST 
und lO (12:00 UTC). Diese sind in Abb. 6.3 dargestellt. Auf die auffälligsten Merkmale der Pro- 
file der Windgeschwindigkeit und Windrichtung anhand Messungen des Wolkenradars in Abb. 
6.3 wurde bereits in Kap. 5.2.2 eingegangen. Unterschiede im Verlauf der Profile sind aufgrund 
der verschiedenen Messzeitpunkte und Standorte zu erkennen. Messungen des Wolkenradars 
und des Sondenaufstiegs am Standort ST zeigen beispielsweise den Beginn der ausgeprägten 
Rechtsdrehung des Windes in einer Höhe knapp oberhalb 2 km, während Messungen des Son- 
denaufstiegs am Standort 10 dies erst knapp unterhalb 2 km zeigen. 
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Abbildung 6.3: Profile der Windgeschwindigkeit (links) und Windrichtung (rechts) aus Mes- 
sungen des Wolkenradars MIRA36-S (rot) und der Radiosondenaufstiege am Standort ST 
(schwarz) und IG (grün). 

Zusammenfassend bleibt festzustellen, dass die Höhenabschnitte, in denen Wellen vermutet 
werden (Abb. 6.1), gekennzeichnet sind durch eine starke Richtungsscherung, sowie eine star- 
ke Geschwindigkeitsscherung. Die höchsten Amplituden der Wellen sind im Höhenabschnitt 
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zwischen 2.5 km und 2.8km zu erkennen (Abb. 6.1). In diesem Höhenbereich ist eine ausge- 
prägte monotone Rechtsdrehung des Windes und eine Vorzeichenänderung des Geschwindig- 
keitsgradienten des Windes, gefolgt von einer starken Abnahme der Geschwindigkeit mit der 
Höhe zu erkennen. Die starken Scherungen, die zuvor beschriebene stabile, nahezu isother- 
me Schichtung in Abschnitten der Scherungsschicht und das Vorhandensein von regelmäßigen 
angeordneten Strukturen in der Reflektivität und in der residualen Geschwindigkeit innerhalb 
einer ausgeprägten Scherungsschicht sprechen für die Annahme, dass interne Schwerewellen 
vorliegen. 

Zur Abschätzung der Wellenlängen und Phasengeschwindigkeiten nach Gl. 6.10 wird das 
Hauptaugenmerk auf die in Kap. 6.1 genannten drei Höhenabschnitte gelegt, in denen peri- 
odische Fluktuationen der residualen Geschwindigkeit detektiert wurden. Diese befinden sich 
in den Höhenintervallen (1) 2.5 km-2.8 km, (2) 2.8 km-3.1 km und (3) 3.15 km-3.35 km. Die für 
diese Höhenintervalle repräsentativen Brunt-Väisälä-Frequenzen und Windgeschwindigkeiten 
auf Basis der gezeigten Profile sind Tab. 6.1 zu entnehmen. 





2.5 km-2.8 km 


2.8 km-3.1 km 


3.15 km— 3.35 km 


NsT [10-3 s-1] 


16.3 ± 1.6 


16.7 ±0.2 


17.1 ±0 


Nio [10-^s-i] 


12.2 ±3.1 


17.6 ±3.1 


20.6 ± 1.3 


Mppi [m/s] 


10.4 ±1.1 


4.2 ±2.0 


6.3 ± 1.3 


«ST [m/s] 


8.8 ±0.1 


9.1 ±0.8 


11.3 ±0.3 


uio [m/s] 


7.6 ± 0.4 


4.4 ±2.1 


4.9 ± 1.3 


aw,ppi [°] 


108 ± 11 


133 ±5 


180 ±6 


Olw,ST [°] 


137 ±7 


170 ± 19 


203 ±0 


aw,io [°] 


76 ± 1 


97 ±24 


180 ± 17 



Tabelle 6.1: Mittlere Werte und deren Standardabweichung von der horizontalen Windge- 
schwindigkeit «ppi und Windrichtung aw,ppi zum Zeitpunkt 14:18 UTC (MIRA36-S: PPI-Scan). 
IVlittlere Werte und deren Standardabweichung von der Brunt-Väisälä-Frequenz N, der horizon- 
talen Windgeschwindigkeit u und Windrichtung zum Zeitpunkt 12:00 UTC (Radiosonden: 
ST, 10). 

Für eine erste Abschätzung der Phasengeschwindigkeit wird diese als Funktion der Wel- 
lenlänge berechnet (Gl. 6.10) und ist in Abb. 6.4 dargestellt. Hierfür wurden drei verschie- 
dene Brunt-Väisälä-Frequenzen, = 12 X 10-3 s-^ (rot), A/" = 16 X 10-3 s-^ (grün), = 
20xl0~3 (blau)) und zwei verschiedene Horizontalgeschwindigkeiten von u = 5 m/s und 
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u = 10 m/s fest gewählt, die die Variabilität der Messwerte aus Tab. 6.1 wiedergeben. Die in 
Abb. 6.4 gezeigten Werte sollen die zu erwarteten Phasengeschwindigkeiten als Funktion der 
Wellenlängen im Bereich der Scherungsschicht repräsentieren. 




i I \ \ I i i i i I \ i i i 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 

Wellenlänge [km] 

Abbildung 6.4: Berechnung der Phasengeschwindigkeit cw interner Schwerewellen nach Gl. 
6.10 als Funktion der Wellenlänge für die Brunt-Väisälä-Frequenzen N = 12xl0~^s"^ (rot), 
= I6xl0~'^s"^ (grün), N = 20x10"^ s"^ (blau) und die Horizontalgeschwindigkeiten von 

u = 5 m/s und u= 10 m/s (für Xw = 0 ist cw = u). 



6.2.2 Abschätzung der Phasengeschwindigkeit und Wellenlänge mittels Wolkenradar 

Eine erste Abschätzung der Phasengeschwindigkeit in Scan-Richtung mittels Wolkenradar lie- 
fert die Kreuzkorrelation der residualen Geschwindigkeiten jedes gleitenden Höhenintervalls 
über drei Höhenstufen zwischen zwei zeitlich aufeinanderfolgenden (At = 30 s) RHI- 
Scans. Dabei wird die horizontale Verteilung der residualen Geschwindigkeiten in einem 
Höhenintervall zum Zeitpunkt t^+i um bis zu 25 Gitterpunkte in x-Richtung (dx = 30 m) 
jeweils nach rechts und links verschoben. Die maximale Anzahl der Gitterpunkte wurde so 
gewählt, dass kein Aliasing der Verlagerungsgeschwindigkeit stattfinden kann. Ein Aliasing der 
Geschwindigkeit wird erst möglich, wenn sich die Strukturen mit mindestens 25 m/s verlagern 
würden. Anhand der zu erwartenden abgeschätzten Phasengeschwindigkeiten für Wellenlängen 
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kleiner 3 km (Abb. 6.4) wird diese Aliasing-Geschwindigkeit nicht überschritten und die Be- 
grenzung der Verschiebung auf 25 Gitterpunkte ist ausreichend. Für jeden Verschiebungsschritt 
dx wird der Korrelationskoeffizient aus dem horizontal verschobenen Höhenintervall zum Zeit- 
punkt ii-i-i und dem unverschobenen zum Zeitpunkt U berechnet. Der höchste Korrelationskoef- 
fizient gibt dann die horizontale Verlagerung der Strukturen relativ zur Richtung des RHI-Scans 
an. Zur Reduktion des Rauschens erfolgt zunächst eine Mittelung. Die Verschiebung für den 
Zeitpunkt ti und der Höhe zi wird durch den maximalen Korrelationskoeffizient gdU, ^i) in ei- 
nem Zeitabschnitt von ti_2 bis ^»+2 und einem Höhenabschnitt zwischen Zi^^ und 2:1+2 bestimmt. 
Aus der Verschiebung in Gitterpunkten, dem Gitterabstand und der Zeitspanne aufeinander 
folgender RHI-Scans kann die horizontale Verlagerungsgeschwindigkeit in Scan-Richtung be- 
stimmt werden. Mit der Annahme, dass sich die Strukturen näherungsweise mit der gleichen 
Geschwindigkeit wie die mögfichen internen Schwerewellen in Richtung der Wellen bewe- 
gen, wird die Verlagerungsgeschwindigkeit der Strukturen cw,r (R steht für Radar) als eine 
Abschätzung der Phasengeschwindigkeit aufgefasst. 

Die horizontale Verlagerungsgeschwindigkeit ist für den Zeitraum von 13:30 UTC bis 
14:17 UTC für einen Höhenbereich zwischen 2.2 km und 3.3 km dargestellt in Abb. 6.5, 
oben. Bereiche negativer Verlagerungsgeschwindigkeiten sind meist unterhalb einer Höhe von 
2.9 km zu finden, d.h. die Strukturen bewegen sich in diesem Fall in Richtung des mittleren Ho- 
rizontalwindes (siehe Abb. 6.6), in Scan-Richtung (RHI-Scan) von links nach rechts. Oberhalb 
2.9 km bewegen sich die Strukturen einerseits in Richtung des Horizontalwindes. Ein Zeitab- 
schnitt, für den dies besonders zutrifft, ist zwischen 14:10 UTC und 14: 17 UTC zu erkennen. 
Andererseits sind auch Bereiche erkennbar, in denen sich die Strukturen zeitweise deutlich 
entgegen dem Horizontalwind bewegen, besonders im Zeitintervall zwischen 13:48 UTC und 
13:54 UTC. In dieser Zeit bewegen sich die Strukturen mit 6 m/s bis 8 m/s von links nach rechts, 
während die Streuer durch dem Horizontalwind mit etwa 3 m/s von rechts nach links bewegt 
werden. Stabile, gleichmäßige Strukturen in den residualen Geschwindigkeiten sind anhand 
hoher Korrelationskoeffizienten (Abb. 6.5, unten), Qc{ti, Zi) > 0.7 zu erkennen, die besonders 
im Zeitraum 13:48 UTC und 13:54 UTC in einer Höhe von 2.5 km bis 2.7 km zu finden sind. 
Bei der Betrachtung des gesamten Zeitraums fällt auf, dass stabile, sich verlagernde Strukturen 
der residualen Geschwindigkeit immer wieder in Höhen zwischen 2.30 km und 2.75 km sowie 
zwischen 3.05 km und 3.30 km zu sehen sind. Instabile Strukturen der residualen Geschwin- 
digkeit sind dagegen hauptsächlich zwischen 2.75 km und 3.05 km augenscheinfich. Vergleicht 
man dieses Ergebnis, speziell die letzten Minuten vor 14:17 UTC, mit den Windprofilen (Abb. 
6.3) des Wolkenradars von 14:18 UTC, befinden sich instabile Strukturen der residualen Ge- 
schwindigkeit in einem Höhenabschnitt minimaler Geschwindigkeiten des Horizontalwindes 
und gleichzeitig im Bereich schnell wechselnder Windrichtungen und somit vermutlich in ei- 
nem turbulenzgeprägten Höhenabschnitt. Die stabilen Strukturen stehen hingegen in Zusam- 
menhang mit einer ausgeprägten Rechtsdrehung des Windes und einem starken Gradienten der 
Windgeschwindigkeit. 
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Abbildung 6.5: Oben: Aus den residualen Geschwindigkeiten (abgeleitet aus Dopplermes- 
sungen mittels RHI-Scans vom 08.01.2010) berechnete Werte der horizontalen Verlagerungs- 
geschwindigkeit cw,R von Strukturen hoher und niedriger residualer Geschwindigkeiten in 
m/s zwischen 13:30 UTC und 14:17 UTC für das Höhenintervall 2.2 km bis 3.3 km. Unten: Der 
für die Berechnung der Verlagerungsgeschwindigkeit zugrundeliegende Korrelationskoeffizient 
Qciti,zi) (Siehe Text). 
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Abbildung 6.6: Aus Dopplermessungen (RHI-Scans) abgeleiteter Horizontalwind u,„ in m/s 
vom 08.01 .201 0 zwischen 1 3:30 UTC und 14:17 UTC, für das Höhenintervall 2.2 km bis 3.3 km. 
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Abbildung 6.7: Wie Abb. 6.1 für das Zeitintervall zwischen 1 3:48 UTC bis 1 3:54 UTC. 
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Abbildung 6.8: Wie Abb. 6.1 für das Zeitintervall zwischen 14:10 UTC bis 14:1 7 UTC 

Die besonderen Zeitabschnitte, in denen sich Strukturen periodischer residualer Geschwindig- 
keiten oberhalb 3 km einmal mit dem Horizontalwind und einmal entgegengesetzt verlagern, 
werden gesondert betrachtet. Die Wellenlänge wird analog zur Detektion von Wellen berech- 
net (Kap. 6.1) für die beiden Zeitfenster von a) 13:48 UTC bis 13:54 UTC und b) 14:10UTC 
bis 14:17 UTC. Für beide Zeitfenster sind zwei Bereiche erhöhter Amplituden zu erkennen. 
Der erste befindet sich im Höhenabschnitt zwischen 2.5 km und 2.8 km (Abb. 6.7, 6.8). Die 
maximalen Amplituden in jeder Höhe in diesem Höhenabschnitt können dabei Wellenlängen 
von ca. 1.6 km zugeordnet werden. Die dafür verantwortlichen Periodizität der residualen Ge- 
schwindigkeit wird als ,, Welle 1" bezeichnet. Der zweite Bereich befindet sich direkt darüber 
und erstreckt sich von 2.8 km bis 3.1 km. Die maximalen Amplituden sind bei Wellenlängen 
zwischen 1.5 km und 1.8 km zu erkennen. Diese Periodizität wird als „Welle 2" bezeichnet. Für 
beide Wellen sind die Amplituden im Zeitfenster b) schwächer ausgeprägt. Zusätzlich sind im 
Zeitfenster b) zwei weitere Bereiche erhöhter Amplituden in den Höhenabschnitten von 3.0 km 
bis 3.1 km und 3.15 km bis 3.35 km für Wellenlängen von etwa 0.9 km zu sehen. Der Bereich 
der erhöhten Amplituden, der sich deutlicher abzeichnet, befindet sich im Höhenabschnitt von 
3.15 km bis 3.35 km und wird als „Welle 3" bezeichnet. 

Um die anschließende Diskussion übersichtlich zu gestalten, werden für die beiden beschriebe- 
nen Zeitfenster Höhenabschnitte und Wellenlängenbereiche betrachtet, die den drei definierten 
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Wellen zugeordnet wurden. Die ermittelten Größen für die jeweilige Welle im jeweiligen Zeit- 
fenster ist in Tab. 6.2 aufgeführt. 

Mit der Annahme, dass sich die Wellen in Richtung des horizontalen Grundstroms bewegen, 
wird die Wellenlänge Ai^ und die Phasengeschwindigkeit cw für die drei Wellenbereiche Welle 
1, Welle 2 und Welle 3 abgeschätzt: 



Xw = Ah/,rhi I cos(a - a^^ppi) \ 

cw,Rm (6.11) 

Cw — 7 r 

cos(ü; - a^^ppi) 



Der Azimuth der RHI-Scans ist a = 58°. 



RHI-Scans 


13:48UTC-13:54UTC 


14:10UTC-14:17UTC 


Welle 1 


Welle 2 


Welle 3 


WeUe 1 


Welle 2 


Welle 3 


Xw,RBi [km] 


1.6 ±0.2 


1.7 ±0.2 




1.8 ±0.2 


1.8 ±0.2 


0.9 ±0.1 


ci^,RHi [m/s] 


-8.1 ± 1.5 


-6.1 ±2.2 


-7.2 ± 1.6 


-9.3 ±1.7 


-2.7 ±4.3 


3.6 ±0.8 


Qc 


0.8 ±0.2 


0.6 ±0.2 


0.7 ±0.2 


0.8 ±0.1 


0.6 ±0.1 


0.8 ±0.1 


Um [m/s] 


-6.3 ± 1.8 


-1.1 ±0.8 


3.6 ±0.8 


-6.4 ± 1.8 


-1.0 ±0.7 


3.9 ±0.5 



Tabelle 6.2: Mittlere Werte und deren Standardabweichung von der Wellenlänge Xw,Rm, ho- 
rizontale Verlagerungsgeschwindigkeit von räumlichen Fluktuationen in den residualen Ge- 
schwindigkeiten cw,RBi, Korrelationskoeffizient der Verlagerungsgeschwindigkeit Qc und ho- 
rizontale Windgeschwindigkeit Um für die Zeiträume 13:48 UTC-1 3:54 UTC und 14:10 UTC- 
14:17 UTC (RHI-Scans). 



Welle 1 Die Wellenlänge in Scanrichtung beträgt im Mittel über das jeweilige Zeitfenster 
><w,Rm = (1.6 ± 0.2) km (Tab. 6.2, erstes Zeitfenster) und Xw,Rm = (1-8 ± 0.2) km (Tab. 6.2, 
zweites Zeitfenster). Der Mittelwert der Verlagerungsgeschwindigkeit, die nach Annahme der 
Phasengeschwindigkeit entsprechen soll, ist cv^.rhi = (—8, 1 ± 1.5) m/s (Tab. 6.2, erstes Zeit- 
fenster) und cw,Rm = (—9.3 ± 1.7) m/s (Tab. 6.2, zweites Zeitfenster). Die Windrichtung ist in 
dieser Schicht zum Zeitpunkt 14: 18 UTC, eine Minute nach dem Ende des ersten Zeitfensters, 
ctw,ppi = (108 ± 11)° (Tab. 6.1). Für den Höhenabschnitt, in dem Welle 1 delektiert wurde, ist 
die Wellenlänge in Ausbreitungsrichtung Xw = (1.0 ± 0.6) km für das erste Zeitfenster, bzw. 
Xw = (1-1 ± 0.4) km für das zweite (Gl. 6.11). Die Phasengeschwindigkeit in Ausbreitungs- 
richtung ist Cw = (12.6 ± 5.7) m/s im ersten Zeitfenster und cw — (14.8 ± 6.5) m/s im zweiten 
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Zeitfenster (Gl. 6.1 1). Im ersten Höhenabschnitt, in der sich Wellel befindet, ist die mittlere ho- 
rizontale Geschwindigkeit «m.ppi = (10.4 ± 1.1) m/s (Tab. 6.1). Etwa zwei Stunden zuvor (Ra- 
diosondenaufstieg) nimmt die Brunt-Väisälä-Frequenz Werte von (16.3 ± 1.6)xl0~^ s~^ (ST) 
und (12.2 ± 3.1) x 10"^ s-^ (10) an (Tab. 6.1). 

Nach den theoretischen Vorgaben auf der Basis der Radiosondenmessungen (Abb. 6.4) müsste 
Welle 1 für diese Werte mit etwa 12 m/s in Richtung des horizontalen Grundstroms, bzw. mit 
etwa 8 m/s entgegen des Grundstroms laufen. Da sich die Strukturen in diesem Höhenabschnitt 
mit dem Grundstrom verlagern, ist der erste Bereich der Entscheidende. D.h. die aus Messungen 
des Wolkenradars abgeleitete Verlagerungsgeschwindigkeit stimmt im ersten Zeitfenster in gu- 
ter Näherung mit der aus Radiosondenmessungen abgeleiteten Phasengeschwindigkeit überein 
und zeigt im zweiten Zeitfenster eine Abweichung von ca. 3 m/s, die wiederum im Bereich der 
Standardabweichung von cw liegt. 



Welle 2 Für die Welle 2 wird gleichermaßen vorgegangen und die Phasengeschwindigkeit in 
Ausbreitungsrichtung ist cvf^ = (23.6±18.2) m/s (erstes Zeitfenster) und c^^ — (10.4±22.5) m/s 
(zweites Zeitfenster), während die Wellenlänge Avi^ = (0.4 ± 0.2) km (erstes Zeitfenster) 
und Xw = (0.5 ± 0.2) km (zweites Zeitfenster) in Richtung des Horizontalwinds beträgt. 
Mit einem Horizontalwind von u = (4.2 ± 2.0) m/s und einer Brunt-Väisälä-Frequenz von 
NsT = (16.7 ± 0.2)xl0-^s-\ bzw. Njo = (17.6 ± 3.1)xl0-^s-^ kann die zu erwartende 
Phasengeschwindigkeit abgeschätzt werden zu 5.5 m/s in Windrichtung und 2.5 m/s entgegen 
der Windrichtung. 

In diesem Fall zeigt sich, dass sich die Welle in diesem Höhenabschnitt nicht mit dem Horizon- 
talwind dieses Höhenabschnitts ausbreitet. Dabei ist zu beachten, dass die Standardabweichung 
der Wellenlänge und Phasengeschwindigkeit in dieser Schicht aufgrund der starken Scherung 
sehr große Werte annimmt. Dennoch ist die Abweichung der abgeschätzten Phasengeschwin- 
digkeit aus Messungen des Wolkenradars gegenüber den zu erwartenden Phasengeschwindig- 
keiten im zweiten Zeitfenster groß. Mit der hypothetischen Annahme, dass im ersten Zeitfen- 
ster Welle 2 von Welle 1 beeinflusst wird, wäre die Annahme der Ausbreitungsrichtung der 
Welle 2 mit der Windrichtung falsch und sie müsste sich in Ausbreitungsrichtung der Welle 
1 bewegen. Dies zeigt sich zum einen darin, dass die Verlagerungsrichtung relativ zur Scan- 
richtung der RHI-Scans im Höhenabschnitt der Welle 2 und darüber hinaus, bis zu einer Höhe 
von 3.3 km, jener von Welle 1 entspricht (Abb. 6.5, oben). Zum anderen ist dies im Verhal- 
ten der Änderungen der Wellenlängen zu erkennen, die im ersten Zeitfenster (13:48UTC bis 
13:54 UTC) synchron verlaufen (Abb. 6.9). Für die in Abb. 6.9 dargestellte zeitliche Änderung 
der Wellenlängen und Amplituden werden im ersten Schritt für jeden RHI-Scan im Abstand 
von 30 s die Wellenlängen und Amplituden jeder Höhenstufe analog zum präsentierten Ver- 
fahren in Kap. 6.1 berechnet. Im zweiten Schritt wird daraus die mittlere Wellenlänge und 



100 



6 Detektion von Wellen in Wolken 



die mittlere Amplitude über alle Höhen im jeweiligen Höhenabschnitt für den jeweiligen Wel- 
lenlängenbereich von Welle 1, Welle 2 und Welle 3 ermittelt. 

Welle 3 Welle 3 ist im ersten Zeitabschnitt nicht zu erkennen, kann aber im zweiten Zeitab- 
schnitt identifiziert werden. Analog zu den zuvor gezeigten Wellen ergibt sich eine Phasenge- 
schwindigkeit von cw = (6.8 zt 2.7) m/s und eine Wellenlänge von Xw = (0.5 ± 0) km, wenn 
sich die Welle in Richtung des horizontalen Windfeldes verlagert. Mit einem Horizontalwind 
von (6.3 dz 1.3) m/s und einer Brunt-Väisälä-Frequenz von 17.1x10"^ s~^ ergibt sich eine zu 
erwartende Phasengeschwindigkeit von ca. 7.5 m/s in Richtung des Horizontalwindes und ca. 
5 m/s entgegen der Windrichtung. Damit liegt die Verlagerungsgeschwindigkeit der Welle 3 im 
Bereich der zu erwartenden Phasengeschwindigkeit. 

Interessant ist, dass Welle 3 die höchsten Amplituden aufweist (Abb. 6.9), wenn die Verla- 
gerungsgeschwindigkeit im Höhenabschnitt dieser Welle näherungsweise der mittleren Wind- 
geschwindigkeit entspricht (Abb. 6.6). Diese Beziehung zusammen mit den synchronen 
Änderungen der Wellenlänge und der Amplitude von Welle 1 und Welle 2 (Abb. 6.9) führt 
zu folgenden Schlussfolgerung: Welle 1 und Welle 3 werden innerhalb ausgeprägter Scherun- 
gen erzeugt und breiten sich näherungsweise in Richtung des mittleren Horizontalwinds im 
jeweiligen Höhenabschnitt aus. Welle 2 hingegen wird von Welle 1 getriggert. 




Welle 3 



13:40 13:50 14:00 14:10 13:40 13:50 14:00 14:10 

Tageszeit ( UTC) Tageszeit ( UTC) 

Abbildung 6.9: Zeitreihen der Wellenlänge und Amplitude dreier ausgewählter periodischer 
Fluktuationen der residualen Geschwindigkeit (siehe Text) 

6.2.3 Kelvin-Helmholtz-Instabilitäten 

Der starke Anstieg der Amplitude der Welle 1 und der Welle 2 zusammen mit einer abnehmen- 
den Wellenlänge direkt vor 13:50 UTC (Abb. 6.9) deutet auf eine sich aufsteilende Welle hin, die 
nachfolgend bricht (wieder abnehmende Amplitude). Dies würde für eine Kelvin-Helmholtz- 
Instabilität sprechen. 
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Ist die Scherung der Windgeschwindigkeit groß und gleichzeitig die Brunt-Väisälä-Frequenz 
niedrig, können im stabilen Fall aufgrund der Geschwindigkeitsscherung Kelvin-Helmholtz- 
Instabilitäten auftreten und die internen Schwerewellen brechen. Ein Maß dafür, ab wann 
Kelvin-Helmholtz-Wellen auftreten können, ist die Richardson-Zahl: 

Ar2 



Ri 



dz 



(6.12) 



Für Ri < 1/4 (Houze, 1993) kann eine stabil geschichtete Scherströmung instabil werden, 
die Turbulenz wird verstärkt und die Kelvin-Helmholtz-Instabilität entsteht. Für Ri > 1/4 tritt 
dagegen Turbulenzabschwächung auf. Für den 08.01.2010 wird die Richardson-Zahl mit der 
oben beschriebenen, aus einer Radiosondenmessung berechneten, Brunt-Väisälä-Frequenz (Gl. 
6.7, Abb. 6.2) und mit dem vertikalen Gradienten des mittleren Horizontalwinds aus Messun- 
gen des Wolkenradars (Abb. 6.3) abgeleitet. Die berechnete Richardson-Zahl in Abhängigkeit 




10 10 
Richardson-Zahl 

Abbildung 6.10: Richardson-Zahlen auf Basis der Brunt-Väisälä-Frequenz für trockene Luft 
aus IVIessungen der in Stuttgart/Schnarrenberg (schwarz) und Idar-Oberstein (rot) gestarteten 
Radiosonden um 12:00 UTC und der Geschwindigkeiten des mittleren Horizontalwinds aus 
Messungen des Wolkenradars MIRA36-S um 14:18 UTC. 

von der Höhe ist in Abb. 6.10 dargestellt. Nach der Theorie liegt die wichtige Grenze der 
Richardson-Zahl für Turbulenzverstärkung, bzw. -abschwächung bei 1/4. Die Basis der Be- 
rechnung der Richardsonzahl bilden zum einen Profile der Temperatur und des Drucks ei- 
nes einzelnen Radiosondenaufstiegs. Diese Werte unterliegen allerdings großen Schwankun- 
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gen, sind abhängig von Messungenauigkeiten und fließen quadratisch in die Berechnung der 
Richardson-Zahl ein. Gleiches gilt für den Beitrag aus dem vertikalen Gradienten der mittleren 
Horizontalgeschwindigkeit. Wird alleine die Brunt-Väisälä-Frequenz minimal um Werte, die 
sich innerhalb ihrer Standardabweichung befinden, verändert, hat dies einen großen Einfluss 
auf die Werte der Richardson-Zahl. Ein Beispiel: Ist der Gradient der Horizontalgeschwindig- 
keit in einem Höhenabschnitt 2.1 m/s /lOO m und die Brunt-Väisälä-Frequenz im einen Fall 
10xlO~^s~^ und im anderen 15xl0~^s~^ , dann ist ist die Richardson-Zahl im ersten Fall 
0.23 (Turbulenzverstärkung) und im zweiten Fall bei 0.51 (Turbulenzabschwächung). Dies ist 
bei der Interpretation der Richardson-Zahl zu berücksichtigen. So bleibt anhand der Abb. 6.10 
festzustellen, dass Werte, die im Bereich der Grenze zwischen Turbulenzabschwächung und 
-Verstärkung liegen, in Höhen von 2.7 km bis 2.9 km, sowie in 3.2 km auftreten. D.h., dass 
im oberen Bereich der Welle 1, zwischen 2.7 km und 2.9 km eine Turbulenzverstärkung auf- 
tritt. Diese ist allerdings nur schwach ausgeprägt. Somit kann ein Aufsteilen der Welle 1 in 
Kombination mit Welle 2 und ein anschließendes Überschlagen der Welle und demzufolge eine 
Kelvin-Heknholtz-Instabilität nicht ausgeschlossen werden. 

Zusammenfassend ergeben sich für den 08.01.2010 für den Zeitraum von 13:30 UTC- 
14:17 UTC folgende Schlussfolgcrungcn mit den in diesem Kapitel gezeigten Verfahren zur 
Detektion von Wellen und der Abschätzung ihrer Wellenlänge und Phasengeschwindigkeit: 

• Eine ausgeprägte horizontal ausgerichtete Welle ist anhand von Periodizitäten der resi- 
dualen Geschwindigkeit in einer Höhe zwischen 2.5 km und 2.8 km zu erkennen. Die- 
se Welle beeinflusst darüberliegende Bereiche dahingehend, dass sich Periodizitäten in 
den residualen Geschwindigkeiten zeitweise unabhängig der horizontalen Windrichtung 

in näherungsweise die gleiche Richtung bewegen wie die Welle in 2.5 km bis 2.8 km. 
In Verbindung mit Richardson-Zahlen, die auf eine Turbulenzverstärkung in 2.7 km bis 
2.9 km schließen lassen, liegt die Vermutung nahe, dass es sich hierbei um Kelvin- 
Helmholtz -Wogen handelt. 

• Eine Welle kürzerer Wellenlänge (ca. 0.9 km) ist in einem Höhenabschnitt zwischen 
3.15 km und 3.35 km zeitweise zu erkennen. Während dieses Zeitraums ist der Einfluss 
der Welle in 2.5 km bis 2.8 km auf die darüber liegenden Schichten deutlich reduziert und 
die Richtung in die sich die Periodizitäten in den residualen Geschwindigkeiten bewegen 
entspricht annähernd der Windrichtung (relativ zur Scanrichtung). 
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Die Besonderheit des scannenden Wolkenradars MIRA36-S, in gewissen Grenzen den Azi- 
mut und die Elevation der Antenne zu variieren, wird genutzt, um u.a. dynamische Prozesse in 
Wolken zu detektieren und zu analysieren. Dafür werden RHI-Scans verwendet, die räumlich 
zweidimensionale Daten liefern. Durch eine schnelle Abfolge von RHI-Scans wird der Ge- 
winn an Informationen über besondere Strukturen in Wolken erhöht, beispielsweise sichtbar in 
der Reflektivität und deren Entwicklung. Mit den durchgeführten zweidimensionalen Messun- 
gen werden neuartige Mess- und Auswertungsstrategien entwickelt. Im Zuge dessen wird zum 
einen für Wolken, die Niederschlag produzieren, ein automatischer Algorithmus zur Detektion 
der Schmelzschicht entwickelt. Damit können beispielsweise Messbereiche, die von Eisparti- 
keln bestimmt werden, von rein flüssigem Niederschlag (Regen) separiert werden. Zum anderen 
werden Verfahren entwickelt, dynamische Prozesse in Wolken und deren Entwicklung zu un- 
tersuchen. 

Die Detektion der Schmelzschicht wird durch ein neuartiges, automatisiertes Verfahren 
ermöglicht. Für dieses Verfahren wird die Schmelzschichthöhe aus Profilen dreier Größen be- 
stimmt, die aus Messungen abgeleitet werden, der Reflektivität C, des linearen Depolarisations- 
verhältnisses L und der Vertikalgeschwindigkeit Wm- Eine Stärke des Verfahrens ist es, nicht 
nur die Höhe, sondern auch die Oberkante und die Unterkante der Schmelzschicht zuverlässig 
zu bestimmen. Dadurch wird die getrennte Betrachtung von drei Höhenabschnitten realisiert: 
Ein oberer Abschnitt, in dem Eiskristalle und Eispartikel vorliegen und in dem oft die Pro- 
duktion von Niederschlag beginnt, ein mittlerer Abschnitt, in dem Eispartikel schmelzen und 
ein unterer Abschnitt, in dem rein flüssige Hydrometeore (Tropfen) vorliegen. Die statistische 
Untersuchung dieser Bereiche führt zu folgenden Ergebnissen: 

• Die Median-Profile zeigen oberhalb der Schmelzschicht eine Zunahme der Partikelgröße 
von oben nach unten (nahezu linearer Anstieg der Reflektivität und der Sedimentations- 
geschwindigkeit der Streuer). 

• Die mittleren Partikelgrößen im Streuvolumen knapp oberhalb der Schmelzschicht, die 
für Schmelzschichtdicken von mehr als etwa 350 m verantwortlich sind, liegen bei D > 
3 mm. 

• Die Schmelzdauer eines Ensembles von Eispartikel ist für Schmelzschichtdicken zwi- 
schen 200 m und 500 m konstant und beträgt etwa 150 s, unabhängig vom Durchmesser 
der Eispartikel direkt oberhalb der Schmelzschicht. 

Weiterhin wurde gezeigt, dass alle drei Größen notwendig sind, die Schmelzschichtdicke bzw. 
die -Oberkante und -unterkante zuverlässig zu bestimmen. So konnte der Einfluss des Schmelz- 
vorgangs auf die drei Größen durch die außergewöhnlich hohe räumliche Auflösung gezeigt 
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werden. Zu Beginn des Schmelzprozesses reagieren die Werte der Reflektivität und der Sedi- 
mentationsgeschwindigkeit mit einem klar erkennbaren Anstieg. Erst wenn die Hydrometeore 
weitere 40 m gefallen sind, ist auch im Profil des LDR ein Anstieg der Werte zu erkennen. 
Zum Ende des Schmelzvorgangs hingegen ist zuerst im Profil der Reflektivität und nach weite- 
ren 50 m in den Profilen des LDR und der Sedimentationsgeschwindigkeit die Unterkante der 
Schmelzschicht erkennbar. 

Mit der Unterscheidung der gezeigten Höhenabschnitte ist die Basis für die Untersuchung dy- 
namischer Prozesse oberhalb der Schmelzschicht geschaffen. Dynamische Abläufe innerhalb 
einer Wolke, bzw. eines Mischbereichs aus Wolkenpartikeln und Niederschlagspartikeln, sind 
anhand der radialen Geschwindigkeiten nicht direkt erkennbar. Deshalb wurde eine neue Größe 
eingeführt, die residuale Geschwindigkeit. Mit ihrer Hilfe ist es möglich, Feinstrukturen der 
radialen Geschwindigkeit aufzulösen. Eine Funktion, die den mittleren Horizontalwind und 
die mitüere Vertikalgeschwindigkeit, sowie deren Trend in einer festen Höhe berücksichtigt, 
wird an die gemessenen radialen Geschwindigkeiten in einer festen Höhe mit der Methode der 
kleinsten Quadrate gefittet. Damit werden Profile der mittleren horizontalen Windgeschwin- 
digkeit und der mittleren Vertikalgeschwindigkeit abgeschätzt. Aus diesen Profilen kann das 
großräumige mittlere Windfeld berechnet und von den radialen Geschwindigkeiten abgezogen 
werden. Als Ergebnis erhält man die residualen Geschwindigkeiten. 

Die Erfahrung mit dieser neuen Größe hat gezeigt, dass die residualen Geschwindigkeiten 
meist vertikale Umlagerungen der Stieuer anzeigen. Dadurch können erstmals gemessene 
kleinräumige Fluktuation der radialen Geschwindigkeiten räumlich hoch aufgelöst in kurzem 
zeitlichen Abstand von 30 s gezeigt werden. Diese Fluktuationen können meist entweder Auf-, 
und Abwindbereichen oder verschiedenen Größen von Stteuem zugeordnet werden. Für die- 
se Unterscheidung ist die kombinierte Betrachtung der residualen Geschwindigkeiten mit der 
mittleren Vertikalgeschwindigkeit der Streuer und der Reflektivität entscheidend. 

Anhand einzelner Messtage konnte gezeigt werden, dass die residualen Geschwindigkeiten 
wertvolle zusätzliche Informationen über die Dynamik in Wolken liefern. Beispielsweise lässt 
sich anhand der residualen Geschwindigkeiten die Entwicklung eines Bereichs erhöhter Re- 
flektivität in einer Wolke abschätzen. Wird ein Aufwind, Konvektion, vorderseitig des Be- 
reichs erhöhter Reflektivitäten in Bewegungsrichtung detektiert, ist in manchen Fällen mit ei- 
ner Verstärkung von einem RHI-Scan zum nächsten zu rechnen. Wird ein starker Abwind vor- 
derseitig detektiert, ist in manchen Fällen mit einer Abschwächung des Reflektivitätsbereichs 
zu rechnen. Ein Aufwindbereich kann einerseits detektiert werden, wenn sich die gemessenen 
Stieuer nach oben bewegen. In diesem Fall ist eine positive residuale Geschwindigkeiten zu 
erkennen, die die mittlere Sedimentationsgeschwindigkeit in der zugehörigen Höhe übertrifft. 
Andererseits ist eine Detektion auch möglich, wenn die Stieuer im Fallen durch den Aufwind 
abgebremst werden. Dies ist deutlich schwieriger, da aus einer Kombination von Messgrößen 
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abgewogen werden muss, ob es sich um einen Aufwind handelt. Sind Bereiche positiver Resi- 
dualgeschwindigkeit kleiner als die Sedimentationsgeschwindigkeit können Aufwindbereiche 
detektiert werden, wenn die Reflektivitäten im Bereich positiver residualer Geschwindigkeiten 
denen links und rechts davon entsprechen. Dies ist der Fall, wenn die Aufwindgeschwindigkeit 
nicht ausreicht, der Sedimentationsgeschwindigkeit der Streuer entgegen zu wirken und die 
Streuer dadurch lediglich im Fallen abbremst, dabei aber Luft nach oben transportiert wird. In 
diesem Zusammenhang wurde gezeigt, dass die Erhöhung der Reflektivität in Zusammenhang 
mit größeren Partikeln steht, die verstärkt ausfallen und zu einer größeren Schmelzschichtdicke 
und einer Erhöhung des Niederschlags beitragen. 

Neben der Analyse räumlich begrenzter Strukturen erhöhter Reflektivität, die durch das Mess- 
gebiet transportiert werden, werden mit Hilfe der residualen Geschwindigkeiten auch Unter- 
suchungen von horizontal ausgerichteten internen Schwerewellen ermöglicht. Ist der Bereich, 
in dem WeUen vermutet werden, bekannt, beispielsweise anhand einer periodischen Abfolge 
wechselnder Vorzeichen der residualen Geschwindigkeit in einer festen Höhe, wird die Wel- 
lenlänge in Scanrichtung aus den Fouriertransformierten residualen Geschwindigkeiten in einer 
festen Höhe berechnet. Mit Hilfe der Kreuzkorrelation der residualen Geschwindigkeiten in 
einer festen Höhe wird die Verlagerungsgeschwindigkeit der Welle in Scanrichtung ermittelt. 
Der Vergleich der Verlagerungsgeschwindigkeit mit der horizontalen Windgeschwindigkeit in 
Scanrichtung liefert die Ausbreitungsgeschwindigkeit der Welle in Scanrichtung. Anhand des 
Vorzeichens kann unterschieden werden ob sich die Welle in bestimmten Höhenbereichen in 
Richtung oder entgegen der Richtung des horizontalen Windes in Scanrichtung ausbreitet. 

Ein Ergebnis dieser Untersuchung ist, dass Wellen unterhalb der durch Scherung induzierten In- 
versionsschicht zeitweise Wellen oberhalb auslösen können. Die Wellen oberhalb der Inversion 
breiten sich dadurch zeitweise entgegen der Windrichtung aus. Durch die schnelle Abfolge von 
RHI-Scans kann der zeitliche Verlauf der Wellenlängenänderung sowie der damit verbundenen 
Amplitude der residualen Geschwindigkeiten verfolgt werden. Dadurch wurde ein Zusammen- 
hang des Verhaltens der Wellenlängen der Welle oberhalb mit jener unterhalb der Inversion 
festgestellt. Das heißt, dass zeitweise die untere Welle den Bereich oberhalb der Inversion be- 
einflusst. Ob es sich in einem solchen Fall um Kelvin-Hehnholtz-Instabilitäten handelt, wird mit 
Hilfe der Richardson-Zahl überprüft, die ein Maß für die Turbulenzzunahme, oder -abnähme in 
der Atmosphäre ist. Das Ergebnis zeigt in diesem Zusammenhang Werte der Richardson-Zahl in 
Höhen, in denen Wellen auftreten, die ein Vorhandensein von Kelvin-Helmholtz-Instabilitäten 
nahelegen. 

Die Untersuchung sowohl dynamischer als auch physikalische Prozesse in Wolken und deren 
Entwicklung werden Dank der gezeigten Verfahren mit einem einzigen Scanmuster ermöglicht: 
mit schnell durchgeführten RHI-Scans. Dies erlaubt die direkte Bestimmung der Schmelz- 
schichtgrenzen, sowie dynamischer Verlagerungen von kleinräumigen Wolkenbereichen im 
operationeilen Betrieb. Zudem lässt sich die räumliche Entwicklung dynamischer Strukturen 
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in einem kurzen zeitlichen Abstand von 30 s ableiten. Die Veränderang von lokalen konvekti- 
ven Bereichen oder horizontalen Wellen in Wolken können dann mit Hilfe der in dieser Arbeit 
definierten residualen Geschwindigkeiten analysiert werden. 

Die residuale Geschwindigkeit ist eine neue Größe, die enormes Potential hat, dynamische Pro- 
zesse in Wolken zu verstehen. Ein interessanter Aspekt für zukünftige Untersuchungen ist, wie 
sich Muster in den residualen Geschwindigkeiten relativ zum Horizontalwind und der Sedimen- 
tationsgeschwindigkeit der Streuer fortbewegen. Dies könnte mit Hilfe angepasster Tracking- 
Algorithmen (COTREC, Li et al., 1995) untersucht werden. Im Zuge dessen sind die residualen 
Geschwindigkeiten auch auf Doppler-Messungen weiterer Messinstrumente, beispielsweise ei- 
nes Lidars, anwendbar. Damit würden zusätzliche parallele Messungen (schnelle RHI-Scans) 
eines Wolkenradars mit einem Lidar die Visualisierung und Entwicklung von Auf-, Abwind- 
bereichen in der Troposphäre vervollständigen und Aufschlüsse geben, beispielsweise über die 
Einmischung von Umgebungsluft in Wolken (Entrainment) oder den konvektiven Transport von 
feuchtwarmen Luftmassen in die Höhe, die zur Entwicklung von Cumulus-Wolken führen. 



107 



A Minimierung der Kostenfunktion 



Die Minimierung der Kostenfunktion (Gl. 5.8) in Kap. 5.1 erfolgt nach dem Prinzip der klein- 
sten Quadrate. Die detaillierte Vorgehensweise wird im Folgenden beschrieben. 



Einsetzten von Gl. 5.6 in Gl. 5.8 liefert: 

n 

J(z) = J2i<^i^ z)-iafi + bh + ch)f 



(1.1) 



i=l 



mit den Substitutionen 



/i = s\u{ez{xi,z)) 

dw 



/2 = COs{€z{Xi,z)) 



1 



d 



u 



b ^ w^{z) - Q^{z)z 



cos{ez{xi,z) ^^2) 

du 



Zur Berechnung der unbekannten Größen a, b, c wird das Minimum der Kostenfunktion 
abhängig von den gesuchten Größen berechnet. Dazu werden die Ableitungen dJ/da, dJ/dh 
und dJ/db gleich Null gesetzt: 



dJ 
da 

dJ 

~db 

dJ_ 
de 



2 Y.{v{xi, z)-{ah + b /2 + c/3))(-/i) = 0 

i=l 
n 

2 z)-iaf, + b h + c/3))(-/2) = 0 

i=l 
n 

2 Y,{v{xi, z)-{ah + b h + cf3))i-h) = 0 



(1.3) 



i=l 



Daraus folgt nach Umformungen 
/ 



i=l /l Z^i=l /l/2 Z^i=l /l/3 

En ff sr^n f2 V^n f f 

i=i /2/1 2^i=i h z^i=i J'iJ^ 

En ff sr^n f r v-vn f2 

^ i=l/3/l 2^i=l/3/2 2^i=i/3 y 

und entspricht in Vektorschreibweise 



a 



\E7=lV{Xi,z)f3j 



(1.4) 



Ax = b 



Für X = (a, b, c) ergibt sich 



(1.5) 



(1.6) 



falls A invertierbar ist. Für det A ^ 0 ist die quadratische (3,3)-Matrix invertierbar und die 
Lösung eindeutig bestimmbar: 



Xi 



detAi 
detA 



(1.7) 



Für Aj wird die z-te Spalte von A durch den Vektor b ersetzt. 
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B Symbolverzeichnis 



a 


[m] 


Kugelradius 




[m/s] 


Koeffizienten der Fourierreihenentwicklung von Vj.^s{oi) 


A 


[m/s] 


Amplitude der residualen Geschwindigkeit interner Schwerewellen 


b 


[m] 


Halbe Dicke der Schmelzschicht 


K 


[m/s] 


Koeffizienten der Fourierreihenentwicklung von Vj.^s{oi) 


B 


[m] 


Dicke der Schmelzschicht 


Br 


[Hz] 


Empfängerbandbreite 


c 


[m/s] 


Vakuum-Lichtgeschwindigkeit 


Cp 


[J/(kgK)] 


spezifische Wärmekapazität bei konstantem Druck 


Cw 


[m/s] 


horizontale Phasengeschwindigkeit einer internen Schwerewelle 


C 


[W/m] 


Radarkonstante 


D 


[m] 


Durchmesser 


er 




Einheitsvektor in radialer Richtung 


f 


[Hz] 


Frequenz des gesendeten Radarsignals 


fü 


[Hz] 


Doppler-Frequenz 


fm 


[Hz] 


Diskrete Frequenzen des Leistungsspektrums Pk 


fp 


[Hz] 


Pulswiederholfrequenz 


F 




Funktion des mexican hat 


Fr 


[dB] 


Rauschzahl des Empfängers 


G 




Antennengewinn 


J 


[(m/s)2] 


Kostenfunktion 


k 




Einheitsvektor in 2:-Richtung 


k 


[m-i] 


Wellenzahl 


kß 


[J/K] 


Boltzmann-Konstante 


< ? 




Dielektrizitätsfaktor von Wasser oder Eis 


L 


[dB] 


Lineares Depolarisationsverhältnis (LDR) 
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Ly [dB] Leistungsverluste 

m komplexer Brechungsindex 

M [g/m''] Flüssigwassergehalt 

Mfc Maximale Anzahl der mit Daten besetzten Gitterpunkte in einer fe- 

sten Höhe 

n Nyquistzahl (gerundete Faltungszahl) 

rij Index in x-Richtung, zentriert im Zenit eines RHI-Scans 

h Faltung s zahl 

npFT FFT-Länge (Anzahl der Pulszyklen, die in die FFT einfließen) 

ns Anzahl der Spektren (FFTs), über die gemittelt wird 

N [s~^] Brunt-Väisälä-Frequenz für trockene Luft 

Np (D) Größenverteilung der Streupartikel in einem Mess volumenausschnitt 

p [hPa] Luftdruck 

Pk [V/s] Leistungsspektrum 

Pr [W] mittlere empfangene Leistung 

Pt [W] Mittlere Sendeleistung 

Pn [W] Rauschleistung des Empfängers 

R [m] radiale Entfernung zum Radar 

Ri Richardson-Zahl 

Rl [J/(kg K)] spezifische Gaskonstante für trockene Luft 

s-r [m] Pulslänge 
S Schmelzschichtindex 

S* [m/s^] spektrale Signalstärke 

Sk [V/m] komplexes Eingangssignal 

S,^ Schmelzschichtindex für ( 

Sl Schmelzschichtindex für L 

Sw Schmelzschichtindex für 

SNR Signal-Rausch-Verhältnis 



III 



t 


[s] 


Zeit 




[s] 


Schmelzdauer eines Ensembles von fallenden Eisteilchen 




[s] 


Mittelungszeit über npFT x ris Pulse 


tp 


[s] 


Pulswiederholdauer 


T 


[K] 


Temperatur 


I 


[°C] 


Feuchttemperatur 


To 


[K] 


Referenztemperatur zur Bestimmung der Rauschzahl 


Tn 


[K] 


Rauschtemperatur des Empfängers 


u 


[m/s] 


Horizontale Komponente der Geschwindigkeit der Streuer in x- 



Ub 



Ur. 



Uc,R 



U 



m/s 
m/s 

m/s 

m/s 

m/s 



Richtung 

Horizontaler Hintergrundwind in Richtung des RHI-Scans 

aus Vr abgeleitete (RHI-Scans) mittlere horizontale Komponente der 
Geschwindigkeit der Streuer 

Kombinierter horizontaler Hintergrundwind aus Radiosonden- und 
C-Band-Radar-Messungen 

lokale Horizontalgeschwindigkeit 

Horizontale Komponente der Geschwindigkeit der Streuer in y- 
Richtung 



Vb 


[m/s] 


Hintergrund-Windfeld 


VN 


[m/s] 


Nyquist-Geschwindigkeit 


Vr 


[m/s] 


entfaltete Radialkomponente der Geschwindigkeit der Streuer 


Vr,S 


[m/s] 


gemessene Radialkomponente der Geschwindigkeit der Streuer 


VRes 


[m/s] 


residuale Geschwindigkeit 


VRes 


[m/s] 


Fouriertransformierte der residualen Geschwindigkeit 


Vs 


[m/s] 


Dreidimensionale Geschwindigkeit der Streuer 


Vw 


[m/s] 


Dreidimensionale Windgeschwindigkeit 


V* 


[m/s] 


Radialkomponente des großskaligen Windfeldes 


w 


[m/s] 


Vertikale Komponente der Geschwindigkeit der Streuer 


Wb 


[m/s] 


Vertikale Hintergrund-Streuergeschwindigkeit 
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[m/s] aus abgeleitete (RHI-Scans) mittlere vertikale Komponente der 
Geschwindigkeit der Streuer 

wp [m/s] Sedimentationsgeschwindigkeit relativ zur Vertikalbewegung der 
Luft 

w* Hamming-Fenster 

w [m/s] lokale Vertikalgeschwindigkeit 

Zh [m] Untergrenze der Schmelzschicht 

Zp [m] Höhe der Schmelzschicht im Maximum von S 

Zt [m] Obergrenze der Schmelzschicht 

Z [mm^/m^] Radarreflektivitätsfaktor 

a [°] Azimut der Radarantenne 

ctc,R [°] Kombinierte horizontale Windrichtung aus Radiosonden- und C- 

Band-Radar-Messungen 

aR Radioelektrische Größe 

aw [°] Ausbreitungsrichtung einer internen Schwerewelle 

ß Korrekturfaktor für die Trägheit von Streukörpem 

AV [m^] Volumenauschnitt des Radarstrahls (Range Gate) 

e [°] Elevation der Radarantenne 

ez [°] Zenitwinkel der Radarantenne 

C [dBz] logarithmierter Radarreflektivitätsfaktor 

T] [1/m] Radarreflektivität 

9 [K] Potentielle Temperatur trockener Luft 

9a [°] Antennenöffnungswinkel 

A [m] Radarwellenlänge 

Xw [m] Wellenlänge einer internen Schwerewelle 

ß [g/m^] Dichte 

qf,s Korrelationskoeffizient der Korrelation zwischen F und S 

CT [m^] Rückstreuquerschnitt eines Einzelstreuers 

(Ta [m^] Absorptionsquerschnitt 
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(7y [m/s] Standardabweichung der Streuergeschwindigkeiten 

iPr Phase zweier aufeinanderfolgender Pulse 

r [s] Pulsdauer 

td [s] Dekorrelationszeit 

[°/s] azimutale Scangeschwindigkeit 

[°/s] zenitale Scangeschwindigkeit 



114 B Symbolverzeichnis 



Literatur 



115 



Literatur 

Ackerman, T. P. und G. Stokes, 2003: The Atmospheric Radiation Measurement Program. Phys. 
Today, 56, 38-45. 

Austin, P. M. und A. C. Bemis, 1950: A Quantitative Study of the "Bright Band" in Radar 
Precipitation Echoes. /. Atmos. Sei., 7, 145-151. 

Battan, L., (Ed.) , 1973: Radar Observation of the Atmosphere. The University of Chicago 
Press. 

Beheng, K. D., 1998: Grundlagen der Wolkenmikrophysik und der Dynamik von Wolken und 
Fronten. Herbstschule Radarmeteorologie, Ann. Met. 38, Deutscher Wetterdienst, 7-24. 

Blahak, U., 2005: Analyse des Extinktionseffektes bei Niederschlagsmessungen mit einem C- 
Band-Radar anhand von Simulationen und Messungen. Dissertation, Universität Karlsruhe. 

Bringi, V. und V. Chandrasekar, 2005: Polarimetrie Doppler Weather Radar. Cambridge Uni- 
versity Press. 

Bronstein, I., K. Semendjajew, G. Musiol, und H. Mühlig, 2008: Tasehenhandbueh der Mathe- 
matik. Verlag Harri Deutsch. 

Browning, K. A. und R. Wexler, 1968: A Determination of Kinematic Properties of a Wind 
Field Using Doppler Radar. J. Appl. Meteor., 7, 105 - 113. 

Carbone, R. E., J. D. Tuttle, D. A. Ahijevych, und S. B. Trier, 2002: Inferences of Predictability 
Associated with Warm Season Precipitation Episodes. /. Atmos. Sei., 59, 2033 - 2056. 

Clothiaux, E. E., M. A. Miller, B. A. Albrecht, T. P Ackerman, J. Verlinde, D. M. Babb, R. M. 
Peters, und W. J. Syrett, 1995: An evaluation of a 94-GHz radar for remote sensing of cloud 
properties. /. Atmos. Oeeanie TeehnoL, 12, 201-228. 

Damiani, R., G. Vali, und S. Haimov, 2006: The structure of thermals in cumulus from airborne 
dual-Doppler radar observations . J. Atmos. Sei., 63, 1432 - 1450. 

Dotzek, N., 1999: Mesoskalige numerische Simulation von Wolken- und Niederschlagsprozes- 
sen über strukturiertem Gelände. Dissertation, Universität Karlsruhe. 

Doviac, R. J. und D. S. Zmic, 1984: Doppler Radar and Weather Observation. Academic Press. 



116 



Literatur 



Evans, A. G., J. D. Locatelli, M. T. Stoelinga, und P. V. Hobbs, 2005: The IMPROVE-1 storm 
of 1-2 February 2001. Part II: Cloud structures and the formation of precipitation. /. Arnos. 
Sei., 62, 3456 - 3473. 

Fabry, F. und I. Zawadzki, 1995: Long-Term Radar Observations of the Melting Layer of Pre- 
cipitation and Their Interpretation. J. Atmos. Sei., 52, 838-851. 

Frisch, A., C. Fairall, und J. Gibson, 1995: Doppler radar measurements of turbulence in marine 
stratiform cloud during ASTEX. /. Atmos. Sei., 52, 2800 - 2808. 

Geerts, B., R. Damiani, und S. Haimov, 2006: Fine-scale vertical structure of a cold front as 
revealed by airborne radar. Mon. Weather Rev., 134, 251 - 272. 

Giangrande, S. E., J. M. Krause, und A. V. Ryzhkov, 2008: Automatic Designation of the Mel- 
ting Layer with a Polarimetrie Prototype of the WSR-88D Radar. J. Appl. Meteor. ClimatoL, 
47, 1354- 1364. 

Grenzhäuser, J. und J. Handwerker, 2010: Residual velocities and wind profiles in clouds as 
measured by a scanning 35.5GHz cloud radar. Proceedings ofthe Sixth European Conference 
on Radar in Meteorology and Hydrology. 6-10 September 2010, Sibiu, Romania. 

Gunn, K. L. und T. Fast, 1954: The microwave properties of precipitation particles. Q. J. R. 
Meteorol. Soe., 80, 522 - 545. 

Gysi, H., R. Hannesen, und K. D. Beheng, 1997: A method for bright-band correction in hori- 
zontal rain intensity distributions. Proceedings ofthe 28th Radar Conference. 7-12 September 
1997, AMS, Austin, USA, 214 - 215. 

Hallett, J. und S. C. Mossop, 1974: Production of secondary ice particles during the riming 
process. Natura, 249, 26-28. 

Handwerker, J. und U. Görsdorf, 2006: A Comparison of Vertical Cloud Radar Profiles with 
RHI Scans. Atmospheric Profiles in Research and Operations, COST 720. 15-18 May 2006, 
Toulouse, France. 

Handwerker, J. und A. Miller, 2008: Intercomparison of measurements obtained by vertically 
pointing coUocated 95 GHz and 35.5 GHz radars. Proceedings ofthe Fifth European Confe- 
rence on Radar in Meteorology and Hydrology. 26-30 May, 2008, Helsinki, Finland. 

Handwerker, J., K. Träumner, J. Grenzhäuser, und A. Wieser, 2008: Simultaneous wind mea- 
surements with lidar and cloud radar - complementarity and quality check. Proceedings of 
the Fifth European Conference on Radar in Meteorology and Hydrology. 26-30 May 2008, 
Helsinki, Finland. 

Hannesen, R., 1998: Analyse konvektiver Niederschlagssysteme mit einem C-Band Dopplerra- 
dar in orographisch gegliedertem Gelände. Dissertation, Universität Karlsruhe. 



Literatur 



117 



Hennington, L., 1981: Reducing the effects of Doppler radar ambiguities. J. Appl. Meteor, 20, 
1543-1546. 

Hogan, R. J. und E. J. O'Connor, 2004: Facilitating cloud radar and lidar algorithms: The cloud- 
net Instrument synergy/target categorization product. cloudnet documentation. Available on- 
line at www.cloud-net.org/data/products/categorize.html. 

Houze, R. A., 1993: Cloud Dynamics. Academic Press. 

Illingworth, A. J., R. J. Hogan, E. J. O'Connor, D. Bouniol, J. Delanoe, J. Pelon, A. Protat, 
M. E. Brooks, N. Gaussiat, D. R. Wilson, D. R Donovan, H. K. Baltink, G.-J. van Zadelhoff, 
J. D. Eastment, J. W. F. Goddard, C. L. Wrench, M. Haeffelin, O. A. Krasnov, H. W. J. Rus- 
schenberg, J.-M. Piriou, F. Vinit, A. Seifert, A. M. Tompkins, und U. WiUen, 2007: Cloudnet 
- Continuous evaluation of cloud profiles in seven operational models using ground-based 
observations. Bull. Amer. Meteor Soc, 88, 883 - 898. 

KoUias, P. und B. Albrecht, 2000: The turbulent structure in a continental stratocumulus cloud 
from millimeter wavelength radar observations. /. Atmos. Sei., 57, 2417 - 2434. 

KoUias, F., B. Albrecht, R. Lhermitte, und A. Savtchenko, 2001: Radar observations of updrafts, 
downdrafts, and turbulence in fair weather cumuü. J. Atmos. Sei., 58, 1750 - 1766. 

Kottmeier, C, N. Kalthoff, C. Barthlott, U. Corsmeier, J. van Baelen, A. Behrendt, R. Beh- 
rendt, A. Blyth, R. Coulter, S. Creweil, P. di Girolamo, M. Dominger, C. Flamant, T. Foken, 
M. Hagen, C. Hauck, H. Höller, H. Konow, M. Kunz, H. Mahlke, S. Mobbs, E. Richard, 
R. Steinacker, T. Weckwerth, A. Wieser, und V. Wulfmeyer, 2008: Mechanisms initiating 
deep convection over complex terrain during COPS . Meteorol. Z., 17, 931-948. 

Leuenberger, D. und A. Rossa, 2007: Revisiting the latent heat nudging scheme for the rainfall 
assimilation of a simulated convective storm. Meteorol. Atmos. Phys., 98, 195-215. 

Li, L., W. Schmid, und J. Joss, 1995: Nowcasting of motion and growth of precipitation with 
radar over a complex orography. /. Appl. Meteor, 34, 1286 - 1300. 

Matrosov, S. M., K. A. Clark, und D. E. Kingsmill, 2007: A Polarimetrie Radar Approach to 
Identify Rain, Melting-Layer, and Snow Regions for Applying Corrections to Vertical Profiles 
of Reflectivity. /. Appl. Meteor. ClimatoL, 46, 154 - 166. 

Melchionna, S., M. Bauer, und G. Peters, 2008: A new algorithm for the extraction of cloud Pa- 
rameters using multipeak analysis of cloud radar data. First application and results. Meteorol. 
Z, 17, 613-620. 

Meneghini, R. und T. Kozu, (Eds.) , 1990: Spaceborne weather radar. Artech House. 



118 



Literatur 



Park, H., A. V. Ryzhkov, D. S. Zmic, und K. Kim, 2008: The Hydrometeor Classification Algo- 
rithm for the Polarimetrie WSR-88D: Description and Application to an MCS. Weather and 
Forecasting, 24, 730 - 748. 

Pazmany, A. L., R. E. Mcintosh, R. Kelly, und G. Vali, 1994: An airborne 95-GHz dual- 
polarized radar for cloud studies. IEEE Trans. Geosci. Remote Sens., 32, 731-739. 

Peters, T., 2008: Ableitung einer Beziehung zwischen der Radarreflektivität, der Niederschlags- 
rate und weiteren aus Radardaten abgeleiteten Parametern unter Verwendung von Methoden 
der multivariaten Statistik. Dissertation, Universität Karlsruhe. 

Pfeifer, M., G. C. Craig, M. Hagen, und C. Keil, 2008: A Polarimetrie Radar Forward Operator 
for Model Evaluation. /. Appl Meteor. Climatol, 47, 3202 - 3220. 

Plank, V., D. Atlas, und W. Paulsen, 1955: The nature and detectability of clouds and precipita- 
tion as determined by 1.25 centimeter radar. /. MeteoroL, 1, 358 - 378. 

Probert- Jones, J. R., 1962: The radar equation in meteorology. Q. J. R. MeteoroL Soc, 88, 485 
-495. 

Pruppacher, H. R. und J. D. Klett, 1980: Microphysics of clouds and precipitation. D. Reidel 
Publishing Company. 

Rinehart, R. E., 1981: A pattern recognition technique for use with conventional weather radar 
to determine internal storm motions. Recent progress in radar meteorology. Atmos. TechnoL, 
119-134. 

Sassen, K., J. R. Campbell, J. Zhu, P. KoUias, M. Shupe, und C. Williams, 2005: Lidar and 
Triple-Wavelength Doppler Radar Measurements of the Melting Layer: A Revised Model for 
Dark- and Brightband Phenomena. /. Appl. Meteor., 44, 301 - 312. 

Sauvageot, H., 1992: Radar Meteorology. Artech House. 

Skolnik, M., 1990: Radar Handbook. McGraw-Hill. 

Stephens, G. L. und P. J. Webster, 1981: Clouds and climate: sensitivity to simple Systems. /. 
Atmos. Sei., 38, 235 - 247. 

Stewart, R. E., J. D. Marwitz, J. C. Pace, und R. E. Carbone, 1984: Characteristics through the 
melting layer of stratiform clouds. /. Atmos. Sei., 41, 3227 - 3237. 

Stokes, G. M. und S. E. Schwartz, 1994: The Atmospheric Radiation Measurement (ARM) 
Program: Programmatic background and design of the Cloud and Radiation Test Bed. Bull. 
Amer. Meteor Soc., 75, 1201-1221. 

Tabary, R, G. Scialom, und U. Germann, 2001: Real-Time Retrieval of the Wind from Aliased 
Velocities Measured by Doppler Radars. J. Atmos. Oceanic TechnoL, 18, 875-882. 



Literatur 



119 



Träumner, K., J. Handwerker, A. Wieser, und J. Grenzhäuser, 2010: A Synergy Approach to 
Estimate Properties of Raindrop Size Distributions Using a Doppler Lidar and Cloud Radar. 
J. Arnos. Oceanic Technol, 27, 1095 - 1100. 

Vivekanandan, J., D. S. Zrnic, S. M. Ellis, R. Oye, A. V. Ryzhkov, und J. M. Straka, 1999: 
Cloud microphysics retrieval using S-band dual-polarization radar measurements. Bull. Amer. 
Meteor. Soc, 80, 381-388. 

Waldteufel, R und H. Corbin, 1979: On the Analysis of Single Doppler Data. /. Appl. Meteor, 
18, 532-542. 

White, A. B., D. J. Cottas, E. T. Strem, F. M. Ralph, und R J. Neiman, 2002: An Automated 
Brightband Height Detection Algorithm for Use with Doppler Radar Spectral Moments. /. 
Atmos. Oceanic Technol., 19, 687 - 697. 

Woods, C. R, M. T. Stoelinga, und J. D. Locatelli, 2008: Size Spectra of Snow Particles Mea- 
sured in Wintertime Precipitation in the Pacific Northwest. /. Atmos. Sei., 65, 189 - 205. 

Wulfmeyer, V., A. Behrendt, H.-S. Bauer, C. Kottmeier, U. Corsmeier, A. Blyth, C. Craig, 
U. Schumann, M. Hagen, S. Creweil, F. Di Girolamo, C. Flamant, M. Miller, A. Montani, 
S. Mobbs, E. Richard, M. W. Rotach, M. Arpagaus, H. Russchenberg, F. Schlüssel, M. König, 
V. Gärtner, R. Steinacker, M. Dorninger, D. D. Turner, T. Weckwerth, A. Hense, und C. Sim- 
mer, 2008: Research campaign: The Convective and Orographically Induced Precipitation 
Study. Bull. Amer. Meteor Soc, 89, 1477-1486. 



Copyrighted maerial 



Danksagung 



An dieser Stelle möchte ich all jenen Danken, die mich während der Zeit der Promotion un- 
terstützt haben und so zum Gelingen dieser Arbeit beigetragen haben. Mein Dank gilt in erster 
Linie Prof. Dr. Klaus D. Beheng, der mir diese Arbeit ermöglichte und mir damit die Gele- 
genheit gab, in die spannende Welt der Radarmeteorologie einzutauchen. Er unterstützte mich 
stets mit kreativen Anregungen und war in Gesprächen immer wieder inspirierend und zugleich 
motivierend. 

Prof. Dr. Sarah Jones danke ich ganz herzlich für die Übernahme des Korreferats und das Inter- 
esse an dieser Arbeit. 

Des Weiteren gilt mein Dank Dr. Jan Handwerker, einerseits für die wertvollen Ratschläge zur 
Arbeit, zur Programmierung (und zu den grammatikalischen Schwächen des hiesigen Dialekts), 
andererseits für ausführliche Diskussionen über wissenschaftliche Themen und Privates, auch 
wenn wir nicht immer einer Meinungen waren. Außerdem möchte ich auch ausdrücklich Dr. 
Ulrich Blahak danken, der mir immer wieder mit seiner positiven Art, an Herausforderungen 
heranzugehen, mit gutem Rat zur Seite stand. 

Der gesamten Arbeitsgruppe danke ich für eine angenehme Arbeitsatmosphäre und die gemein- 
same Zeit auch außerhalb der Arbeit, beispielsweise beim Ski- bzw. Snowboardfahren. 

Weiterer Dank gebührt Gabi Klinck, die immer zur Stelle war, wenn es Rechnerprobleme gab, 
sowie Dr. Heike Noppel, Dr. Tim Peters, Dr. Winfried Straub, Isabelle Wolff und Malte Neuper, 
die mich durch Einblicke in ihre Arbeit und Diskussionen zur Wolkenphysik und Radarme- 
teorologie bereicherten und auch außerhalb der Arbeit den ein oder anderen Abend zu einem 
Erlebnis werden ließen. 

Im privaten Umfeld möchte ich mich für die jahrelange mentale Unterstützung während des 
Studiums und der Dissertation bedanken bei meiner Familie, ausdrücklich bei Heidi Keller und 
Gertrud Lenhard und meinen Freunden, insbesondere bei der Familie Beyle, Andrea Grohmer, 
Robert Heinke, Dr. Romi Sasse, Peter Hiller sowie Dagmar Fichtner und bei all jenen, die an 
dieser Stelle nicht explizit erwähnt sind. 

Ganz besonderer Dank gebührt meiner Frau Viviane Herrmann-Grenzhäuser. Sie hat mir in die- 
ser Zeit nicht nur eine wundervolle Tochter, Johanna Sophia Grenzhäuser, geschenkt, sondern 
stand in jeder Situation hinter mir und war und ist eine bedeutende Stütze für mich, auch wenn 
ich mal den einen oder anderen Abend der Arbeit widmete. 



Copyrighted maerial 



Wissenschaftliche Berichte des 

Instituts für Meteorologie und Klimaforschung 

des Karlsruher Instituts für Technologie 

(0179-5619) 



Bisher erschienen: 



Nr. 1: Fiedler, F. I Prenosil, T. 

Das MESOKLIP-Experiment. (Mesosicaliges Klimaprogramm im 
Oberrheintal). 

August 1980 

Nr. 2: Tangermann-Dlugi, G. 

Numerische Simulationen atmosphärischer Grenzschichtströmun- 
gen über langgestreckten mesoskaligen Hügelketten bei neutraler 
thermischer Schichtung. 

August 1982 

Nr. 3: Witte, N. 

Ein numerisches Modell des Wärmehaushalts fließender Gewässer 
unter Berücksichtigung thermischer Eingriffe. 

Dezember 1982 



Nr. 4: Fiedler, F I Höschele, K. (Hrsg.) 

Prof. Dr. Max Diem zum 70. Geburtstag. 
Februar 1 983 (vergriffen) 



Nr. 5: Adrian, 6. 

Ein Initialisierungsverfahren für numerische mesoskalige Strömungs- 
modelle. 

Juli 1985 



Nr. 6: Dorwarth, G. 

Numerische Berechnung des Druckwiderstandes typischer Gelände- 
formen. 

Januar 1986 



Nr. 7: Vogel, B.; Adrian, G. I Fiedler, F 

MESOKLIP-Analysen der meteorologischen Beobachtungen von 
mesoskaligen Phänomenen im Oberrheingraben. 

November 1987 



Nr. 8: Hugelmann, C-P. 

Differenzenverfahren zur Behandlung der Advektion. 

Februar 1988 



Copyrighted maienal 



Nr. 9: Hafner, T. 

Experimentelle Untersuchung zum Druckwiderstand der Alpen. 

April 1988 

Nr. 10: Corsmeier, U. 

Analyse turbulenter Bewegungsvorgänge in der maritimen 
atmosphärischen Grenzschicht. 

Mai 1988 

Nr. 1 1 : Walk, O. I Wieringa, J. (eds) 

Tsumeb Studies of the Tropical Boundary-Layer Climate. 

Juli 1988 

Nr. 12: Degrazia, G. A. 

Anwendung von Ähnlichkeitsverfahren auf die turbulente Diffusion 
in der konvektiven und stabilen Grenzschicht. 

Januar 1989 

Nr. 13: Schädler, G. 

Numerische Simulationen zur Wechselwirkung zwischen Landober- 
flächen und atmophärischer Grenzschicht. 
November 1990 

Nr. 14: Heidt, K. 

Untersuchungen zur Überströmung eines mikroskaligen Hindernisses 
in der Atmosphäre. 

Juli 1991 

Nr. 15: Vogel, H. 

Verteilungen reaktiver Luftbeimengungen im Lee einer Stadt - 
Numerische Untersuchungen der relevanten Prozesse. 

Juli 1991 

Nr. 16: Höschele, K.(ed.) 

Flanning Applications of Urban and Building Climatology - Proceedings 
of the IFHP/CIB-Symposium Berlin, October 14-15, 1991. 
März 1992 

Nr. 17: Frank, H. P. 

Grenzschichtstruktur in Fronten. 

März 1992 

Nr. 18: Müller, A. 

Parallelisierung numerischer Verfahren zur Beschreibung von 
Ausbreitungs- und chemischen Umwandlungsprozessen in der 
atmosphärischen Grenzschicht. 

Februar 1996 



Copyrlghted maienal 



Nr. 19: Lenz, C-J. 

Energieumsetzungen an der Erdoberfläche in gegliedertem Gelände. 

Juni 1996 

Nr. 20: Schwartz, A. 

Numerische Simulationen zur Massenbilanz chemisch reaktiver 
Substanzen im mesoskaligen Bereich. 

November 1996 

Nr. 21: Beheng, K. D. 

Professor Dr. Franz Fiedler zum 60. Geburtstag. 

Januar 1998 

Nr. 22: Niemann, V. 

Numerische Simulation turbulenter Scherströmungen mit einem 
Kaskadenmodell. 

April 1998 

Nr. 23: Koßmann, M. 

Einfluß orographisch induzierter Transportprozesse auf die Struktur 
der atmosphärischen Grenzschicht und die Verteilung von 
Spurengasen. 

April 1998 

Nr. 24: Baldauf, M. 

Die effektive Rauhigkeit über komplexem Gelände - Ein Störungs- 
theoretischer Ansatz. 

Juni 1998 

Nr. 25: Noppel, H. 

Untersuchung des vertikalen Wärmetransports durch die Hangwind- 
zirkulation auf regionaler Skala. 

Dezember 1999 

Nr. 26: Kuntze, K. 

Vertikaler Austausch und chemische Umwandlung von Spurenstoffen 
über topographisch gegliedertem Gelände. 

Oktober 2001 

Nr. 27: l/l//7ms- Grabe, W. 

Vierdimensionale Datenassimilation als Methode zur Kopplung zweier 
verschiedenskaliger meteorologischer Modellsysteme. 

Oktober 2001 



Copyrighted maienal 



Nr. 28: Grabe, F. 

Simulation der Wechselwiricung zwischen Atmosphäre, Vegetation 
und Erdoberfläche bei Verwendung unterschiedlicher Parametri- 
sierungsansätze. 

Januar 2002 

Nr. 29: Riemer, N. 

Numerische Simulationen zur Wirkung des Aerosols auf die tropos- 
phärische Chemie und die Sichtweite. 

Mai 2002 

Nr. 30: Braun, F. J. 

Mesoskalige Modellierung der Bodenhydrologie. 

Dezember 2002 

Nr. 31: Kunz, M. 

Simulation von Starkniederschlägen mit langer Andauer über 
Mittelgebirgen. 

März 2003 

Nr. 32: Bäumer, D. 

Transport und chemische Umwandlung von Luftschadstoffen im 
Nahbereich von Autobahnen - numerische Simulationen. 

Juni 2003 

Nr. 33: Barthlott, C. 

Kohärente Wirbelstrukturen in der atmosphärischen Grenzschicht. 

Juni 2003 

Nr. 34: Wieser, A. 

Messung turbulenter Spurengasflüsse vom Flugzeug aus. 

Januar 2005 

Nr. 35: Blalial<, U. 

Analyse des Extinktionseffektes bei Niederschlagsmessungen mit einem 
C-Band Radar anhand von Simulation und Messung. 

Februar 2005 

Nr. 36: Bertram, I. 

Bestimmung der Wasser- und Eismasse hochreichender konvektiver 
Wolken anhand von Radardaten, Modellergebnissen und konzep- 
tioneller Betrachtungen. 
Mai 2005 

Nr. 37: Schmoecl<el, J. 

Orographischer Einfluss auf die Strömung abgeleitet aus Sturm- 
schäden im Schwarzwald während des Orkans „Lothar". 

Mai 2006 



Copyrighted maierial 



Nr. 38: Schmitt, C. 

Interannual Variability in Antarctic Sea Ice Motion: Interannuelle 
Variabilität antarictischer Meereis-Drift. 

Mai 2006 

Nr. 39: Hasel, M. 

Strukturmerkmale und Modelldarstellung der Konvektion über 
Mittelgebirgen. 

Juli 2006 



Ab Band 40 erscheinen die Wissenschaftlichen Berichte des Instituts für Meteo- 
rologie und Klimaforschung bei KIT Scientific Publishing (ISSN 0179-5619). 

Die Bände sind unter www.ksp.kit.edu als PDF frei verfügbar oder als 
Druckausgabe bestellbar. 

Nr. 40: Lux, R. 

Modellsimulationen zur Strömungsverstärkung von orographischen 
Grundstrukturen bei Sturmsituationen. (2007) 
ISBN 978-3-86644-140-8 

Nr. 41: Straub, W. 

Der Einfluss von Gebirgswellen auf die Initiierung und Entwicklung 
konvektiver Wolken. (2008) 
ISBN 978-3-86644-226-9 

Nr. 42: Meißner, C. 

High-resolution sensitivity studies with the regional climate model 
COSMO-CLM. (2008) 
ISBN 978-3-86644-228-3 

Nr. 43: Höpfner, M. 

Charakterisierung polarer stratosphärischer Wolken mittels hochauflö- 
sender Infrarotspektroskopie. (2008) 
ISBN 978-3-86644-294-8 

Nr. 44: Rings, J. 

Monitoring the water content evolution of dikes. (2009) 
ISBN 978-3-86644-321-1 

Nr. 45: Riemer, M. 

Außertropische Umwandlung tropischer Wirbelstürme: Einfluss auf 
das Strömungsmuster in den mittleren Breiten. (2012) 
ISBN 978-3-86644-766-0 



Copyrighted maienal 



Nr. 4fi: Anwender, D. 

Extratropical Transition in the Ensemble Prediction System of the ECMWF: 
Gase Studies and Experiments. (201 2) 
ISBN 978-3-86644-767-7 

Nr. 4Z: Rinke, R. 

Parametrisierung des Auswaschens von Aerosolpartikeln durch 
Niederschlag. (2012) 
ISBN 978-3-86644-768-4 

Nr. 4& Stanelle, T. 

Wechselwirkungen von Mineralstaubpartikeln mit thermodynamischen 
und dynamischen Prozessen in der Atmosphäre über Westafrika. (2012) 
ISBN 978-3-86644-769-1 

Nr. 49: Peters, T. 

Ableitung einer Beziehung zwischen der Radarreflektivität, der 
Niederschlagsrate und weiteren aus Radardaten abgeleiteten 
Parametern unter Verwendung von Methoden der multivariaten 
Statistik. (2012) 
ISBN 978-3-86644-323-5 

Nr. 50; Khodayar Pardo, S. 

High-resolution analysis of the Initiation of deep convection forced 
by boundary-layer processes. (2012) 
ISBN 978-3-86644-770-7 

Nr. Sil Träum ner, K. 

Einmischprozesse am Oberrand der konvektiven atmosphärischen 
Grenzschicht. (2012) 
ISBN 978-3-86644-771-4 

Nr. 52: Schwendike, J. 

Convection in an African Easterly Wave over West Africa and the 
Eastern Atlantic: A Model Gase Study of Hurricane Helene (2006) 
and its Interaction with the Saharan Air Layer. (201 2) 
ISBN 978-3-86644-772-1 

Nr. 53: Lundgren, K. 

Direct Radiative Effects of Sea Salt on the Regional Scale. (201 2) 
ISBN 978-3-86644-773-8 

Nr. 54: Sasse, R. 

Analyse des regionalen atmosphärischen Wasserhaushalts unter 
Verwendung von GOSMO-Simulationen und GPS-Beobachtungen. (2012) 
ISBN 978-3-86644-774-5 



Copyrighted maUnal 



Nr. 55: Grenzhäuser, J. 

Entwicklung neuartiger Mess- und Auswertungsstrategien für ein 
scannendes Wolkenradar und deren Anwendungsbereiche. (2012) 
ISBN 978-3-86644-775-2 



Copyrighted maienal 



55 



JENS GRENZHÄUSER 

Entwicklung neuartiger IVIess- und Auswertungsstrategien 

für ein scannendes Woll<enradar und deren Anwendungsbereiche 

Die Beobachtung wolkenphysikalischer Prozesse wird erst durch räumlich 
und zeitlich hoch aufgelöste Messungen mit Wolkenradaren möglich. Nach- 
dem vertikal zeigende Geräte weit verbreitet sind, erweitern scannende 
Wolkenradare die Messmöglichkeiten, denn sie erlauben es, einen bestimm- 
ten Bereich der Wolke mehrfach nacheinander zu erfassen. Nur so kann 
die zeitliche Entwicklung innerhalb der Wolke getrennt von der räumlichen 
Variabilität erfasst werden. 

In der vorliegenden Arbeit werden Messungen mit dem scannenden 35.5 GHz 
Wolkenradar MIRA36-S vorgestellt. Der Schwerpunkt liegt dabei auf der Ent- 
wicklung geeigneter Scan- und Auswerte verfahren, die von der neuen tech- 
nischen Möglichkeit Gebrauch machen. Bei den meisten hier vorgestellten 
Messungen handelt es sich um eine schnelle Abfolge von RHI-Scans, die in 
Richtung des vorherrschenden Windes durchgeführt wurden. Die Auswer- 
teverfahren sind so konzipiert, dass sie auch automatisiert werden können, 
also fijr den operationeilen Einsatz geeignet sind. 




